Science summaries

Below are some summaries of selected recent publications, written for non-experts. These highlight some of the range of research done by the group.

The LIGO Scientific Collaboration also puts together science summaries, written for the general public, for all their publications.

Model-independent inference on compact-binary observations

Mean number density for different bins in mass space

The increasing number of gravitational-wave observations promises to lead to insights into stellar binary evolution. One approach is to compare specific evolution models against observations; however, given the uncertainty in the models, a model-independent approach to astrophysical inference is desirable. In an earlier paper, we advocated clustering in a multi-dimensional parameter space. Here, we propose a practical procedure to clustering observations which are subject to significant measurement errors. We apply our method to a mock data set of population-synthesis predictions for compact binaries incorporating realistic measurement uncertainties. We demonstrate that a few tens of observations are sufficient to accurately distinguish subpopulations of sources! The example to the left shows what we could do with 160 observations, you can make out clusters for binary neutron stars (in the bottom left corner), neutron star–black hole binaries (along the edges) and binary black holes (in the centre).

Authors: Mandel, I., Farr, W. M., Colonna, A., Stevenson, S., Tiňo, P. & Veitch, J.
Journal: MNRAS
arXiv: 1608.08223 [astro-ph.HE]

Importance of transient resonances in extreme-mass-ratio inspirals

Jump in orbital energy across transient resonance

Extreme-mass-ratio inspirals occur when a compact object, like a solar-mass black hole, orbits around a much larger one, such as the supermassive black holes found in the centres of galaxies. The gravitational waves emitted by these inspirals are a prime target for space-bourne detectors like LISA; they can complete tens of thousands of cycles, allowing us to make exquisite measurements of the system's properties. To detect and analyse these signals, we need waveform templates that are accurate over these long signals.

Waveforms have been previously calculated using the adiabatic approximation, which averages the rate of inspiral over many orbits. This approximation breaks down if the radial and polar orbital frequencies are multiples of each other: the inspiral will pass through many such transient resonances as it evolves.

We studied the effects of transient resonances on extreme-mass-ratio inspirals. The plot to the left shows the kick received (compared to the adiabatic approximation) when passing through the resonance. We found that for most inspirals the effect of resonance is small enough not to impact detectability, we'd still be able to find about 96% of signals with our adiabatic models! Resonances could still have a big impact on parameter estimation (possibly even improving measurements if we can model them properly), but this is something still to be investigated.

Authors: Berry, C.P.L, Cole, R.H., Cañizares, P. & Gair, J.R.
Journal: PRD
arXiv: 1608.08951 [gr-qc]

N-body dynamics of intermediate mass-ratio inspirals in globular clusters

Cluster interactions

Stellar-mass black holes of a few to a few tens of solar masses are known to exist, as are million to billion solar-mass behemoths. What about intermediate-mass black holes weighing in at a few hundred solar masses? It is difficult to convincingly observe these electromagnetically, but it may be possible to detect them using gravitational waves if they merge with other intermediate-mass black holes or stellar-mass black holes.

In this study, we numerically investigate the dynamics of putative globular clusters that harbour intermediate-mass black holes. We include relativistic dynamics (for the first time incorporating terms up to 2.5 post-Newtonian order) in the vicinity of an intermediate-mass black hole. Our simulations show that the intermediate-mass black hole readily forms a binary with a stellar-mass black hole companion. This binary is gradually hardened (becoming more tightly bound) by transient 3-body or 4-body encounters. These also cause frequent substitutions of the black hole companion (as shown in the figure on the left). We find that one of the dozen clusters we simulated will have an intermediate-mass-ratio merger between an intermediate-mass and a stellar-mass black hole within 100 million years.

Authors: Haster, C.-J., Antonini, F., Kalogera, V. & Mandel, I.
Journal: ApJ
arXiv: 1606.07097 [astro-ph.HE]

On tests of general relativity with binary radio pulsars

Probabilities for different constraintson gravity

The timing of radio pulsars in binary systems provides one of the best testing grounds of the predictions of general relativity. Traditional methods are based on least-squares fits that do not account for the actual structure of the likelihood function, and this can potentially bias the inference of the parameters of the binary system. In this paper, we propose a Bayesian approach to carry out these tests of general relativity. Using simulated data sets that are representative of current observations, we demonstrate the power of this framework by (i) using the ratio of the evidences of two models to compare different theories of gravity; (ii) combining constraints from multiple systems to improve our constraints on any theory of gravity. The plot above shows probability distributions of a set of orbital parameters assuming different theories of gravity for a simulated system similar to the double pulsar. Applying our approach to real data will allow us to extract as much information as possible about the nature of gravity.

Authors: Del Pozzo, W., & Vecchio, A.
Journal: MNRASL
arXiv: 1606.02852 [gr-qc]

Testing general relativity using golden black-hole binaries

Consistency of final mass and spin from inspiral and merger

A golden black-hole binary is one where we can measure both the inspiral and post-inspiral (merger and ringdown) portions of the gravitational-wave signal. In this paper, we suggest using golden binaries as tests of general relativity, since they provide a glimpse of the strongest gravitational fields. We measure properties of the binary (we use the mass and spin of the final black hole), from the inspiral and post-inspiral parts of the waveform, and check for consistency. If our general relativity waveform models are correct, we should get compatible results from both, but a difference might hint at new physics. The test could be thought of as making measurements with two different instruments (the low frequency inspiral detector, and the high frequency post-inspiral detector), and seeing that they agree. The plots on the left show an application of the test, the left panels show a general relativity signal with compatible estimates, and the the right shows a modified theory where the estimates are inconsistent. We find that the test is sensitive to some deviations that could otherwise be missed, and also has the advantage that we can combine results from many binaries to get a tougher test of Einstein's theory!

This was one of the tests used for GW150914, and no inconsistencies were found.

Authors: Ghosh, A., Ghosh, A., Johnson-McDaniel, N. K., Kant Mishra, C., Ajith, P., Del Pozzo, W., Nichols, D. A., Chen, Y., Nielsen, A.B., Berry, C. P. L. & London, L.
Journal: PRD
arXiv: 1602.02453 [gr-qc]

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

Uncertainty is mass measurements with different spin cut-offs

Binary neutron-star coalescences are one of the most exciting sources for Advanced LIGO. It is important that we understand how to analyse these signals and measure their parameters. In this study, we analysed a set of fake binary neutron-star signals and investigated how including the neutron stars' spins effected things. We expect that neutrons stars will rotate, but that their spins should be small. Including spins complicates the analysis and makes it much more computationally expensive (taking around 20 times as long), so it is often ignored. We find that allowing for spins increases uncertainty in measurements of the masses, as shown in the plot. Typical uncertainty of the total binary mass is around 6%, and on the chirp mass it is 0.06%, so even with the extra uncertainty, we can still make precise measurements! The sky localization is unaffected by the spin, so using the quicker results without spin is absolutely fine.

Authors: Farr, B., Berry, C. P. L., Farr, W. M., Haster, C.-J., Middleton, H., Cannon, K., Graff, P. B., Hanna, C., Mandel, I., Pankow, C., Price, L. R., Sidery, T., Singer, L. P., Urban, A. L., Vecchio, A., Veitch, J. & Vitale, S.
Journal: ApJ
arXiv: 1508.05336 [astro-ph.HE]
Authorea: Parameter estimation on gravitational waves from neutron-star binaries with spinning components

Prospects for multiband gravitational-wave astronomy after GW150914

Multi-band gravitational-wave astronomy and the evolution of binary black holes

The first observation of gravitational waves from a black hole binary, GW150914, has important consequences for the future of gravitational-wave astronomy beyond ground-based detectors. Because of its high mass (~30+35 solar masses), GW150914 would have been observable by a space-based gravitational-wave detector such as the planned eLISA. With eLISA, the system would have been detectable from several year up to few weeks before coalescence, as shown by the black line in the plot to the left. Merger rates inferred from the first aLIGO detections imply that up to few hundreds similar events would be observable both with eLISA and aLIGO in the future, opening the way to multi-band gravitational-wave astronomy. This has a number of outstanding consequences for astrophysics, tests of gravity and multimessenger astronomy. In particular, eLISA observations can inform aLIGO and electromagnetic telescopes about when and where in the sky black hole merger events are going to occur with several weeks' notice, allowing us to realise the full potential of coincident gravitational-wave and electromagnetic observations.

Author: Sesana, A.
Journal: PRL
arXiv: 1602.06951 [gr-qc]

Dynamical formation of the GW150914 binary black hole

How to form GW150914's source

Binary black holes, like the source of GW150914, are believed to be able to form in two different ways: isolated binary evolution and dynamically interacting binaries. The latter is expected to dominate in dense stellar environments such as globular clusters, and we show in this paper that GW150914 is consistent with coming from this formation channel.

By simulating the gravitational interaction between hundreds of thousands of stars and black holes in a typical globular cluster, while also keeping track of each star's evolution, we can find how many binary black holes of a certain mass can form per cluster. One example of such a formation chain is shown in the figure on the left. In this paper, we used 48 such cluster models, and after comparing our simulation with observations of globular clusters across the Universe we can calculate how many binary black hole mergers are expected to be produced in globular clusters. Combining this result with how well Advanced LIGO is able to detect gravitational waves from these mergers tells us how often we should be able to observe a binary black hole merger from a globular cluster.

We show that the detection of GW150914 is consistent with the masses and merger rates expected for binary black hole from globular clusters. Dynamical processes provide a significant (and well understood) pathway for forming binary black hole mergers in the local Universe. Understanding how binary black hole systems form is one of the key questions for gravitational-wave astronomy.

Authors: Rodriguez, C. L., Haster, C.-J., Chatterjee, S., Kalogera, V. & Rasio, F. A.
Journal: ApJL
arXiv: 1604.04254 [astro-ph.HE]

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

Chemically homogenous evolution of a binary star

Merging stellar-mass binary black holes, such as those recently observed by LIGO, could arise from several possible formation channels. These can be broadly split into two categories: isolated binaries and dynamically interacting binaries (in dense stellar environments such as globular clusters). The isolated binary channel was previously expected to require several episodes of mass transfer, where material is transfered from a star to its partner as it puffs up towards the end of its life, and a common envelope phase where so much material is shared by the stars that their cores orbit inside a combined outer layer.

Stars normally expand as they evolve, but this is not so if the stars are rapidly rotating and subject to strong tidal forces. Massive, low-metallicity stars in compact binaries will be efficiently mixed, leading them to evolve chemically homogeneously. As a consequence, they shrink as they age and avoid mass transfer altogether, forming black holes in situ.

In this paper and its follow-up, we show that chemically homogeneous evolution represents a viable formation channel for heavy black hole binaries, with plausible local merger rates of order 10 per Gpc3 per year. The resulting black holes are nearly equal mass, each about 25 to 55 solar masses, with time delays of 4 to 10 Gyr between formation and merger. We show that the advanced gravitational-wave detectors at design sensitivity could observe hundreds of events per year from this channel, and that GW150914's source is consistent with having been formed this way.

Authors: Mandel, I. & de Mink, S. E.
Journal: MNRAS
arXiv: 1601.00007 [astro-ph.HE]
Paper II: 1603.02291 [astro-ph.HE]
Press release: Merging binary black holes formed through chemically homogeneous evolution

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

Probability for the existance of an intermediate-mass black hole

Intermediate-mass black holes are an elusive class of black holes with masses between stellar-mass black holes (which have masses up to a few tens of solar masses) and supermassive black holes (ten of thousands of solar masses). There is no reason why there should not be any black holes in this intermediate-mass regime; however, we are yet to find any conclusive proof for their existence.

If an intermediate-mass black hole and a stellar-mass black hole could form a binary, which is expected to commonly occur in dense stellar environments (such as globular clusters), and this binary merged from emission of gravitational waves then it could be detectable by LIGO. This type of signal would be expected to be very short and only exist at the low-frequency end of LIGO's sensitivity. Therefore, it might be expected that we could not measure much about the source binary; in this paper we show that this is not the case!

We find that we can confirm the existence of an intermediate-mass black hole (with mass greater than 100 solar masses) with 95% confidence provided the massive body exceeds 130 solar masses, as shown in the figure on the left. Like for intermediate-mass black hole binaries, we find that systems with total mass above ~200 solar masses, the best measured parameter is the frequency of the ringdown of the black hole formed after the merger, as it's settling down into its final state. Consequently, these binaries have a well constrained total mass, whereas lower mass binaries are better described by their chirp mass. The mass of the stellar-mass companion is only measured well during the inspiral phase, which further highlights the importance on reliable low-frequency detector sensitivity.

Authors: Haster, C.-J., Wang, Z., Berry, C. P. L., Stevenson, S., Veitch, J. & Mandel, I.
Journal: MNRAS
arXiv: 1511.01431 [astro-ph.HE]

Improving gravitational-wave parameter estimation using Gaussian process regression

Likelihoods using different waveform models

When we analyse a gravitational-wave signal with our parameter-estimation codes, we need accurate templates for the waveforms. If the template is imperfect (because it does not include all the relevant physics, it relies upon an approximation or it is from a simulation with finite numerical accuracy) it could produce errors in our inferred parameters. This theoretical uncertainty is the dominant source of error for loud signals.

In this paper we introduce a means to include waveform uncertainty in parameter estimation. This could be useful, for example, if you have a few highly accurate waveforms, but they are too expensive to calculate for all the parameter values needed in the analysis. Our technique is based on Gaussian process regression (GPR), which is a flexible way of interpolating functions—in our case waveforms. We give a detailed introduction to GPR. The plot above shows how using GPR can account for the effects of waveform uncertainty: the black curves show results using a perfect waveform, the orange curves show results with an inaccurate waveform and the blue curve shows the results with the inaccurate waveform and GPR. GPR can remove bias and incorporate theoretical uncertainty (broadening distributions).

Authors: Moore, C. J., Berry, C. P. L., Chua, A. J. K. & Gair, J. R.
Journal: PRD
arXiv: 1509.04066 [gr-qc]

Fast simulation of Gaussian-mode scattering for precision interferometry

Time taken to calculate model of LIGO interferometer. Laser beams are used in some of the most precise displacement measurements possible. The laser light is manipulated in devices known as interferometers designed for these measurements. Although precisely constructed, small defects in the components can greatly alter the measurement result. Such defects act to scatter the laser light from an ideal shape and so degrades the performance. Numerical models are used to construct virtual interferometers to analyse how these effects can be corrected or reduced in new instrument designs. However, the computational cost of calculating this scattering can severely limit the usefulness of the models. Our paper outlines a method that provides a thousand-fold improvement in the speed of this computation.

The dominant fraction of the computational cost comes from the evaluation of integrals describing the distortion of the laser beam. These are evaluated whenever the laser changes shape, for example, when optics heat up and expand. Our new, faster computational method can solve these integrals by making use of a technique known as reduced-order quadrature (ROQ). It generates a near-optimal application-specific integration rule for Hermite–Gauss beam modes (for a user-defined accuracy). The thousand-fold increase in performance demonstrated in this work (shown in the plot) opens up new areas of design explorations and optimisation which were practically off-limits before.

Authors: Brown, D., Smith, R. & Freise, A.
Journal: JOpt
arXiv: 1507.03806 [physics.optics]
Code: FINESSE

Astrophysical constraints on massive black hole binary evolution from pulsar timing arrays

Comparison of results to model predictions At the centre of many galaxies there is a massive black hole, millions to billions of times the mass of our Sun. When galaxies merge with each other, their black holes will form a binary and emit gravitational waves. Galaxies are merging all over the Universe, so we expect to see many overlapping signals forming a gravitational-wave background at nanohertz frequencies. How can we search for this background?

Pulsars are the lighthouses of the Universe. They are spinning neutron stars from which we see a pulse of radiation each rotation, and can be used as accurate clocks. Gravitational waves change the arrival times of the pulses from a set of pulsars in a characteristic way. Many international groups are searching for this evidence. So far no signal has been seen, but as pulsars are timed for longer and longer, our sensitivity improves. As a signal has not been found yet, we know that the gravitational-wave background must be smaller than the level we could detect and this sets an upper limit on the strength of the background.

But what can we learn about the massive black-hole binaries themselves using pulsar timing? Are there lots of large black holes merging far away from Earth or a few (relatively) small black holes merging close by? We use a simple model which can describe a plethora of possibilities for the distribution of binaries in mass and redshift to find out what we can learn about the population. Using the best result from pulsar timing, and some mock results for possible future upper limits and detection, we investigate what information pulsar timing holds.

Our main results are shown in the figure which shows results for three different pulsar timing scenarios (the current upper limit, a possible future upper limit, and detection) and some model predictions. We show distributions of sources in chirp mass M and in redshift z. The grey bands in the plot indicate what can be learnt from pulsar timing alone—without relying on any particular model prediction (green). The wide bands in the plot mean that the distribution is poorly constrained—we are unable to tell the distribution of masses or the redshift of the sources as many different possibilities are contained within these bands (dashed lines). However, what we can learn is the number of mergers taking place. For an upper limit we can place a maximum cap on the number of mergers, and for a detection we will be able to place a range on the number because we see that these binaries do indeed merge!

Authors: Middleton, H., Del Pozzo, W., Farr, W. M., Sesana, A. & Vecchio, A.
Journal: MNRASL
arXiv: 1507.00992 [astro-ph.CO]

Enhancing the bandwidth of gravitational-wave detectors with unstable optomechanical filters

Technical detials of Micelson-interferometer configurations According to Einstein's general theory of relativity, gravitational waves change the physical distance between spatially separated objects as they passing by. Such an effect is directly proportional to the initial distance of these objects—the larger the distance between the objects before the gravitational wave, the larger the change. That is why current gravitational-wave detectors are large-scale laser interferometers with kilometre arm lengths. To further increase the effective arm length, we introduce optical cavities, where the light bounces back and forth many times, to the interferometer that Michelson first invented in late 19 century. A schematic is shown on the left. However, the light is not able to instantaneously follow the arm length change induced by gravitational waves because it travels at a finite speed. Particularly for high-frequency gravitational-wave signals, the arm length of the interferometer changes so fast that the light will not able to follow, which results in a reduced signal strength from the interferometer. This means that the detector has lower sensitivity to these high-frequency gravitational waves. In this paper, we discuss a new approach to recover the high-frequency sensitivity by embedding an additional optical component, called an unstable filter (also shown in the figure). Qualitatively, the unstable filter enhances the detector response to high-frequency signals, resulting in an overall enhancement of the detector sensitivity across a broad frequency band.

Authors: Miao, H., Ma, Y., Zhao, C. & Chen, Y.
Journal: PRL
arXiv: 1506.00117 [quant-ph]

European Pulsar Timing Array limits on an isotropic stochastic gravitational-wave background

Gravitational-wave limits from EPTA Millisecond pulsars are the most stable clocks in the Universe, and can be used to probe extremely tiny variations in the spacetime metric caused by to gravitational waves. The passage of a gravitational wave squeezes and stretches the space so that the distance between the Earth and the pulsar is changing in time, causing a variation in the time of arrival of the pulses. By cross-correlating timing data of a ensemble of pulsars (a pulsar timing array), gravitational waves in the nanohertz frequency regime emitted by supermassive black hole binaries or cosmic strings and exotic processes in the early Universe, can in principle be observed. The European Pulsar Timing Array is one of the three international collaborations searching for these waves, and in this paper we employ the most advanced data analysis techniques on high quality pulsar timing data to search for gravitaional waves.

The figure on the left shows the results of our search, expressed in terms of energy density in gravitational wave (the Omega parameter) versus frequency. The black dashed curve is the limit we can put on the gravitational-wave energy density at each frequency. The straight lines show how this converts into a limit on the stochastic background generated by different sources (relic waves from the early Universe, cosmic strings, and supermassive black hole binaries). For black hole binaries, the limit is starting to skim through the region predicted by cosmological models (the red area). By observing longe) and adding pulsars to the array, our sensitivity is bound to increase, and pulsar timing arrays are expected to soon detect gravitational waves from these cosmological monsters.

Authors: Lentati, L., Taylor, S. R., Mingarelli, C. M. F., Sesana, A., Sanidas, S. A., Vecchio, A., Caballero, R. N., Lee, K. J., van Haasteren, R., Babak, S., Bassa, C. G., Brem, P., Burgay, M., Champion, D. J., Cognard, I., Desvignes, G., Gair, J. R., Guillemot, L., Hessels, J. W. T., Janssen, G. H., Karuppusamy, R., Kramer, M., Lassus, A., Lazarus, P., Liu, K., Osłowski, S., Perrodin, D., Petiteau, A., Possenti, A., Purver, M. B., Rosado, P. A., Smits, R., Stappers, B., Theureau, G., Tiburzi, C. & Verbiest, J. P. W.
Journal: MNRAS
arXiv: 1504.03692 [astro-ph.CO]

Measuring intermediate-mass black-hole binaries with advanced gravitational wave detectors

Lower bound on mass measurement Current electromagnetic observations provide evidence for black holes a few times more massive than the Sun, as well as behemoths millions of times more massive than the Sun. Intermediate-mass black holes weighing in at hundreds of solar masses are suspected to exist, but are extremely difficult to find electromagnetically, and none have been confirmed to date. In this paper, we used some of the analysis techniques we have developed within the LIGO Scientifitc Collaboration to demonstrate that gravitational-wave observations could be used to accurately measure masses of black holes in merging binaries and probe the existence of intermediate-mass black holes. We found that (i) at total masses below ~200 solar masses, where the signal-to-noise-ratio is dominated by the inspiral portion of the signal, the chirp mass 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 masses and more extreme mass ratios, and (iv) spin cannot be accurately measured for non-spinning components. Most importantly, we found that for moderately loud signals from systems with non-spinning components, an analysis that allows for aligned spins confirms the existance of an intermediate-mass black hole (with mass >100 solar masses) with 95% confidence in any binary that includes a component with a mass of 130 solar masses or greater (as illustrated in the figure on the left).

Authors: Veitch, J., Pürrer, M. & Mandel, I.
Journal: PRL
arXiv: 1503.05953 [astro-ph.HE]

Distinguishing compact binary population synthesis models using gravitational wave observations of coalescing binary black holes

Checkerborad plot Advanced LIGO is due to begin the hunt for gravitational waves from colliding black holes, with Advanced Virgo joining later next year. Black holes are the remains of stars with masses tens of times the mass of the Sun. For two of these black holes to collide, they need to live and die as a pair, or a binary. These stars lead spectacular lives as they constantly shed mass in stellar winds, violently strip mass from one another, explode as supernova, and even engulf one another (what is known as a common envelope event). Most of these processes have so far been studied using systems we can see, such as X-ray binaries (a black hole eating a normal star), but remain poorly understood. However, these processes also leave an imprint on the masses of the black holes. Therefore, by comparing the distribution of masses of binary black holes that Advanced LIGO and Advanced Virgo find with a prediction from a model, we can learn about these extreme evolutionary processes and the lives of massive stars. The chequerboard type plot shown on the left compares simulated observations from each model in a set of models to all of the others in turn. Darker squares indicate that models are similar, and lighter squares indicate that we can tell models apart. The dark strip on the diagonal shows that models do look like themselves. Since there are lots of lighter squares, it is possible to distinguish between many models. We should be able to learn lots about binary evolution over the next couple of years!

Authors: Stevenson S., Ohme, F. & Fairhurst, S.
Journal: ApJ
arXiv: 1504.07802 [astro-ph.HE]

Double tidal disruptions in galactic nuclei

Lightcurve Stars that get too close to massive black holes get shredded by the strong gravitational fields in spectacular tidal disruption events. In this paper, we explored what happens to stars in a binary which approaches a massive black hole on a deeply penetrating trajectory where the binary can move directly into the region where individual stars are tidally disrupted. We found that the interaction with the massive black hole often leads to both stars being disrupted in sequence. Such events could produce light curves that are substantially different from those of the single disruptions, with possible features such as two peaks (a possible light curve is shown on the left). Tidal forces from the massive black hole 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, something that has been observed but is difficult to explain with existing models.

Authors: Mandel, I. & Levin, I.
Journal: ApJL
arXiv: 1504.02787 [astro-ph.SR]

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

Sky map The first observing run (O1) of Advanced LIGO is expected later in 2015; we must have our data analysis algorithms prepared in anticipation of this. To check the how well we would be able to measure things in O1, we created a set of binary neutron-star signals (which are expected to be the main source for Advanced LIGO and hid these in noise. This noise was created using real data from initial LIGO to capture all the interesting (and troublesome) features. Our fake set of data was analysed as real data would be, passing through the detection pipeline before detected events were analysed with parameter-estimation codes. We checked how well we could measure position on the sky (which is important for pointing telescopes in case the merger of the neutron stars is accompanied by an explosion), distance to the source and the chirp mass (the combination of component masses we're most sensitive to). Uncertainty in sky position is typically about 600 deg2 (an example sky map is shown on the left), which will make it difficult to cover all the area. However, our fastest code (BAYESTAR) can calculate this in under a minute, so it is possible to start looking straight away! The distance is also quite uncertain, with typical uncertainties being tens of percent, but the chirp mass is well measured, with an accuracy typically about 0.00025 of a solar mass. To work out the individual masses, we need to include the effects of spin, which slows down the analysis: we are currently working on these results.

Authors: Berry, C. P. L., Mandel, I., Middleton, H., Singer, L. P., Urban, A. L., Vecchio, A., Vitale, S., Cannon, K., Farr, B., Farr, W. M., Graff, P. B., Hanna, C., Haster, C.-J., Mohapatra, S., Pankow, C., Price, L. R., Sidery, T. & Veitch, J.
Journal: ApJ
arXiv: 1411.6934 [astro-ph.HE]
Data release: The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo

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

Probability distribution for paramters of a black-hole binary How do extract information from a gravitational-wave signal about its source? We are often interested in what we can learn about astrophysical systems from gravitational-wave observations. Over many years, we have worked as part of the LIGO Scientific Collaboration to develop a set of algorithms that can analyse gravitational-wave signals from compact binary coalescences (the inspiral and merger of binaries made up of neutron stars or black holes), and infer the parameters that describe the source (such as the components' masses and spins, the orientation of the binary and its location). This software library is known as LALInference.

Given a measured signal, we can calculate how well it is matched by a template gravitational wave calculated with a particular set of parameters. The better the template fits, the more probably that the true signal had those parameters. We also use Bayes' theorem to fold in any prior information we have about the parameters. We are interested in mapping the probability distribution. LALInference does this by randomly sampling parameter space (using techniques such as Markov-chain Monte Carlo and nested sampling), at each point matching a template to the data. Doing this allows us to find the most probable parameter values, estimate our uncertainty and build up a complete picture of possible parameter combinations; an example probability distribution for the mass paramters of a black-hole binary is shown on the left. This process is computationally expensive, so LALInference features many techniques to help speed things up: we are still working on improvements here! Hopefully, Advanced LIGO and Advanced Virgo will soon detect a gravitational wave from a compact binary coalescence, LALInference will then tell us what we can learn about the source system.

Authors: Veitch, J., Raymond, V., Farr, B., Farr, W., Graff, P., Vitale, S., Aylott, B., Blackburn, K., Christensen, N., Coughlin, M., Del Pozzo, W., Feroz, F., Gair, J., Haster, C.-J., Kalogera, V., Littenberg, T., Mandel, I., O'Shaughnessy, R., Pitkin, M., Rodriguez, C., Röver, C., Sidery, T., Smith, R., Van Der Sluys, M., Vecchio, A., Vousden, W. & Wade, L.
Journal: PRD
arXiv: 1409.7215 [gr-qc]
Code: LALSuite

Gravitational-wave sensitivity curves

Sensitivity-curve plot This article explains some of the conventions used in gravitational-wave astronomy to characterise the sensitivities of detectors and the amplitudes of signals. Several different conventions are commonly in use (which can lead to confusion), each with their own beneficial properties. (i) The characteristic strain (as shown in the figure to the left), is closely related to the signal-to-noise ratio (how loud a signal is). On a log–log plot, the signal-to-noise ratio is the area between the detector's sensitivity curve and the top of the source box. (ii) The power spectral density is a measure of the sensitivity of the detector. It gives a measure of the average noise in the detector. (iii) The energy density is directly related to the energy carried by a gravitational wave. This quantity is often used by cosmologists to measure the contents of the Universe. The paper explains how to convert between these different quantities, and gives approximate values for the sensitivities of a range of detectors as well as the amplitudes for a range of possible sources. Accompanying this article is a website where you can create your own sensitivity-curve plots.

Authors: Moore, C. J., Cole, R. H. & Berry, C. P. L.
Journal: CQG
arXiv: 1408.0740 [gr-qc]
Website: Gravitational-wave sensitivity curve plotter

Sensitivity of intracavity filtering schemes for detecting gravitational waves

Michelson interferometer Advanced LIGO is currently the most sensitive gravitational-wave detector in the world. We are now at a stage where the quantum noise, arising from vacuum fluctuations of the electromagnetic field, sets the sensitivity limit over most of the detection band. Can we bring the kilometre-scale detectors beyond the quantum limit of light while retaining their dual-recycled Michelson configuration (combining power-recycling and signal-recycling mirrors)?

We investigated a new configuration that introduces additional optical filters inside the signal-recycling cavity (between the main interferometer and the signal-recycling mirror), as shown on the left. Optical filters can coherently feed the signal back to the interferometer with the proper frequency-dependent phase, resulting in different quantum-noise behaviours depending on the filter-cavity parameters. By adding these intracavity filters, a gravitational-wave detector can be tuned from being a pure position meter, to being partially a speed meter. Additionally, our optimisation of the intracavity filtering scheme suggested that it is as good as frequency-dependent input filtering in reducing quantum noise over a large frequency band.

Authors: Wang, M., Miao, H., Freise, A. & Chen, Y.
Journal: PRD
arXiv: 1310.3120 [gr-qc]

Realistic polarising Sagnac topology with DC readout for the Einstein Telescope

Sagnac sensitivity curve

The configuration of the Michelson interferometer was initially chosen for gravitational-wave detectors as it naturally sees gravitational-wave signals, directly outputting displacement differences between two perpendicular arms. The Sagnac interferometer was not as appealing to the gravitational-wave community until it was demonstrated that a zero-area Sagnac could potentially reach the required sensitivity to measure gravitational waves. The Sagnac interferometer, functions as a speed meter, which has excellent quantum-noise behaviour at low frequencies with only a little susceptibility to optical losses. As illustrated in the sensitivity-curve plot on the left, the low-frequency quantum noise of a Sagnac interferometer can be below the standard quantum limit (SQL), which would be unachievable for a Michelson interferometer without using more advanced techniques like frequency-dependent squeezing! This advantage encourages considering the Sagnac interferometer as a practical alternative for future gravitational-wave detectors such as the Einstein Telescope.

Authors: Wang, M., Bond, C., Brown, D., Brückner, F., Carborne, L., Palmer, R. & Freise, A.
Journal: PRD
arXiv: 1303.5236 [physics.ins-det]