Source author record

Ilya Mandel

Ilya Mandel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE gr-qc astro-ph.SR astro-ph.IM astro-ph.CO physics.data-an astro-ph.GA astro-ph astro-ph.EP physics.ins-det

Catalog footprint

What is connected

77works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2025arXiv

Converged simulations of the nozzle shock in tidal disruption events

When debris from a star that experienced a tidal disruption events (TDE) after passing too close to a massive black hole returns to pericenter on the second passage, it is compressed, leading to the formation of nozzle shocks (in the orbital plane) and pancake shocks (perpendicular to the orbital plane). Resolving these shocks is a long-standing problem in the hydrodynamic simulations of parabolic TDEs. Excessive numerical energy dissipation or heating unrealistically expands the stream. In this Letter, we apply adaptive particle refinement to our 3D general relativistic smoothed particle simulations to locally increase the resolution near the pericenter. We achieve resolutions equivalent to $6.55\times10^{11}$ particles, allowing us to converge on the true energy dissipation. We conclude that only $4\times10^{-5}$ of the orbital energy is dissipated in nozzle shocks for a Sun-like star tidally disrupted by a $10^6$ solar-mass black hole, therefore the nozzle shocks are unlikely to be important in the evolution of TDEs.

preprint2022arXiv

Building better spin models for merging binary black holes: Evidence for non-spinning and rapidly spinning nearly aligned sub-populations

Recent work paints a conflicting portrait of the distribution of black hole spins in merging binaries measured with gravitational waves. Some analyses find that a significant fraction of merging binaries contain at least one black hole with a spin tilt $>90^\circ$ with respect to the orbital angular momentum vector, which has been interpreted as a signature for dynamical assembly. Other analyses find the data are consistent with a bimodal population in which some binaries contain black holes with negligible spin while the rest contain black holes with spin vectors preferentially aligned with the orbital angular momentum vector. In this work, we scrutinize models for the distribution of black hole spins to pinpoint possible failure modes in which the model yields a faulty conclusion. We reanalyze data from the second LIGO--Virgo gravitational-wave transient catalog (GWTC-2) using a revised spin model, which allows for a sub-population of black holes with negligible spins. In agreement with recent results by Roulet et al., we show that the GWTC-2 detections are consistent with two distinct sub-populations. We estimate that $29-75\%$ (90\% credible interval) of merging binaries contain black holes with negligible spin $χ\approx 0$. The remaining binaries are part of a second sub-population in which the spin vectors are preferentially (but not exactly) aligned to the orbital angular momentum. The black holes in this second sub-population are characterized by spins of $χ\sim0.45$. We suggest that the inferred spin distribution is consistent with the hypothesis that all merging binaries form via the field formation scenario.

preprint2022arXiv

Common envelopes in massive stars II: The distinct roles of hydrogen and helium recombination

The role of recombination during a common-envelope event has been long debated. Many studies have argued that much of hydrogen recombination energy, which is radiated in relatively cool and optically-thin layers, might not thermalise in the envelope. On the other hand, helium recombination contains 30% of the total recombination energy, and occurs much deeper in the stellar envelope. We investigate the distinct roles played by hydrogen and helium recombination in a common-envelope interaction experienced by a 12 solar mass red supergiant donor. We perform adiabatic, 3D hydrodynamical simulations that (i) include hydrogen, helium, and molecular hydrogen recombination, (ii) include hydrogen and helium recombination, (iii) include only helium recombination, and (iv) do not include recombination energy. By comparing these simulations, we find that the addition of helium recombination energy alone ejects 30% more envelope mass, and leads to a 16% larger post-plunge-in separation. Under the adiabatic assumption, adding hydrogen recombination energy increases the amount of ejected mass by a further 40%, possibly unbinding the entire envelope, but does not affect the post-plunge separation. Most of the ejecta becomes unbound at relatively high (>70%) degrees of hydrogen ionisation, where the hydrogen recombination energy is likely to expand the envelope instead of being radiated away.

preprint2022arXiv

Common envelopes in massive stars: Towards the role of radiation pressure and recombination energy in ejecting red supergiant envelopes

We perform 3D hydrodynamical simulations of a common-envelope event involving a 12 solar mass red supergiant donor. Massive stars are expected to be qualitatively different from low-mass stars as their envelopes have significant support from radiation pressure, which increases both the final separation and amount of mass ejected through the common-envelope interaction. We perform adiabatic simulations that include radiation energy through the equation of state, which results in ejecting 60 per cent more mass (up to two thirds of the total envelope mass becoming unbound, or more) and yield a 10 per cent larger final separation compared to simulations that assume an ideal gas. When also including recombination energy, we find that at least three quarters of the envelope, and possibly the entire envelope, may be unbound. The final separation further increases by almost 20 per cent. The additional amount of ejected material is mainly due to energy injected from helium recombination. Hydrogen recombination plays a comparatively small role, as it mainly occurs in gas that has already become unbound. We conclude that the internal energy of the envelope can be a significant energy source for ejecting the common envelope, but ultimately radiation transport and convection need to be included.

preprint2022arXiv

Extreme eccentricities of triple systems: Analytic results

Triple stars and compact objects are ubiquitously observed in nature. Their long-term evolution is complex; in particular, the von-Zeipel-Lidov-Kozai (ZLK) mechanism can potentially lead to highly eccentric encounters of the inner binary. Such encounters can lead to a plethora of interacting binary phenomena, as well as stellar and compact-object mergers. Here we find explicit analytical formulae for the maximal eccentricity, $e_{\rm max}$, of the inner binary undergoing ZLK oscillations, where both the test particle limit (parametrised by the inner-to-outer angular momentum ratio $η$) and the double-averaging approximation (parametrised by the period ratio, $ε_{\rm SA}$) are relaxed, for circular outer orbits. We recover known results in both limiting cases (either $η$ or $ε_{\rm SA} \to 0$) and verify the validity of our model using numerical simulations. We test our results with two accurate numerical N-body codes, $\texttt{Rebound}$ for Newtonian dynamics and $\texttt{Tsunami}$ for general-relativistic (GR) dynamics, and find excellent correspondence. We discuss the implications of our results for stellar triples and both stellar and supermassive triple black hole mergers.

preprint2022arXiv

Merging stellar-mass binary black holes

The LIGO and Virgo detectors have directly observed gravitational waves from mergers of pairs of stellar-mass black holes, along with a smaller number of mergers involving neutron stars. These observations raise the hope that compact object mergers could be used as a probe of stellar and binary evolution, and perhaps of stellar dynamics. This colloquium-style article summarises the existing observations, describes theoretical predictions for formation channels of merging stellar-mass black-hole binaries along with their rates and observable properties, and presents some prospects for gravitational-wave astronomy.

preprint2022arXiv

Rates of Compact Object Coalescences

Gravitational-wave detections are enabling measurements of the rate of coalescences of binaries composed of two compact objects -- neutron stars and/or black holes. The coalescence rate of binaries containing neutron stars is further constrained by electromagnetic observations, including Galactic radio binary pulsars and short gamma-ray bursts. Meanwhile, increasingly sophisticated models of compact objects merging through a variety of evolutionary channels produce a range of theoretically predicted rates. Rapid improvements in instrument sensitivity, along with plans for new and improved surveys, make this an opportune time to summarise the existing observational and theoretical knowledge of compact-binary coalescence rates.

preprint2022arXiv

Stellar response after stripping as a model for common-envelope outcomes

Binary neutron stars have been observed as millisecond pulsars, gravitational-wave sources, and as the progenitors of short gamma-ray bursts and kilonovae. Massive stellar binaries that evolve into merging double neutron stars are believed to experience a common-envelope episode. During this episode, the envelope of a giant star engulfs the whole binary. The energy transferred from the orbit to the envelope by drag forces or from other energy sources can eject the envelope from the binary system, leading to a stripped short-period binary. In this paper, we use one-dimensional single stellar evolution to explore the final stages of the common-envelope phase in progenitors of neutron star binaries. We consider an instantaneously stripped donor star as a proxy for the common-envelope phase and study the star's subsequent radial evolution. We determine a range of stripping boundaries which allow the star to avoid significant rapid re-expansion and which thus represent plausible boundaries for the termination of the common-envelope episode. We find that these boundaries lie above the maximum compression point, a commonly used location of the core/envelope boundary. We conclude that stars may retain fractions of a solar mass of hydrogen-rich material even after the common-envelope episode. If we consider orbital energy as the only energy source available, all of our models would overfill their Roche lobe after ejecting the envelope, whose binding energy includes gravitational, thermal, radiation, and recombination energy terms.

preprint2021arXiv

Cygnus X-1 contains a 21-solar mass black hole -- implications for massive star winds

The evolution of massive stars is influenced by the mass lost to stellar winds over their lifetimes. These winds limit the masses of the stellar remnants (such as black holes) that the stars ultimately produce. We use radio astrometry to refine the distance to the black hole X-ray binary Cygnus X-1, which we find to be $2.22^{+0.18}_{-0.17}$ kiloparsecs. When combined with previous optical data, this implies a black hole mass of $21.2\pm2.2$ solar masses, higher than previous measurements. The formation of such a high-mass black hole in a high-metallicity system constrains wind mass loss from massive stars.

preprint2021arXiv

Impact of Massive Binary Star and Cosmic Evolution on Gravitational Wave Observations II: Double Compact Object Rates and Properties

Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, $Z$, and redshift $z, \mathcal{S}(Z,z)$. Considering these uncertainties we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal{S}(Z,z)$ variations can impact the predicted intrinsic and detectable merger rates by factors $10^2$-$10^4$. We find that BHBH rates are dominantly impacted by $\mathcal{S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal{S}(Z,z)$ for BHBH, BHNS and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95% of BHBH detections to contain a BH $\gtrsim 8\,\rm{M}_{\odot}$ and have mass ratios $\lesssim 4$. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties.

preprint2021arXiv

Rapid stellar and binary population synthesis with COMPAS

Compact Object Mergers: Population Astrophysics and Statistics (COMPAS; https://compas.science) is a public rapid binary population synthesis code. COMPAS generates populations of isolated stellar binaries under a set of parametrized assumptions in order to allow comparisons against observational data sets, such as those coming from gravitational-wave observations of merging compact remnants. It includes a number of tools for population processing in addition to the core binary evolution components. COMPAS is publicly available via the github repository https://github.com/TeamCOMPAS/COMPAS/, and is designed to allow for flexible modifications as evolutionary models improve. This paper describes the methodology and implementation of COMPAS. It is a living document which will be updated as new features are added to COMPAS; the current document describes COMPAS v02.21.00.

preprint2021arXiv

Wind mass-loss rates of stripped stars inferred from Cygnus X-1

Recent observations of the high-mass X-ray binary Cygnus X-1 have shown that both the companion star (41 solar masses) and the black hole (21 solar masses) are more massive than previously estimated. Furthermore, the black hole appears to be nearly maximally spinning. Here we present a possible formation channel for the Cygnus X-1 system that matches the observed system properties. In this formation channel, we find that the orbital parameters of Cygnus X-1, combined with the observed metallicity of the companion, imply a significant reduction in mass loss through winds relative to commonly used prescriptions for stripped stars.

preprint2020arXiv

Age dissection of the Milky Way discs: red giants in the Kepler field

[Abridged] Ensemble studies of red-giant stars with exquisite asteroseismic, spectroscopic, and astrometric constraints offer a novel opportunity to recast and address long-standing questions concerning the evolution of stars and of the Galaxy. Here, we infer masses and ages for nearly 5400 giants with available Kepler light curves and APOGEE spectra, and discuss some of the systematics that may affect the accuracy of the inferred stellar properties. First, we look at age-chemical-abundances relations. We find a dearth of young, metal-rich stars, and the existence of a significant population of old (8-9 Gyr), low-[$α$/Fe], super-solar metallicity stars, reminiscent of the age and metallicity of the well-studied open cluster NGC6791. The age-chemo-kinematic properties of these stars indicate that efficient radial migration happens in the thin disk. We find that ages and masses of the nearly 400 $α$-element-rich red-giant-branch (RGB) stars in our sample are compatible with those of an old (~11 Gyr), nearly coeval, chemical-thick disk population. Using a statistical model, we show that 95% of the population was born within ~1.5 Gyr. Moreover, we find a difference in the vertical velocity dispersion between low- and high-[$α$/Fe] populations, confirming their different chemo-dynamical histories. We then exploit the almost coeval $α$-rich population to gain insight into processes that may have altered the mass of a star along its evolution, which are key to improve the mapping of the observed stellar mass to age. We find evidence for a mean integrated RGB mass loss <$Δ$M>= 0.10 $\pm$ 0.02 Msun and that the occurrence of massive (M $\gtrsim$ 1.1 Msun) $α$-rich stars is of the order of 5% on the RGB, and significantly higher in the RC, supporting the scenario in which most of these stars had undergone interaction with a companion.

preprint2020arXiv

An alternative interpretation of GW190412 as a binary black hole merger with a rapidly spinning secondary

The LIGO-Virgo collaboration recently reported the properties of GW190412, a binary black hole merger with unequal component masses (mass ratio $0.25^{+0.06}_{-0.04}$ when using the EOBNR PHM approximant) and a non-vanishing effective spin aligned with the orbital angular momentum. They used uninformative priors to infer that the more massive black hole had a dimensionless spin magnitude between 0.17 and 0.59 at 90\% confidence. We argue that, within the context of isolated binary evolution, it is more natural to assume a priori that the first-born, more massive black hole has a negligible spin, while the spin of the less massive black hole is preferentially aligned with the orbital angular momentum if it is spun up by tides. Under this astrophysically motivated prior, we conclude that the lower mass black hole had a dimensionless spin component between 0.64 and 0.99 along the orbital angular momentum.

preprint2020arXiv

Be X-ray binaries in the SMC as indicators of mass transfer efficiency

Be X-ray binaries (BeXRBs) consist of rapidly rotating Be stars with neutron star companions accreting from the circumstellar emission disk. We compare the observed population of BeXRBs in the Small Magellanic Cloud with simulated populations of BeXRB-like systems produced with the COMPAS population synthesis code. We focus on the apparently higher minimal mass of Be stars in BeXRBs than in the Be population at large. Assuming that BeXRBs experienced only dynamically stable mass transfer, their mass distribution suggests that at least 30% of the mass donated by the progenitor of the neutron star is typically accreted by the B-star companion. We expect these results to affect predictions for the population of double compact object mergers. A convolution of the simulated BeXRB population with the star formation history of the Small Magellanic Cloud shows that the excess of BeXRBs is most likely explained by this galaxy's burst of star formation around 20--40 Myr ago.

preprint2020arXiv

Constraining the masses of microlensing black holes and the mass gap with Gaia DR2

Context: Gravitational microlensing is sensitive to compact-object lenses in the Milky Way, including white dwarfs, neutron stars or black holes, and could potentially probe a wide range of stellar remnant masses. However, the mass of the lens can be determined only in very limited cases, due to missing information on both source and lens distances and their proper motions. Aims: We aim at improving the mass estimates in the annual parallax microlensing events found in the 8 years of OGLE-III observations towards the Galactic Bulge (Wyrzykowski et al. 2016) with the use of Gaia Data Release 2 (DR2). Methods: We use Gaia DR2 data on distances and proper motions of non-blended sources and recompute the masses of lenses in parallax events. We also identify new events in that sample which are likely to have dark lens; the total number of such events is now 18. Results: The derived distribution of masses of dark lenses is consistent with a continuous distribution of stellar remnant masses. A mass gap between neutron-star and black-hole masses in the range between 2 and 5 solar masses is not favoured by our data, unless black holes receive natal-kicks above 20-80 km/s. We present 8 candidates for objects with masses within the putative mass gap, including a spectacular multi-peak parallax event with mass of $2.4^{+1.9}_{-1.3}\ M_\odot$ located just at 600 pc. The absence of an observational mass gap between neutron stars and black holes, or, conversely, the evidence for black hole natal kicks if a mass gap is assumed, can inform future supernova modelling efforts.

preprint2020arXiv

Detecting Double Neutron Stars with LISA

We estimate the properties of the double neutron star (DNS) population that will be observable by the planned space-based interferometer LISA. By following the gravitational radiation driven evolution of DNSs generated from rapid population synthesis of massive binary stars, we estimate that around 35 DNSs will accumulate a signal-to-noise ratio above 8 over a four-year LISA mission. The observed population mainly comprises Galactic DNSs (94 per cent), but detections in the LMC (5 per cent) and SMC (1 per cent) may also be expected. The median orbital frequency of detected DNSs is expected to be 0.8 mHz, and many of them will be eccentric (median eccentricity of $0.11$). The orbital properties will provide insights into DNS progenitors and formation channels. LISA is expected to localise these DNSs to a typical angular resolution of $2^\circ$, with best-constrained sources localised to a few arcminutes. These localisations may allow neutron star natal kick magnitudes to be constrained through the Galactic distribution of DNSs, and make it possible to follow up the sources with radio pulsar searches. However, LISA is also expected to resolve $\sim 10^4$ Galactic double white dwarfs, many of which may have binary parameters that resemble DNSs; we discuss how the combined measurement of binary eccentricity, chirp mass, and sky location may aid the identification of a DNS. We expect the best-constrained DNSs to have eccentricities known to a few parts in a thousand, chirp masses measured to better than 1 per cent fractional uncertainty, and sky localisation at the level of a few arcminutes.

preprint2020arXiv

Luminous Red Novae: population models and future prospects

A class of optical transients known as Luminous Red Novae (LRNe) have recently been associated with mass ejections from binary stars undergoing common-envelope evolution. We use the population synthesis code COMPAS to explore the impact of a range of assumptions about the physics of common-envelope evolution on the properties of LRNe. In particular, we investigate the influence of various models for the energetics of LRNe on the expected event rate and light curve characteristics, and compare with the existing sample. We find that the Galactic rate of LRNe is $\sim 0.2$ yr$^{-1}$, in agreement with the observed rate. In our models, the luminosity function of Galactic LRNe covers multiple decades in luminosity and is dominated by signals from stellar mergers, consistent with observational constraints from iPTF and the Galactic sample of LRNe. We discuss how observations of the brightest LRNe may provide indirect evidence for the existence of massive ($> 40$ M$_\odot$) red supergiants. Such LRNe could be markers along the evolutionary pathway leading to the formation of double compact objects. We make predictions for the population of LRNe observable in future transient surveys with the Large Synoptic Survey Telescope and the Zwicky Transient Facility. In all plausible circumstances, we predict a selection-limited observable population dominated by bright, long-duration events caused by common envelope ejections. We show that the Large Synoptic Survey Telescope will observe $20$--$750$ LRNe per year, quickly constraining the luminosity function of LRNe and probing the physics of common-envelope events.

preprint2020arXiv

Merger rates in primordial black hole clusters without initial binaries

Primordial black holes formed through the collapse of cosmological density fluctuations have been hypothesised as contributors to the dark matter content of the Universe. At the same time, their mergers could contribute to the recently observed population of gravitational-wave sources. We investigate the scenario in which primordial black holes form binaries at late times in the Universe. Specifically, we re-examine the mergers of primordial black holes in small clusters of ~30 objects in the absence of initial binaries. Binaries form dynamically through Newtonian gravitational interactions. These binaries act as heat sources for the cluster, increasing the cluster's velocity dispersion, which inhibits direct mergers through gravitational-wave two-body captures. Meanwhile, three-body encounters of tight binaries are too rare to tighten binaries sufficiently to allow them to merge through gravitational-wave emission. We conclude that in the absence of initial binaries, merger rates of primordial black holes in the Bird et al. (2016) initial cluster configuration are at least an order of magnitude lower than previously suggested, which makes gravitational-wave detections of such sources improbable.

preprint2020arXiv

The origin of spin in binary black holes: Predicting the distributions of the main observables of Advanced LIGO

We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of the main observables of Advanced LIGO. We used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. By assuming efficient angular momentum transport, we are able to present a model that is capable of simultaneously predicting the three main gravitational-wave observables: the effective inspiral spin parameter $χ_{eff}$, the chirp mass $M_{chirp}$, and the cosmological redshift of merger $z_{merger}$. We find an excellent agreement between our model and the ten events from the first two advanced detector observing runs. We make predictions for the third observing run O3 and for Advanced LIGO design sensitivity. We expect around 80% of events with $χ_{eff} < 0.1$, while the remaining 20% of events with $χ_{eff} \ge 0.1$ are split into ~10% with $M_{chirp} < 15$ M$_\odot$ and ~10% with $M_{chirp} \ge 15$ M$_\odot$. In conclusion, the favorable comparison of the existing LIGO/Virgo observations with our model predictions gives support to the idea that the majority, if not all of the observed mergers, originate from the evolution of isolated binaries. The first-born black hole has negligible spin because it lost its envelope after it expanded to become a giant star, while the spin of the second-born black hole is determined by the tidal spin up of its naked helium star progenitor by the first-born black hole companion after the binary finished the common-envelope phase.

preprint2020arXiv

The spectral evolution of AT 2018dyb and the presence of metal lines in tidal disruption events

We present light curves and spectra of the tidal disruption event (TDE) ASASSN-18pg / AT 2018dyb spanning a period of one year. The event shows a plethora of strong emission lines, including the Balmer series, He II, He I and metal lines of O III $λ$3760 and N III $λλ$ 4100, 4640 (blended with He II). The latter lines are consistent with originating from the Bowen fluorescence mechanism. By analyzing literature spectra of past events, we conclude that these lines are common in TDEs. The spectral diversity of optical TDEs is thus larger than previously thought and includes N-rich events besides H- and He-rich events. We study how the spectral lines evolve with time, by means of their width, relative strength, and velocity offsets. The velocity width of the lines starts at $\sim$ 13000 km s$^{-1}$ and decreases with time. The ratio of He II to N III increases with time. The same is true for ASASSN-14li, which has a very similar spectrum to AT 2018dyb but its lines are narrower by a factor of $>$2. We estimate a black hole mass of $M_{\rm BH}$ = $3.3^{+5.0}_{-2.0}\times 10^6$ $M_{\odot}$ by using the $M$-$σ$ relation. This is consistent with the black hole mass derived using the MOSFiT transient fitting code. The detection of strong Bowen lines in the optical spectrum is an indirect proof for extreme ultraviolet and (reprocessed) X-ray radiation and favors an accretion origin for the TDE optical luminosity. A model where photons escape after multiple scatterings through a super-Eddington thick disk and its optically thick wind, viewed at an angle close to the disk plane, is consistent with the observations.

preprint2019arXiv

Astrophysical science metrics for next-generation gravitational-wave detectors

The second generation of gravitational-wave detectors are being built and tuned all over the world. The detection of signals from binary black holes is beginning to fulfill the promise of gravitational-wave astronomy. In this work, we examine several possible configurations for third-generation laser interferometers in existing km-scale facilities. We propose a set of astrophysically motivated metrics to evaluate detector performance. We measure the impact of detector design choices against these metrics, providing a quantitative cost-benefit analyses of the resulting scientific payoffs.

preprint2019arXiv

Common Envelope Wind Tunnel: The Effects of Binary Mass Ratio and Implications for the Accretion-Driven Growth of LIGO Binary Black Holes

We present three-dimensional local hydrodynamic simulations of flows around objects embedded within stellar envelopes using a "wind tunnel" formalism. Our simulations model the common envelope dynamical inspiral phase in binary star systems in terms of dimensionless flow characteristics. We present suites of simulations that study the effects of varying the binary mass ratio, stellar structure, equation of state, relative Mach number of the object's motion through the gas, and density gradients across the gravitational focusing scale. For each model, we measure coefficients of accretion and drag experienced by the embedded object. These coefficients regulate the coupled evolution of the object's masses and orbital tightening during the dynamical inspiral phase of the common envelope. We extrapolate our simulation results to accreting black holes with masses comparable to that of the population of LIGO black holes. We demonstrate that the mass and spin accrued by these black holes per unit orbital tightening are directly related to the ratio of accretion to drag coefficients. We thus infer that the mass and dimensionless spin of initially non-rotating black holes change by of order $1\%$ and 0.05, respectively, in a typical example scenario. Our prediction that the masses and spins of black holes remain largely unmodified by a common envelope phase aids in the interpretation of the properties of the growing observed population of merging binary black holes. Even if these black holes passed through a common envelope phase during their assembly, features of mass and spin imparted by previous evolutionary epochs should be preserved.

preprint2019arXiv

Eclipses of continuous gravitational waves as a probe of stellar structure

Although gravitational waves only interact weakly with matter, their propagation is affected by a gravitational potential. If a gravitational wave source is eclipsed by a star, measuring these perturbations provides a way to directly measure the distribution of mass throughout the stellar interior. We compute the expected Shapiro time delay, amplification and deflection during an eclipse, and show how this can be used to infer the mass distribution of the eclipsing body. We identify continuous gravitational waves from neutron stars as the best candidates to detect this effect. When the Sun eclipses a far-away source, depending on the depth of the eclipse the time-delay can change by up to $\sim 0.034$ ms, the gravitational-wave strain amplitude can increase by $\sim 4$%, and the apparent position of the source in the sky can vary by $4''$. Accreting neutron stars with Roche-lobe filling companion stars have a high probability of exhibiting eclipses, producing similar time delays but undetectable changes in amplitude and sky location. Even for the most rapidly rotating neutron stars, this time delay only corresponds to a few percent of the phase of the gravitational wave, making it an extremely challenging measurement. However, if sources of continuous gravitational waves exist just below the limit of detection of current observatories, next-generation instruments will be able to observe them with enough precision to measure the signal of an eclipsing star. Detecting this effect would provide a new direct probe to the interior of stars, complementing asteroseismology and the detection of solar neutrinos.

preprint2019arXiv

Explosions Driven by the Coalescence of a Compact Object with the Core of a Massive-Star Companion Inside a Common Envelope: Circumstellar Properties, Light Curves, and Population Statistics

We model explosions driven by the coalescence of a black hole or neutron star with the core of its massive-star companion. Upon entering a common envelope phase, a compact object may spiral all the way to the core. The concurrent release of energy is likely to be deposited into the surrounding common envelope, powering a merger-driven explosion. We use hydrodynamic models of binary coalescence to model the common envelope density distribution at the time of coalescence. We find toroidal profiles of material, concentrated in the binary's equatorial plane and extending to many times the massive star's original radius. We use the spherically-averaged properties of this circumstellar material (CSM) to estimate the emergent light curves that result from the interaction between the blast wave and the CSM. We find that typical merger-driven explosions are brightened by up to three magnitudes by CSM interaction. From population synthesis models we discover that the brightest merger-driven explosions, $M_V \sim -18$ to $-19$, are those involving black holes because they have the most massive and extended CSM. Black hole coalescence events are also common; they represent about 50% of all merger-driven explosions and approximately 0.3% of the core-collapse rate. Merger-driven explosions offer a window into the highly-uncertain physics of common envelope interactions in binary systems by probing the properties of systems that merge rather than eject their envelopes.

preprint2018arXiv

Comment on "An excess of massive stars in the local 30 Doradus starburst"

Schneider et al. (Science, 2018) used an ad hoc statistical method in their calculation of the stellar initial mass function. Adopting an improved approach, we reanalyse their data and determine a power law exponent of $2.05_{-0.14}^{+0.13}$. Alternative assumptions regarding data set completeness and the star formation history model can shift the inferred exponent to $2.11_{-0.19}^{+0.17}$ and $2.15_{-0.13}^{+0.13}$, respectively.

preprint2016arXiv

$N-$body dynamics of Intermediate mass-ratio inspirals in globular clusters

The intermediate mass-ratio inspiral of a stellar compact remnant into an intermediate mass black hole (IMBH) can produce a gravitational wave (GW) signal that is potentially detectable by current ground-based GW detectors (e.g., Advanced LIGO) as well as by planned space-based interferometers (e.g., eLISA). Here, we present results from a direct integration of the post-Newtonian $N$-body equations of motion describing stellar clusters containing an IMBH and a population of stellar-mass black holes (BHs) and solar mass stars. We take particular care to simulate the dynamics closest to the IMBH, including post-Newtonian effects up to order $2.5$. Our simulations show that the IMBH readily forms a binary with a BH companion. This binary is gradually hardened by transient 3-body or 4-body encounters, leading to frequent substitutions of the BH companion, while the binary's eccentricity experiences large amplitude oscillations due to the Lidov-Kozai resonance. We also demonstrate suppression of these resonances by the relativistic precession of the binary orbit. We find an intermediate mass-ratio inspiral in one of the 12 cluster models we evolved for $\sim 100$ Myr. This cluster hosts a $100 M_\odot$ IMBH embedded in a population of 32 $10M_\odot$ BH and 32,000 $1M_\odot$ stars. At the end of the simulation, after $\sim 100$ Myr of evolution, the IMBH merges with a BH companion. The IMBH--BH binary inspiral starts in the eLISA frequency window ($\gtrsim 1\rm mHz$) when the binary reaches an eccentricity $1-e\simeq 10^{-3}$. After $\simeq 10^5$ years the binary moves into the LIGO frequency band with a negligible eccentricity. We comment on the implications for GW searches, with a possible detection within the next decade.

preprint2016arXiv

Detectability of Gravitational Waves from High-Redshift Binaries

Recent non-detection of gravitational-wave backgrounds from pulsar timing arrays casts further uncertainty on the evolution of supermassive black hole binaries. We study the capabilities of current gravitational-wave observatories to detect individual binaries and demonstrate that, contrary to conventional wisdom, some are in principle detectable throughout the Universe. In particular, a binary with rest-frame mass $\gtrsim10^{10}\,M_\odot$ can be detected by current timing arrays at arbitrarily high redshifts. The same claim will apply for less massive binaries with more sensitive future arrays. As a consequence, future searches for nanohertz gravitational waves could be expanded to target evolving high-redshift binaries. We calculate the maximum distance at which binaries can be observed with pulsar timing arrays and other detectors, properly accounting for redshift and using realistic binary waveforms.

preprint2016arXiv

Dynamic temperature selection for parallel-tempering in Markov chain Monte Carlo simulations

Modern problems in astronomical Bayesian inference require efficient methods for sampling from complex, high-dimensional, often multi-modal probability distributions. Most popular methods, such as Markov chain Monte Carlo sampling, perform poorly on strongly multi-modal probability distributions, rarely jumping between modes or settling on just one mode without finding others. Parallel tempering addresses this problem by sampling simultaneously with separate Markov chains from tempered versions of the target distribution with reduced contrast levels. Gaps between modes can be traversed at higher temperatures, while individual modes can be efficiently explored at lower temperatures. In this paper, we investigate how one might choose the ladder of temperatures to achieve more efficient sampling, as measured by the autocorrelation time of the sampler. In particular, we present a simple, easily-implemented algorithm for dynamically adapting the temperature configuration of a sampler while sampling. This algorithm dynamically adjusts the temperature spacing to achieve a uniform rate of exchanges between chains at neighbouring temperatures. We compare the algorithm to conventional geometric temperature configurations on a number of test distributions and on an astrophysical inference problem, reporting efficiency gains by a factor of 1.2-2.5 over a well-chosen geometric temperature configuration and by a factor of 1.5-5 over a poorly chosen configuration. On all of these problems a sampler using the dynamical adaptations to achieve uniform acceptance ratios between neighbouring chains outperforms one that does not.

preprint2016arXiv

Kernel regression estimates of time delays between gravitationally lensed fluxes

Strongly lensed variable quasars can serve as precise cosmological probes, provided that time delays between the image fluxes can be accurately measured. A number of methods have been proposed to address this problem. In this paper, we explore in detail a new approach based on kernel regression estimates, which is able to estimate a single time delay given several datasets for the same quasar. We develop realistic artificial data sets in order to carry out controlled experiments to test of performance of this new approach. We also test our method on real data from strongly lensed quasar Q0957+561 and compare our estimates against existing results.

preprint2016arXiv

Limits of Astrophysics with Gravitational-Wave Backgrounds

The recent Advanced LIGO detection of gravitational waves from the binary black hole GW150914 suggests there exists a large population of merging binary black holes in the Universe. Although most are too distant to be individually resolved by advanced detectors, the superposition of gravitational waves from many unresolvable binaries is expected to create an astrophysical stochastic background. Recent results from the LIGO and Virgo collaborations show that this astrophysical background is within reach of Advanced LIGO. In principle, the binary black hole background encodes interesting astrophysical properties, such as the mass distribution and redshift distribution of distant binaries. However, we show that this information will be difficult to extract with the current configuration of advanced detectors (and using current data analysis tools). Additionally, the binary black hole background also constitutes a foreground that limits the ability of advanced detectors to observe other interesting stochastic background signals, for example from cosmic strings or phase transitions in the early Universe. We quantify this effect.

preprint2016arXiv

Merging binary black holes formed through chemically homogeneous evolution in short-period stellar binaries

We explore a newly proposed channel to create binary black holes of stellar origin. This scenario applies to massive, tight binaries where mixing induced by rotation and tides transports the products of hydrogen burning throughout the stellar envelopes. This slowly enriches the entire star with helium, preventing the build-up of an internal chemical gradient. The stars remain compact as they evolve nearly chemically homogeneously, eventually forming two black holes, which, we estimate, typically merge 4--11 Gyr after formation. Like other proposed channels, this evolutionary pathway suffers from significant theoretical uncertainties, but could be constrained in the near future by data from advanced ground-based gravitational-wave detectors. We perform Monte Carlo simulations of the expected merger rate over cosmic time to explore the implications and uncertainties. Our default model for this channel yields a local binary black hole merger rate of about $10$ Gpc$^{-3}$ yr$^{-1}$ at redshift $z=0$, peaking at twice this rate at $z=0.5$. This means that this channel is competitive, in terms of expected rates, with the conventional formation scenarios that involve a common-envelope phase during isolated binary evolution or dynamical interaction in a dense cluster. The events from this channel may be distinguished by the preference for nearly equal-mass components and high masses, with typical total masses between 50 and 110 $\textrm{M}_\odot$. Unlike the conventional isolated binary evolution scenario that involves shrinkage of the orbit during a common-envelope phase, short time delays are unlikely for this channel, implying that we do not expect mergers at high redshift.

preprint2016arXiv

Parameter estimation on gravitational waves from neutron-star binaries with spinning components

Inspiraling binary neutron stars are expected to be one of the most significant sources of gravitational-wave signals for the new generation of advanced ground-based detectors. We investigate how well we could hope to measure properties of these binaries using the Advanced LIGO detectors, which began operation in September 2015. We study an astrophysically motivated population of sources (binary components with masses $1.2~\mathrm{M}_\odot$--$1.6~\mathrm{M}_\odot$ and spins of less than $0.05$) using the full LIGO analysis pipeline. While this simulated population covers the observed range of potential binary neutron-star sources, we do not exclude the possibility of sources with parameters outside these ranges; given the existing uncertainty in distributions of mass and spin, it is critical that analyses account for the full range of possible mass and spin configurations. We find that conservative prior assumptions on neutron-star mass and spin lead to average fractional uncertainties in component masses of $\sim 16\%$, with little constraint on spins (the median $90\%$ upper limit on the spin of the more massive component is $\sim 0.7$). Stronger prior constraints on neutron-star spins can further constrain mass estimates, but only marginally. However, we find that the sky position and luminosity distance for these sources are not influenced by the inclusion of spin; therefore, if LIGO detects a low-spin population of BNS sources, less computationally expensive results calculated neglecting spin will be sufficient for guiding electromagnetic follow-up.

preprint2015arXiv

An Efficient Interpolation Technique for Jump Proposals in Reversible-Jump Markov Chain Monte Carlo Calculations

Selection among alternative theoretical models given an observed data set is an important challenge in many areas of physics and astronomy. Reversible-jump Markov chain Monte Carlo (RJMCMC) is an extremely powerful technique for performing Bayesian model selection, but it suffers from a fundamental difficulty: it requires jumps between model parameter spaces, but cannot efficiently explore both parameter spaces at once. Thus, a naive jump between parameter spaces is unlikely to be accepted in the MCMC algorithm and convergence is correspondingly slow. Here we demonstrate an interpolation technique that uses samples from single-model MCMCs to propose inter-model jumps from an approximation to the single-model posterior of the target parameter space. The interpolation technique, based on a kD-tree data structure, is adaptive and efficient in modest dimensionality. We show that our technique leads to improved convergence over naive jumps in an RJMCMC, and compare it to other proposals in the literature to improve the convergence of RJMCMCs. We also demonstrate the use of the same interpolation technique as a way to construct efficient "global" proposal distributions for single-model MCMCs without prior knowledge of the structure of the posterior distribution, and discuss improvements that permit the method to be used in higher-dimensional spaces efficiently.

preprint2015arXiv

Distinguishing types of compact-object binaries using the gravitational-wave signatures of their mergers

We analyze the distinguishability of populations of coalescing binary neutron stars, neutron-star black-hole binaries, and binary black holes, whose gravitational-wave signatures are expected to be observed by the advanced network of ground-based interferometers LIGO and Virgo. We consider population-synthesis predictions for plausible merging binary distributions in mass space, along with measurement accuracy estimates from the main gravitational-wave parameter-estimation pipeline. We find that for our model compact-object binary mass distribution, we can always distinguish binary neutron stars and black-hole--neutron-star binaries, but not necessarily black-hole--neutron-star binaries and binary black holes; however, with a few tens of detections, we can accurately identify the three subpopulations and measure their respective rates.

preprint2015arXiv

Double tidal disruptions in galactic nuclei

A star on a nearly radial trajectory approaching a massive black hole (MBH) gets tidally disrupted if it comes sufficiently close to the MBH. Here we explore what happens to binary stars whose centers of mass approach the MBH on nearly radial orbits. The interaction with the MBH often leads to both stars being disrupted in sequence. We argue that such events could produce light curves that are substantially different from those of the single disruptions, with possible features such as two local maxima. Tidal forces from the MBH can also lead the binary components to collide; these merger products can form highly magnetized stars, whose subsequent tidal disruption may enable prompt jet formation.

preprint2015arXiv

Efficient method for measuring the parameters encoded in a gravitational-wave signal

Once upon a time, predictions for the accuracy of inference on gravitational-wave signals relied on computationally inexpensive but often inaccurate techniques. Recently, the approach has shifted to actual inference on noisy signals with complex stochastic Bayesian methods, at the expense of significant computational cost. Here, we argue that it is often possible to have the best of both worlds: a Bayesian approach that incorporates prior information and correctly marginalizes over uninteresting parameters, providing accurate posterior probability distribution functions, but carried out on a simple grid at a low computational cost, comparable to the inexpensive predictive techniques.

preprint2015arXiv

Inference on gravitational waves from coalescences of stellar-mass compact objects and intermediate-mass black holes

Gravitational waves from coalescences of neutron stars or stellar-mass black holes into intermediate-mass black holes (IMBHs) of $\gtrsim 100$ solar masses represent one of the exciting possible sources for advanced gravitational-wave detectors. These sources can provide definitive evidence for the existence of IMBHs, probe globular-cluster dynamics, and potentially serve as tests of general relativity. We analyse the accuracy with which we can measure the masses and spins of the IMBH and its companion in intermediate-mass ratio coalescences. We find that we can identify an IMBH with a mass above $100 ~ M_\odot$ with $95\%$ confidence provided the massive body exceeds $130 ~ M_\odot$. For source masses above $\sim200 ~ M_\odot$, the best measured parameter is the frequency of the quasi-normal ringdown. Consequently, the total mass is measured better than the chirp mass for massive binaries, but the total mass is still partly degenerate with spin, which cannot be accurately measured. Low-frequency detector sensitivity is particularly important for massive sources, since sensitivity to the inspiral phase is critical for measuring the mass of the stellar-mass companion. We show that we can accurately infer source parameters for cosmologically redshifted signals by applying appropriate corrections. We investigate the impact of uncertainty in the model gravitational waveforms and conclude that our main results are likely robust to systematics.

preprint2015arXiv

Measuring intermediate mass black hole binaries with advanced gravitational wave detectors

We perform a systematic study to explore the accuracy with which the parameters of intermediate-mass black-hole binary systems can be measured from their gravitational wave (GW) signatures using second-generation GW detectors. We make use of the most recent reduced-order models containing inspiral, merger and ringdown signals of aligned-spin effective-one-body waveforms (SEOBNR) to significantly speed up the calculations. We explore the phenomenology of the measurement accuracies for binaries with total masses between 50 and 500 $M_\odot$ and mass ratios between 0.1 and 1. We find that (i) at total masses below ~200 $M_\odot$, where the signal-to-noise-ratio is dominated by the inspiral portion of the signal, the chirp mass parameter can be accurately measured; (ii) at higher masses, the information content is dominated by the ringdown, and total mass is measured more accurately; (iii) the mass of the lower-mass companion is poorly estimated, especially at high total mass and more extreme mass ratios; (iv) spin cannot be accurately measured for our injection set with non-spinning components. Most importantly, we find that for binaries with non-spinning components at all values of the mass ratio in the considered range and at network signal-to-noise ratio of 15, analyzed with spin-aligned templates, the presence of an intermediate-mass black hole with mass >100 $M_\odot$ can be confirmed with 95% confidence in any binary that includes a component with a mass of 130 $M_\odot$ or greater.

preprint2015arXiv

Parameter estimation for binary neutron-star coalescences with realistic noise during the Advanced LIGO era

Advanced ground-based gravitational-wave (GW) detectors begin operation imminently. Their intended goal is not only to make the first direct detection of GWs, but also to make inferences about the source systems. Binary neutron-star mergers are among the most promising sources. We investigate the performance of the parameter-estimation \edit{(PE)} pipeline that will be used during the first observing run of the Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) in 2015: we concentrate on the ability to reconstruct the source location on the sky, but also consider the ability to measure masses and the distance. Accurate, rapid sky-localization is necessary to alert electromagnetic (EM) observatories so that they can perform follow-up searches for counterpart transient events. We consider PE accuracy in the presence of \edit{non-stationary}, non-Gaussian noise. We find that the character of the noise makes negligible difference to the PE performance \edit{at a given signal-to-noise ratio}. The source luminosity distance can only be poorly constrained, the median $90\%$ ($50\%$) credible interval scaled with respect to the true distance is $0.85$ ($0.38$). However, the chirp mass is well measured. Our chirp-mass estimates are subject to systematic error because we used gravitational-waveform templates without component spin to carry out inference on signals with moderate spins, but the total error is typically less than $10^{-3} M_\odot$. The median $90\%$ ($50\%$) credible region for sky localization is $\sim600~\mathrm{deg^{2}}$ ($\sim150~\mathrm{deg^{2}}$), with $3\%$ ($30\%$) of detected events localized within $100~\mathrm{deg^{2}}$. Early aLIGO, with only two detectors, will have a sky-localization accuracy for binary neutron stars of hundreds of square degrees; this makes EM follow-up challenging, but not impossible.

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.

preprint2014arXiv

Counting And Confusion: Bayesian Rate Estimation With Multiple Populations

We show how to obtain a Bayesian estimate of the rates or numbers of signal and background events from a set of events when the shapes of the signal and background distributions are known, can be estimated, or approximated; our method works well even if the foreground and background event distributions overlap significantly and the nature of any individual event cannot be determined with any certainty. We give examples of determining the rates of gravitational-wave events in the presence of background triggers from a template bank when noise parameters are known and/or can be fit from the trigger data. We also give an example of determining globular-cluster shape, location, and density from an observation of a stellar field that contains a non-uniform background density of stars superimposed on the cluster stars.

preprint2014arXiv

Parameter estimation on compact binary coalescences with abruptly terminating gravitational waveforms

Gravitational-wave astronomy seeks to extract information about astrophysical systems from the gravitational-wave signals they emit. For coalescing compact-binary sources this requires accurate model templates for the inspiral and, potentially, the subsequent merger and ringdown. Models with frequency-domain waveforms that terminate abruptly in the sensitive band of the detector are often used for parameter-estimation studies. We show that the abrupt waveform termination contains significant information that affects parameter-estimation accuracy. If the sharp cutoff is not physically motivated, this extra information can lead to misleadingly good accuracy claims. We also show that using waveforms with a cutoff as templates to recover complete signals can lead to biases in parameter estimates. We evaluate when the information content in the cutoff is likely to be important in both cases. We also point out that the standard Fisher matrix formalism, frequently employed for approximately predicting parameter-estimation accuracy, cannot properly incorporate an abrupt cutoff that is present in both signals and templates; this observation explains some previously unexpected results found in the literature. These effects emphasize the importance of using complete waveforms with accurate merger and ringdown phases for parameter estimation.

preprint2014arXiv

Testing general relativity with compact coalescing binaries: comparing exact and predictive methods to compute the Bayes factor

The second generation of gravitational-wave detectors is scheduled to start operations in 2015. Gravitational-wave signatures of compact binary coalescences could be used to accurately test the strong-field dynamical predictions of general relativity. Computationally expensive data analysis pipelines, including TIGER, have been developed to carry out such tests. As a means to cheaply assess whether a particular deviation from general relativity can be detected, Cornish et al. and Vallisneri recently proposed an approximate scheme to compute the Bayes factor between a general-relativity gravitational-wave model and a model representing a class of alternative theories of gravity parametrised by one additional parameter. This approximate scheme is based on only two easy-to-compute quantities: the signal-to-noise ratio of the signal and the fitting factor between the signal and the manifold of possible waveforms within general relativity. In this work, we compare the prediction from the approximate formula against an exact numerical calculation of the Bayes factor using the lalinference library. We find that, using frequency-domain waveforms, the approximate scheme predicts exact results with good accuracy, providing the correct scaling with the signal-to-noise ratio at a fitting factor value of 0.992 and the correct scaling with the fitting factor at a signal-to-noise ratio of 20, down to a fitting factor of $\sim$ 0.9. We extend the framework for the approximate calculation of the Bayes factor which significantly increases its range of validity, at least to fitting factors of $\sim$ 0.7 or higher.

preprint2014arXiv

The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo

We anticipate the first direct detections of gravitational waves (GWs) with Advanced LIGO and Virgo later this decade. Though this groundbreaking technical achievement will be its own reward, a still greater prize could be observations of compact binary mergers in both gravitational and electromagnetic channels simultaneously. During Advanced LIGO and Virgo's first two years of operation, 2015 through 2016, we expect the global GW detector array to improve in sensitivity and livetime and expand from two to three detectors. We model the detection rate and the sky localization accuracy for binary neutron star (BNS) mergers across this transition. We have analyzed a large, astrophysically motivated source population using real-time detection and sky localization codes and higher-latency parameter estimation codes that have been expressly built for operation in the Advanced LIGO/Virgo era. We show that for most BNS events the rapid sky localization, available about a minute after a detection, is as accurate as the full parameter estimation. We demonstrate that Advanced Virgo will play an important role in sky localization, even though it is anticipated to come online with only one-third as much sensitivity as the Advanced LIGO detectors. We find that the median 90% confidence region shrinks from ~500 square degrees in 2015 to ~200 square degrees in 2016. A few distinct scenarios for the first LIGO/Virgo detections emerge from our simulations.

preprint2014arXiv

The Formation and Gravitational-Wave Detection of Massive Stellar Black-Hole Binaries

If binaries consisting of two 100 Msun black holes exist they would serve as extraordinarily powerful gravitational-wave sources, detectable to redshifts of z=2 with the advanced LIGO/Virgo ground-based detectors. Large uncertainties about the evolution of massive stars preclude definitive rate predictions for mergers of these massive black holes. We show that rates as high as hundreds of detections per year, or as low as no detections whatsoever, are both possible. It was thought that the only way to produce these massive binaries was via dynamical interactions in dense stellar systems. This view has been challenged by the recent discovery of several stars with mass above 150 Msun in the R136 region of the Large Magellanic Cloud. Current models predict that when stars of this mass leave the main sequence, their expansion is insufficient to allow common envelope evolution to efficiently reduce the orbital separation. The resulting black-hole--black-hole binary remains too wide to be able to coalesce within a Hubble time. If this assessment is correct, isolated very massive binaries do not evolve to be gravitational-wave sources. However, other formation channels exist. For example, the high multiplicity of massive stars, and their common formation in relatively dense stellar associations, opens up dynamical channels for massive black hole mergers (e.g., via Kozai cycles or repeated binary-single interactions). We identify key physical factors that shape the population of very massive black-hole--black-hole binaries. Advanced gravitational-wave detectors will provide important constraints on the formation and evolution of very massive stars.

preprint2014arXiv

The Occurrence of Earth-Like Planets Around Other Stars

The quantity $η_\oplus$, the number density of planets per star per logarithmic planetary radius per logarithmic orbital period at one Earth radius and one year period, describes the occurrence of Earth-like extrasolar planets. Here we present a measurement of $η_\oplus$ from a parameterised forward model of the (correlated) period-radius distribution and the observational selection function in the most recent (Q17) data release from the Kepler satellite. We find $η_\oplus = 3.9_{-1.6}^{+2.2}\%$ (90% CL). We conclude that each star hosts $3.83_{-0.62}^{+0.76}$ planets with $P \lesssim 3 \mathrm{yr}$ and $R \gtrsim 0.2 R_\oplus$. Our empirical model for false-positive contamination is consistent with the dominant source being background eclipsing binary stars. The distribution of planets we infer is consistent with a highly-stochastic planet formation process producing many correlated, fractional changes in planet sizes and orbits.

preprint2013arXiv

Characterising gravitational wave stochastic background anisotropy with Pulsar Timing Arrays

Detecting a stochastic gravitational wave background, particularly radiation from individually unresolvable super-massive black hole binary systems, is one of the primary targets for Pulsar Timing Arrays. Increasingly more stringent upper limits are being set on these signals under the assumption that the background radiation is isotropic. However, some level of anisotropy may be present and the characterisation of the power at different angular scales carries important information. We show that the standard analysis for isotropic backgrounds can be generalised in a conceptually straightforward way to the case of generic anisotropic background radiation by decomposing the angular distribution of the gravitational wave power on the sky into multipole moments. We introduce the concept of generalised overlap reduction functions which characterise the effect of the anisotropy multipoles on the correlation of the timing residuals from the pulsars timed by a Pulsar Timing Array. In a search for a signal characterised by a generic anisotropy, the generalised overlap reduction functions play the role of the so-called Hellings and Downs curve used for isotropic radiation. We compute the generalised overlap reduction functions for a generic level of anisotropy and Pulsar Timing Array configuration. We also provide an order of magnitude estimate of the level of anisotropy that can be expected in the background generated by super-massive black hole binary systems.

preprint2013arXiv

Constraining the luminosity function parameters and population size of radio pulsars in globular clusters

Studies of the Galactic population of radio pulsars have shown that their luminosity distribution appears to be log-normal in form. We investigate some of the consequences that occur when one applies this functional form to populations of pulsars in globular clusters. We use Bayesian methods to explore constraints on the mean and standard deviation of the luminosity function, as well as the total number of pulsars, given an observed sample of pulsars down to some limiting flux density, accounting for measurements of flux densities of individual pulsars as well as diffuse emission from the direction of the cluster. We apply our analysis to Terzan 5, 47 Tucanae and M 28, and demonstrate, under reasonable assumptions, that the number of potentially observable pulsars should be within 95% credible intervals of $147^{+112}_{-65}$, $83^{+54}_{-35}$ and $100^{+91}_{-52}$, respectively. Beaming considerations would increase the true population size by approximately a factor of two. Using non-informative priors, however, the constraints are not tight due to the paucity and quality of flux density measurements. Future cluster pulsar discoveries and improved flux density measurements would allow this method to be used to more accurately constrain the luminosity function, and to compare the luminosity function between different clusters.

preprint2013arXiv

Double Compact Objects II: Cosmological Merger Rates

The development of advanced gravitational wave (GW) observatories, such as Advanced LIGO and Advanced Virgo, provides impetus to refine theoretical predictions for what these instruments might detect. In particular, with the range increasing by an order of magnitude, the search for GW sources is extending beyond the "local" Universe and out to cosmological distances. Double compact objects (neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS) and black hole-black hole (BH-BH) systems) are considered to be the most promising gravitational wave sources. In addition, NS-NS and/or BH-NS systems are thought to be the progenitors of gamma ray bursts (GRBs), and may also be associated with kilonovae. In this paper we present the merger event rates of these objects as a function of cosmological redshift. We provide the results for four cases, each one investigating a different important evolution parameter of binary stars. Each case is also presented for two metallicity evolution scenarios. We find that (i) in most cases NS-NS systems dominate the merger rates in the local Universe, while BH-BH mergers dominate at high redshift; (ii) BH-NS mergers are less frequent than other sources per unit volume, for all time; and (iii) natal kicks may alter the observable properties of populations in a significant way, allowing the underlying models of binary evolution and compact object formation to be easily distinguished. This is the second paper in a series of three. The third paper will focus on calculating the detection rates of mergers by gravitational wave telescopes.

preprint2013arXiv

Electromagnetic transients as triggers in searches for gravitational waves from compact binary mergers

The detection of an electromagnetic transient which may originate from a binary neutron star merger can increase the probability that a given segment of data from the LIGO-Virgo ground-based gravitational-wave detector network contains a signal from a binary coalescence. Additional information contained in the electromagnetic signal, such as the sky location or distance to the source, can help rule out false alarms, and thus lower the necessary threshold for a detection. Here, we develop a framework for determining how much sensitivity is added to a gravitational-wave search by triggering on an electromagnetic transient. We apply this framework to a variety of relevant electromagnetic transients, from short GRBs to signatures of r-process heating to optical and radio orphan afterglows. We compute the expected rates of multi-messenger observations in the Advanced detector era, and find that searches triggered on short GRBs --- with current high-energy instruments, such as Fermi --- and nucleosynthetic `kilonovae' --- with future optical surveys, like LSST --- can boost the number of multi-messenger detections by 15% and 40%, respectively, for a binary neutron star progenitor model. Short GRB triggers offer precise merger timing, but suffer from detection rates decreased by beaming and the high a priori probability that the source is outside the LIGO-Virgo sensitive volume. Isotropic kilonovae, on the other hand, could be commonly observed within the LIGO-Virgo sensitive volume with an instrument roughly an order of magnitude more sensitive than current optical surveys. We propose that the most productive strategy for making multi-messenger gravitational-wave observations is using triggers from future deep, optical all-sky surveys, with characteristics comparable to LSST, which could make as many as ten such coincident observations a year.

preprint2013arXiv

Inadequacies of the Fisher Information Matrix in gravitational-wave parameter estimation

The Fisher Information Matrix (FIM) has been the standard approximation to the accuracy of parameter estimation on gravitational-wave signals from merging compact binaries due to its ease-of-use and rapid computation time. While the theoretical failings of this method, such as the signal-to-noise ratio (SNR) limit on the validity of the lowest-order expansion and the difficulty of using non-Gaussian priors, are well understood, the practical effectiveness compared to a real parameter estimation technique (e.g. Markov-chain Monte Carlo) remains an open question. We present a direct comparison between the FIM error estimates and the Bayesian probability density functions produced by the parameter estimation code lalinference_mcmc. In addition to the low-SNR issues usually considered, we find that the FIM can greatly overestimate the uncertainty in parameter estimation achievable by the MCMC. This was found to be a systematic effect for systems composed of binary black holes, with the disagreement increasing with total mass. In some cases, the MCMC search returned standard deviations on the marginalized posteriors that were smaller by several orders of magnitude than the FIM estimates. We conclude that the predictions of the FIM do not represent the capabilities of real gravitational-wave parameter estimation.

preprint2013arXiv

Rapidly evaluating the compact binary likelihood function via interpolation

Bayesian parameter estimation on gravitational waves from compact binary coalescences (CBCs) typically requires millions of template waveform computations at different values of the parameters describing the binary. Sampling techniques such as Markov chain Monte Carlo and nested sampling evaluate likelihoods, and hence compute template waveforms, serially; thus, the total computational time of the analysis scales linearly with that of template generation. Here we address the issue of rapidly computing the likelihood function of CBC sources with non-spinning components. We show how to efficiently compute the continuous likelihood function on the three-dimensional subspace of parameters on which it has a non-trivial dependence (the chirp mass, symmetric mass ratio and coalescence time) via interpolation. Subsequently, sampling this interpolated likelihood function is a significantly cheaper computational process than directly evaluating the likelihood; we report improvements in computational time of two to three orders of magnitude while keeping likelihoods accurate to < 0.025 %. Generating the interpolant of the likelihood function over a significant portion of the CBC mass space is computationally expensive but highly parallelizable, so the wall time can be very small relative to the time of a full parameter-estimation analysis.

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

The Transient Gravitational-Wave Sky

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

preprint2013arXiv

Utility of galaxy catalogs for following up gravitational waves from binary neutron star mergers with wide-field telescopes

The first detections of gravitational waves from binary neutron star mergers with advanced LIGO and Virgo observatories are anticipated in the next five years. These detections could pave the way for multi-messenger gravitational-wave (GW) and electromagnetic (EM) astronomy if GW triggers are successfully followed up with targeted EM observations. However, GW sky localization is relatively poor, with expected localization areas of ~10--100 square degrees; this presents a challenge for following up GW signals from compact binary mergers. Prioritizing wide-field follow-ups based on the probability of successful imaging is important when it is impossible to tile the entire gravitational-wave localization region. Galaxy catalogs were effective at narrowing down regions of the sky to search in initial attempts at joint GW/EM observations. The next generation of GW detectors will have a ten-fold increase in range thereby increasing the expected number of galaxies per unit solid angle a thousand-fold. As an additional complication, catalogs may be highly incomplete. We show how to quantify the advantages of using galaxy catalogs as a function of only two parameters: the three-dimensional volume within the field of view of a telescope after accounting for the GW distance measurement uncertainty, and the fraction of the GW sky localization uncertainty region that can be followed up. We find that the use of galaxy catalogs can improve success rates by ~10% to a factor of 4 relative to follow-up strategies that do not utilize such catalogs for the scenarios we considered. We determine that catalogs with a 75% completeness perform comparably to complete catalogs in most cases, while 33%-complete catalogs can lead to lower follow-up success rates than complete catalogs for small fields of view, though still providing an advantage over strategies that do not use a catalog at all.

preprint2012arXiv

Constraining the luminosity function parameters and population size of radio pulsars in globular clusters

The luminosity distribution of Galactic radio pulsars is believed to be log-normal in form. Applying this functional form to populations of pulsars in globular clusters, we employ Bayesian methods to explore constraints on the mean and standard deviation of the function, as well as the total number of pulsars in the cluster. Our analysis is based on an observed number of pulsars down to some limiting flux density, measurements of flux densities of individual pulsars, as well as diffuse emission from the direction of the cluster. We apply our analysis to Terzan 5 and demonstrate, under reasonable assumptions, that the number of potentially observable pulsars is in a 95.45% credible interval of 133$^{+101}_{-58}$. Beaming considerations would increase the true population size by approximately a factor of two.

preprint2012arXiv

Cyg X-3: a Galactic double black hole or black hole-neutron star progenitor

There are no known double black hole (BH-BH) or black hole-neutron star (BH-NS) systems. We argue that Cyg X-3 is a very likely BH-BH or BH-NS progenitor. This Galactic X-ray binary consists of a compact object, wind-fed by a Wolf-Rayet (WR) type companion. Based on a comprehensive analysis of observational data, it was recently argued that Cyg X-3 harbors a 2-4.5 Msun BH and a 7.5-14.2 Msun WR companion. We find that the fate of such a binary leads to the prompt (<1 Myr) formation of a close BH-BH system for the high end of the allowed WR mass (M_WR>13 Msun). For the low- to mid-mass range of the WR star (M_WR=7-10 Msun) Cyg X-3 is most likely (probability 70%) disrupted when WR ends up as a supernova. However, with smaller probability, it may form a wide (15%) or a close (15%) BH-NS system. The advanced LIGO/VIRGO detection rate for mergers of BH-BH systems from the Cyg X-3 formation channel is 10 per year, while it drops down to 0.1 per year for BH-NS systems. If Cyg X-3 in fact hosts a low mass BH and massive WR star, it lends additional support for the existence of BH-BH/BH-NS systems.

preprint2012arXiv

Double Compact Objects I: The Significance Of The Common Envelope On Merger Rates

The last decade of observational and theoretical developments in stellar and binary evolution provides an opportunity to incorporate major improvements to the predictions from populations synthesis models. We compute the Galactic merger rates for NS-NS, BH-NS, and BH-BH mergers with the StarTrack code. The most important revisions include: updated wind mass loss rates (allowing for stellar mass black holes up to $80 \msun$), a realistic treatment of the common envelope phase (a process that can affect merger rates by 2--3 orders of magnitude), and a qualitatively new neutron star/black hole mass distribution (consistent with the observed "mass gap"). Our findings include: (i) The binding energy of the envelope plays a pivotal role in determining whether a binary merges within a Hubble time. (ii) Our description of natal kicks from supernovae plays an important role, especially for the formation of BH-BH systems. (iii) The masses of BH-BH systems can be substantially increased in the case of low metallicities or weak winds. (iv) Certain combinations of parameters underpredict the Galactic NS-NS merger rate, and can be ruled out. {\em (v)} Models incorporating delayed supernovae do not agree with the observed NS/BH "mass gap", in accordance with our previous work. This is the first in a series of three papers. The second paper will study the merger rates of double compact objects as a function of redshift, star formation rate, and metallicity. In the third paper we will present the detection rates for gravitational wave observatories, using up-to-date signal waveforms and sensitivity curves.

preprint2012arXiv

Effect of calibration errors on Bayesian parameter estimation for gravitational wave signals from inspiral binary systems in the advanced detectors era

By 2015 the advanced versions of the gravitational-wave detectors Virgo and LIGO will be online. They will collect data in coincidence with enough sensitivity to potentially deliver multiple detections of gravitation waves from inspirals of compact-object binaries. This work is focused on understanding the effects introduced by uncertainties in the calibration of the interferometers. We consider plausible calibration errors based on estimates obtained during LIGO's fifth and Virgo's third science runs, which include frequency-dependent amplitude errors of $\sim 10%$ and frequency-dependent phase errors of $\sim 3$ degrees in each instrument. We quantify the consequences of such errors estimating the parameters of inspiraling binaries. We find that the systematics introduced by calibration errors on the inferred values of the chirp mass and mass ratio are smaller than 20% of the statistical measurement uncertainties in parameter estimation for 90% of signals in our mock catalog. Meanwhile, the calibration-induced systematics in the inferred sky location of the signal are smaller than $\sim 50%$ of the statistical uncertainty. We thus conclude that calibration-induced errors at this level are not a significant detriment to accurate parameter estimation.

preprint2012arXiv

Hubble without the Hubble: cosmology using advanced gravitational-wave detectors alone

We investigate a novel approach to measuring the Hubble constant using gravitational-wave (GW) signals from compact binaries by exploiting the narrowness of the distribution of masses of the underlying neutron-star population. Gravitational-wave observations with a network of detectors will permit a direct, independent measurement of the distance to the source systems. If the redshift of the source is known, these inspiraling double-neutron-star binary systems can be used as standard sirens to extract cosmological information. Unfortunately, the redshift and the system chirp mass are degenerate in GW observations. Thus, most previous work has assumed that the source redshift is obtained from electromagnetic counterparts. In this paper, we explore what we can learn about the background cosmology and the mass distribution of neutron stars from the set of neutron-star (NS) mergers detected by such a network. We use a Bayesian formalism to analyze catalogs of NS-NS inspiral detections. We find that it is possible to constrain the Hubble constant, H_0, and the parameters of the NS mass function using gravitational-wave data alone, without relying on electromagnetic counterparts. Under reasonable assumptions, we will be able to determine H_0 to +/- 10% using ~100 observations, provided the Gaussian half-width of the underlying double NS mass distribution is less than 0.04 solar masses. The expected precision depends linearly on the intrinsic width of the NS mass function, but has only a weak dependence on H_0 near the default parameter values. Finally, we consider what happens if, for some fraction of our data catalog, we have an electromagnetically measured redshift. The detection, and cataloging, of these compact-object mergers will allow precision astronomy, and provide a determination of H_0 which is independent of the local distance scale.

preprint2012arXiv

Verifying the no-hair property of massive compact objects with intermediate-mass-ratio inspirals in advanced gravitational-wave detectors

The detection of gravitational waves from the inspiral of a neutron star or stellar-mass black hole into an intermediate-mass black hole (IMBH) promises an entirely new look at strong-field gravitational physics. Gravitational waves from these intermediate-mass-ratio inspirals (IMRIs), systems with mass ratios from ~10:1 to ~100:1, may be detectable at rates of up to a few tens per year by Advanced LIGO/Virgo and will encode a signature of the central body's spacetime. Direct observation of the spacetime will allow us to use the "no-hair" theorem of general relativity to determine if the IMBH is a Kerr black hole (or some more exotic object, e.g. a boson star). Using modified post-Newtonian (pN) waveforms, we explore the prospects for constraining the central body's mass-quadrupole moment in the advanced-detector era. We use the Fisher information matrix to estimate the accuracy with which the parameters of the central body can be measured. We find that for favorable mass and spin combinations, the quadrupole moment of a non-Kerr central body can be measured to within a ~15% fractional error or better using 3.5 pN order waveforms; on the other hand, we find the accuracy decreases to ~100% fractional error using 2 pN waveforms, except for a narrow band of values of the best-fit non-Kerr quadrupole moment.

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

The Mass Distribution of Stellar-Mass Black Holes

We perform a Bayesian analysis of the mass distribution of stellar-mass black holes using the observed masses of 15 low-mass X-ray binary systems undergoing Roche lobe overflow and five high-mass, wind-fed X-ray binary systems. Using Markov Chain Monte Carlo calculations, we model the mass distribution both parametrically---as a power law, exponential, gaussian, combination of two gaussians, or log-normal distribution---and non-parametrically---as histograms with varying numbers of bins. We provide confidence bounds on the shape of the mass distribution in the context of each model and compare the models with each other by calculating their relative Bayesian evidence as supported by the measurements, taking into account the number of degrees of freedom of each model. The mass distribution of the low-mass systems is best fit by a power-law, while the distribution of the combined sample is best fit by the exponential model. We examine the existence of a "gap" between the most massive neutron stars and the least massive black holes by considering the value, M_1%, of the 1% quantile from each black hole mass distribution as the lower bound of black hole masses. The best model (the power law) fitted to the low-mass systems has a distribution of lower-bounds with M_1% > 4.3 Msun with 90% confidence, while the best model (the exponential) fitted to all 20 systems has M_1% > 4.5 Msun with 90% confidence. We conclude that our sample of black hole masses provides strong evidence of a gap between the maximum neutron star mass and the lower bound on black hole masses. Our results on the low-mass sample are in qualitative agreement with those of Ozel, et al (2010).

preprint2011arXiv

Towards improving the prospects for coordinated gravitational-wave and electromagnetic observations

We discuss two approaches to searches for gravitational-wave (GW) and electromagnetic (EM) counterparts of binary neutron star mergers. The first approach relies on triggering archival searches of GW detector data based on detections of EM transients. We introduce a quantitative approach to evaluate the improvement to GW detector reach due to the extra information gained from the EM transient and the increased confidence in the presence of a signal from a binary merger. We also advocate utilizing other transients in addition to short gamma ray bursts. The second approach involves following up GW candidates with targeted EM observations. We argue for the use of slower but optimal parameter-estimation techniques to localize the source on the sky, and for a more sophisticated use of astrophysical prior information, including galaxy catalogs, to find preferred followup locations.

preprint2010arXiv

Compact Binary Coalescences in the Band of Ground-based Gravitational-Wave Detectors

As the ground-based gravitational-wave telescopes LIGO, Virgo, and GEO 600 approach the era of first detections, we review the current knowledge of the coalescence rates and the mass and spin distributions of merging neutron-star and black-hole binaries. We emphasize the bi-directional connection between gravitational-wave astronomy and conventional astrophysics. Astrophysical input will make possible informed decisions about optimal detector configurations and search techniques. Meanwhile, rate upper limits, detected merger rates, and the distribution of masses and spins measured by gravitational-wave searches will constrain astrophysical parameters through comparisons with astrophysical models. Future developments necessary to the success of gravitational-wave astronomy are discussed.

preprint2010arXiv

Detecting coalescences of intermediate-mass black holes in globular clusters with the Einstein Telescope

We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope to detect mergers of intermediate-mass black holes that may have formed through runaway stellar collisions in globular clusters. We find that detection rates of 500 events per year are plausible.

preprint2010arXiv

Exploring intermediate and massive black-hole binaries with the Einstein Telescope

We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope (ET) to enhance our astrophysical knowledge through detections of gravitational waves emitted by binaries including intermediate-mass and massive black holes. The design target for such instruments calls for improved sensitivity at low frequencies, specifically in the ~ 1-10 Hz range. This will allow the detection of gravitational waves generated in binary systems containing black holes of intermediate mass, ~ 100-1000 solar masses. We primarily discuss two different source types -- mergers between two intermediate mass black holes (IMBHs) of comparable mass, and intermediate-mass-ratio inspirals (IMRIs) of smaller compact objects with mass ~ 1-10 solar masses into IMBHs. IMBHs may form via two channels: (i) in dark matter halos at high redshift through direct collapse or the collapse of very massive metal-poor Population III stars, or (ii) via runaway stellar collisions in globular clusters. In this paper, we will discuss both formation channels, and both classes of merger in each case. We review existing rate estimates where these exist in the literature, and provide some new calculations for the approximate numbers of events that will be seen by a detector like the Einstein Telescope. These results indicate that the ET may see a few to a few thousand comparable-mass IMBH mergers and as many as several hundred IMRI events per year. These observations will significantly enhance our understanding of galactic black-hole growth, of the existence and properties of IMBHs and of the astrophysics of globular clusters. We finish our review with a discussion of some more speculative sources of gravitational waves for the ET, including hypermassive white dwarfs and eccentric stellar-mass compact-object binaries.

preprint2010arXiv

The Distribution of Coalescing Compact Binaries in the Local Universe: Prospects for Gravitational-Wave Observations

Merging compact binaries are the most viable and best studied candidates for gravitational wave (GW) detection by the fully operational network of ground-based observatories. In anticipation of the first detections, the expected distribution of GW sources in the local universe is of considerable interest. Here we investigate the full phase space distribution of coalescing compact binaries at $z = 0$ using dark matter simulations of structure formation. The fact that these binary systems acquire large barycentric velocities at birth ("kicks") results in merger site distributions that are more diffusely distributed with respect to their putative hosts, with mergers occurring out to distances of a few Mpc from the host halo. Redshift estimates based solely on the nearest galaxy in projection can, as a result, be inaccurate. On the other hand, large offsets from the host galaxy could aid the detection of faint optical counterparts and should be considered when designing strategies for follow-up observations. The degree of isotropy in the projected sky distributions of GW sources is found to be augmented with increasing kick velocity and to be severely enhanced if progenitor systems possess large kicks as inferred from the known population of pulsars and double compact binaries. Even in the absence of observed electromagnetic counterparts, the differences in sky distributions of binaries produced by disparate kick-velocity models could be discerned by GW observatories, within the expected accuracies and detection rates of advanced LIGO--in particular with the addition of more interferometers.

preprint2010arXiv

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise.

preprint2010arXiv

Unraveling Binary Evolution from Gravitational-Wave Signals and Source Statistics

The next generation of ground-based gravitational-wave detectors are likely to observe gravitational waves from the coalescences of compact-objects binaries. We describe the state of the art for predictions of the rate of compact-binary coalescences and report on initial efforts to develop a framework for converting gravitational-wave observations into improved constraints on astrophysical parameters.

preprint2009arXiv

Massive Black Hole Binary Inspirals: Results from the LISA Parameter Estimation Taskforce

The LISA Parameter Estimation (LISAPE) Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models, and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large, and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

preprint2009arXiv

Observing gravitational waves from the first generation of black holes

The properties of the first generation of black-hole seeds trace and distinguish different models of formation of cosmic structure in the high-redshift universe. The observational challenge lies in identifying black holes in the mass range ~100-1000 solar masses at redshift z~10. The typical frequencies of gravitational waves produced by the coalescence of the first generation of light seed black-hole binaries fall in the gap between the spectral ranges of low-frequency space-borne detectors (e.g., LISA) and high-frequency ground-based detectors (e.g., LIGO, Virgo and GEO 600). As such, these sources are targets for proposed third-generation ground-based instruments, such as the Einstein Telescope which is currently in design study. Using galaxy merger trees and four different models of black hole accretion - which are meant to illustrate the potential of this new type of source rather than to yield precise event-rate predictions - we find that such detectors could observe a few to a few tens of seed black-hole merger events in three years and provide, possibly unique, information on the evolution of structure in the corresponding era. We show further that a network of detectors may be able to measure the luminosity distance to sources to a precision of ~40%, allowing us to be confident of the high-redshift nature of the sources.

preprint2009arXiv

Parameter estimation on gravitational waves from multiple coalescing binaries

Future ground-based and space-borne interferometric gravitational-wave detectors may capture between tens and thousands of binary coalescence events per year. There is a significant and growing body of work on the estimation of astrophysically relevant parameters, such as masses and spins, from the gravitational-wave signature of a single event. This paper introduces a robust Bayesian framework for combining the parameter estimates for multiple events into a parameter distribution of the underlying event population. The framework can be readily deployed as a rapid post-processing tool.

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.

preprint2008arXiv

Rates and Characteristics of Intermediate Mass Ratio Inspirals Detectable by Advanced LIGO

Gravitational waves (GWs) from the inspiral of a neutron star (NS) or stellar-mass black hole (BH) into an intermediate-mass black hole (IMBH) with mass between ~50 and ~350 solar masses may be detectable by the planned advanced generation of ground-based GW interferometers. Such intermediate mass ratio inspirals (IMRIs) are most likely to be found in globular clusters. We analyze four possible IMRI formation mechanisms: (1) hardening of an NS-IMBH or BH-IMBH binary via three-body interactions, (2) hardening via Kozai resonance in a hierarchical triple system, (3) direct capture, and (4) inspiral of a compact object from a tidally captured main-sequence star; we also discuss tidal effects when the inspiraling object is an NS. For each mechanism we predict the typical eccentricities of the resulting IMRIs. We find that IMRIs will have largely circularized by the time they enter the sensitivity band of ground-based detectors. Hardening of a binary via three-body interactions, which is likely to be the dominant mechanism for IMRI formation, yields eccentricities under 10^-4 when the GW frequency reaches 10 Hz. Even among IMRIs formed via direct captures, which can have the highest eccentricities, around 90% will circularize to eccentricities under 0.1 before the GW frequency reaches 10 Hz. We estimate the rate of IMRI coalescences in globular clusters and the sensitivity of a network of three Advanced LIGO detectors to the resulting GWs. We show that this detector network may see up to tens of IMRIs per year, although rates of one to a few per year may be more plausible. We also estimate the loss in signal-to-noise ratio that will result from using circular IMRI templates for data analysis and find that, for the eccentricities we expect, this loss is negligible.

Ilya Mandel

What is connected

Connect this record

See the researcher in context

Building this map preview

77 published item(s)

Converged simulations of the nozzle shock in tidal disruption events

Building better spin models for merging binary black holes: Evidence for non-spinning and rapidly spinning nearly aligned sub-populations

Common envelopes in massive stars II: The distinct roles of hydrogen and helium recombination

Common envelopes in massive stars: Towards the role of radiation pressure and recombination energy in ejecting red supergiant envelopes

Extreme eccentricities of triple systems: Analytic results

Merging stellar-mass binary black holes

Rates of Compact Object Coalescences

Stellar response after stripping as a model for common-envelope outcomes

Cygnus X-1 contains a 21-solar mass black hole -- implications for massive star winds

Impact of Massive Binary Star and Cosmic Evolution on Gravitational Wave Observations II: Double Compact Object Rates and Properties

Rapid stellar and binary population synthesis with COMPAS

Wind mass-loss rates of stripped stars inferred from Cygnus X-1

Age dissection of the Milky Way discs: red giants in the Kepler field

An alternative interpretation of GW190412 as a binary black hole merger with a rapidly spinning secondary

Be X-ray binaries in the SMC as indicators of mass transfer efficiency

Constraining the masses of microlensing black holes and the mass gap with Gaia DR2

Detecting Double Neutron Stars with LISA

Luminous Red Novae: population models and future prospects

Merger rates in primordial black hole clusters without initial binaries

The origin of spin in binary black holes: Predicting the distributions of the main observables of Advanced LIGO

The spectral evolution of AT 2018dyb and the presence of metal lines in tidal disruption events

Astrophysical science metrics for next-generation gravitational-wave detectors

Common Envelope Wind Tunnel: The Effects of Binary Mass Ratio and Implications for the Accretion-Driven Growth of LIGO Binary Black Holes

Eclipses of continuous gravitational waves as a probe of stellar structure

Explosions Driven by the Coalescence of a Compact Object with the Core of a Massive-Star Companion Inside a Common Envelope: Circumstellar Properties, Light Curves, and Population Statistics

Comment on "An excess of massive stars in the local 30 Doradus starburst"

$N-$body dynamics of Intermediate mass-ratio inspirals in globular clusters

Detectability of Gravitational Waves from High-Redshift Binaries

Dynamic temperature selection for parallel-tempering in Markov chain Monte Carlo simulations

Kernel regression estimates of time delays between gravitationally lensed fluxes

Limits of Astrophysics with Gravitational-Wave Backgrounds

Merging binary black holes formed through chemically homogeneous evolution in short-period stellar binaries

Parameter estimation on gravitational waves from neutron-star binaries with spinning components

An Efficient Interpolation Technique for Jump Proposals in Reversible-Jump Markov Chain Monte Carlo Calculations

Distinguishing types of compact-object binaries using the gravitational-wave signatures of their mergers

Double tidal disruptions in galactic nuclei

Efficient method for measuring the parameters encoded in a gravitational-wave signal

Inference on gravitational waves from coalescences of stellar-mass compact objects and intermediate-mass black holes

Measuring intermediate mass black hole binaries with advanced gravitational wave detectors

Parameter estimation for binary neutron-star coalescences with realistic noise during the Advanced LIGO era

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

Counting And Confusion: Bayesian Rate Estimation With Multiple Populations

Parameter estimation on compact binary coalescences with abruptly terminating gravitational waveforms

Testing general relativity with compact coalescing binaries: comparing exact and predictive methods to compute the Bayes factor

The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo

The Formation and Gravitational-Wave Detection of Massive Stellar Black-Hole Binaries

The Occurrence of Earth-Like Planets Around Other Stars

Characterising gravitational wave stochastic background anisotropy with Pulsar Timing Arrays

Constraining the luminosity function parameters and population size of radio pulsars in globular clusters

Double Compact Objects II: Cosmological Merger Rates

Electromagnetic transients as triggers in searches for gravitational waves from compact binary mergers

Inadequacies of the Fisher Information Matrix in gravitational-wave parameter estimation

Rapidly evaluating the compact binary likelihood function via interpolation

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

The Transient Gravitational-Wave Sky

Utility of galaxy catalogs for following up gravitational waves from binary neutron star mergers with wide-field telescopes

Constraining the luminosity function parameters and population size of radio pulsars in globular clusters

Cyg X-3: a Galactic double black hole or black hole-neutron star progenitor

Double Compact Objects I: The Significance Of The Common Envelope On Merger Rates

Effect of calibration errors on Bayesian parameter estimation for gravitational wave signals from inspiral binary systems in the advanced detectors era

Hubble without the Hubble: cosmology using advanced gravitational-wave detectors alone

Verifying the no-hair property of massive compact objects with intermediate-mass-ratio inspirals in advanced gravitational-wave detectors

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

The Mass Distribution of Stellar-Mass Black Holes

Towards improving the prospects for coordinated gravitational-wave and electromagnetic observations

Compact Binary Coalescences in the Band of Ground-based Gravitational-Wave Detectors

Detecting coalescences of intermediate-mass black holes in globular clusters with the Einstein Telescope

Exploring intermediate and massive black-hole binaries with the Einstein Telescope

The Distribution of Coalescing Compact Binaries in the Local Universe: Prospects for Gravitational-Wave Observations

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

Unraveling Binary Evolution from Gravitational-Wave Signals and Source Statistics

Massive Black Hole Binary Inspirals: Results from the LISA Parameter Estimation Taskforce

Observing gravitational waves from the first generation of black holes

Parameter estimation on gravitational waves from multiple coalescing binaries