Saturday, September 30, 2023

Artificial white water

 A few months ago I joined the Cambridge Canoe Club and went for my first paddle on a river in over a decade. But, that was on the rather flat River Cam. Since then, I've been out on the Cam a couple more times; once with other club members heading down river towards Grantchester and another time back up towards Jesus Lock on my own. Other than having to dodge low hanging branches and some ducks and swans, there's not been much of a challenge. However, over the summer the club was running some trips down to the Lee Valley White Water Centre, so I signed up to one to get some white water experience.

A week ago, I went to one of the club's pool sessions to try and see if I could remember how to roll. In my previous kayaking experience I'd become a dab hand at rolling (in a pool or on flat water at least - I don't think I ever successfully rolled on white water!), but it turned out that I needed a bit of practice to get the skill back. My first few attempts in the pool resulted in swims, but after some good coaching and patience from one of the other club members, I was able to roll again. This practice turned out to be useful for the Lee Valley course.

The Lee Valley White Water Centre had been built for the London 2012 Olympics, and I had actually been there before; not to paddle, but to watch some of the Olympic kayaking as shown in the picture below:


For the trip, we weren't going to be on the longer and harder Olympic course, but were instead on the shorter Legacy Loop course. After arriving and getting into my kit (I'd just bought a new Cag as the previous one I'd owned had suffered from some rubber decay around the neck) we got onto the small lake that you head across onto the Legacy Loop. As I wasn't quite sure where to go or what to do, I waited for some others to head over first. Once at the bottom of the course, I asked someone what I had to do and they suggested that I might want to get out and carry my boat halfway up the course and just do the easier bottom section to start with. On looking at the course, and feeling rather cocky as it looked nowhere near as intimidating as a Scottish river, I decided I'd just go straight to the top. To do this you had to paddle on to a short conveyor belt. Once at the top, it looked a bit more tricky, but I went for it and got myself down to the bottom in one piece - the main shock had been the first few times getting a face full of water after going over the drops, but it was fairly mild compared to what I've previously endured. I went straight back up for another go, but this time about halfway down I went over - I managed not to panic and decided that I may as well attempt to roll, so I put the practice from last week into action and lo-and-behold got myself back upright! The first time I'd actually rolled on white water!

I went again, capsized again, but also successfully rolled again getting a thumbs up from one of the more experienced club members who was on safety. I went a couple more times without any incident, but I was getting quite tired - I had the start of a cold coming on and hadn't actually felt great all day. After about an hour and a quarter (into the two hour session we had) I decided that I'd give it one more go before calling it a day. This time when I capsized I rushed my roll attempt a bit (I'll use my tiredness/cold as an excuse!) and failed, I tried once more, but then had to pop my deck and swim - my shoe came off and I didn't manage to hold on to either my boat or paddle, which headed most of the way down the course.

I couldn't end the session on a swim, so after emptying out the boat I went back for another go and managed to successfully complete another run. The swim had somehow invigorated me though and I decided to go again for another clean run. That didn't work out though and I swam again - I kept hold of my paddle and boat (and shoes) this time though, so was able to get back in and complete the course. This time I did decide that it was time to get out - especially as the last swim has caused me to get a big gash on my knuckle and bruises on both knees (mainly from trying to get out the water).

After getting dried and dressed (and luckily having plasters in my bag to dress my wound), I had the drive back up the M11 in the dark to look forward too. My overall experience at Lee Valley was a fun one (the facilities there are great) - it was great to get back on some white water and I look forward to doing it again, but hopefully not when nursing the start of a cold.

Wednesday, June 21, 2023

Punt slalom

According to this blog, the last time I went kayaking on a river was just over 11 years ago. Today I got back in a boat and did some actual river paddling for the first time since then.

For some additional context, last summer my family moved to Cambridge. Recently, my wife suggested to me that I should try out kayaking again, pointing me to the Cambridge Canoe Club website. Despite having just got rid of much of my kayaking gear, which had been sat mouldering in a plastic bag in lofts/garages over several house moves (my kag just had to be binned as the rubber seal around the neck had disintegrated), I decided that this was good idea. A few weeks ago I signed up for, and took, the club's induction session, which they required new members with some previous experience to do - they show you around the club, what boats and equipment they have, and then briefly take you out on the river (the Cam) to check you can paddle about at a reasonable proficiency and perform a wet exit (I attempted a roll, but didn't quite make it).

Today was my first proper paddle with the club, which was a nice sedate trip on the Cam from the club house (by Lammas land), up past The Backs and to Jesus lock. It didn't have quite the same thrill level as, say, taking on right angle on the Etive, but provided a good refresher on handling a boat on the water. It was a lovely evening, with the sun shining (I didn't have that too often on Scottish rivers in the autumn and winter) and the Cam was busy with people out on punts, so the main task was to avoid crashing into them.  Along the way I has to rescue a pigeon that for some reason was sat floating on the water, perhaps after getting itself confused and thinking it was one of the many ducks, swans or geese on the river. I successfully navigated around all these obstacles and manage to not capsize my boat. 

The trip was very pleasant and also provided a nice view into some of the Colleges that you can't get to other than via the river. However, I can already feel that my back is going to be protesting tomorrow about the exertions it had to perform. I expect I'll keep up the trips with the club and hopefully get in some white water paddling down at the Lee Valley White Water Centre over the summer.

Friday, June 08, 2018

Do millisecond pulsars have a minimum ellipticity?

tl;dr - The current millisecond pulsar population is suggestive that the answer is yes (with caveats), which would be very nice for the potential for observing gravitational wave from these pulsar with future gravitational wave observatories.

Millisecond pulsars (MSPs) are neutron stars, the extremely dense remnant cores of massive stars, that we observe to rotate very, very rapidly. In fact they spin with periods of a few milliseconds (meaning that they have surface velocities that are a few percent of the speed of light!), hence the name. They also have very strong surface magnetic fields of around 10,000 Tesla (the Earth's surface field is around 0.0001 T). So, to use the word again, they are pretty extreme. But, while they are characterised by some very large numbers, they also must be very smooth (and symmetric throughout the star). Or, to put it another slightly more technical way, they have a very small ellipticity, which is roughly speaking the difference in the radius of a the star along two perpendicular axes on a plane cutting through the star's equator, divided by it's polar radius. It's essentially the size of the biggest "mountain" or "bump" on a star compared to its total size. For the Earth the 'ellipticity' caused by Mount Everest (at 8.85 km) compared to the Earth's radius (6371 km) is ~0.001. By looking at how fast millisecond pulsars1 are observed to slow down (more on which is described below) we can work out that they can't have ellipticities greater than ~0.00000001, which translates to maximum mountains sizes of  ~0.1 mm given their radii of ~10 km.  So, we can infer that they can't have mountains above a certain size, but is there a minimum size to those mountains too?

Some colleagues (Graham Woan, Bryn Haskell, Ian Jones & Paul Lasky) and I have just released a paper looking to answer this question, which at least suggests that the answer might be yes. How this paper came about is interesting. In April myself (and the above colleagues) were all attending a workshop in Seattle at the Institute for Nuclear Theory. For this workshop Graham was putting together a talk titled "What can we learn from joint EM/GW observations of pulsars?" in which he produced a standard plot of the pulsar population often called a P-Pdot diagram. This is a plot that shows the base-10 logarithm of pulsar period derivatives (the rate at which their period is increasing, i.e., how quickly they are slowing down as they lose energy) against their period. (As a bit of a plug for myself) To create this he used a bit of software written by me, called PSRQpy, that downloads the latest pulsar information from the ATNF Pulsar Catalogue. On this plot you often see straight line contours running diagonally across the plot with slopes of -1, which represent the path that pulsars would evolve along over their lifetime if they had specific surface magnetic fields and lost energy primarily through magnetic dipole radiation. You can also plot similar lines for pulsars with certain ellipticities assuming they are losing energy primarily through the gravitational radiation cause by this ellipticity. These lines are a bit steeper with gradients of -3. When Graham plotted these ellipticity contours, in particular the line corresponding to an ellipticity of 10-9, he noticed that it looked like it provided a cut-off in the population of millisecond pulsars in the lower-left of the diagram (see slide 3 of Graham's talk). This was intriguing, but what could it mean? Well, it would suggest that there was some period-dependent process that dictated the minimum Pdot that a pulsar could have, and that that process might be gravitational wave emission from a common minimum ellipticity. But was this cut-off that we saw real (i.e. are there really no pulsars below the line), and if so what could cause it?

The first thing we needed to check to see whether the cut-off was real was that the period derivatives (Pdot) we were plotting were the true values. Observed Pdots are not the true, or intrinsic, Pdots of the pulsars. The observed values get contaminated by any relative motions and accelerations of the pulsar with respect to the Earth. In particular the proper motion, or relative velocity perpendicular to the line-of-sight, and differential galactic rotation between the Earth and pulsar, will alter the observed Pdot value. If you know the proper motion and distance to the pulsar, then you can correct for this to get the intrinsic Pdot values. So, this is what we did. The accelerations of pulsars in globular clusters can effect the observed Pdot value so much that it is completely swamped. In the majority of cases it is difficult to correct for this, so we ignore these globular cluster-based pulsars. But, as shown in the plot below, which is a close-up P-Pdot diagram for just the MSPs, a cut-off still seems to be present when using intrinsic Pdot values.

The period vs. period derivative diagram showing MSPs (with those in globular clusters removed). The small black stars show their observed period derivative, and the larger blue circles show their intrinsic period derivatives. Note the lack of pulsars in the lower-left part of the diagram. The shaded area on the lower-right part of the diagram is the "death line" below which neutron stars are expected to no longer be observable as pulsars.
Could this just be an observational selection effect, i.e., is there some reason why astronomers are biased against seeing pulsars in this region? Well that is a possibility, and by putting the paper out we are interested in hearing whether this is the case, but we could not think of any reasons why a selection effect causing such a cut-off would be present.

But, can we assign some significance to the reality of the cut-off and compare the evidence for different types of cut-off? If we make some assumptions we can. The details can be found in the paper, but we found that having a cut-off with a negative slope with a gradient2 of ~-3 (the best fit was around to at -3.6) was far better at explaining the observed P-Pdot diagram than having no cut-off at all. And a slope of -3 (which would suggest a process related to gravitational wave emission given a common minimum ellipticity) was ~35 times more likely than a a slope of -1 (which would suggest a process related to magnetic dipole emission and some common minimum surface magnetic field). If we assume a slope of -3 then the best fit cut-off position corresponds to an ellipticity of ~5×10-10 (this assumes a moment of inertia for the star of ~2×1038 kg m2, which is twice the value of the canonical moment of inertia used for the contour lines in the above plot, but means that the same gravitational wave amplitude would be produced as with an ellipticity of 10-9).

But, are there any known processes that mean MSPs might have a minimum ellipticity? Well, for quite a while it has been assumed that there was a minimum ellipticity cause by internal magnetic fields within the stars. But, this minimum value was very uncertain, depended on the internal field strength, and could have been several orders of magnitude below 10-9. However, we know that younger pulsars, which must be the progenitors of millisecond pulsars, have far larger (~1000 times) external surface magnetic field strengths. If similar field strengths were buried inside the stars during their accretion phase, and the interior of that stars are superconducting, then these could well be enough to give the required ellipticities. But, maybe we don't have to even invoke internal magnetic fields and superconductivity. Centrifugal forces mean that these rapidly rotating neutron stars are actually oblate spheroids (indeed their oblateness means they're 10s of metres wider around their middle than at their poles, which is a much large deformation that the ellipticities we have mentioned, but as the oblateness is spherically symmetric it does not induce gravitational waves). As MSPs are spun-up during accretion (or as they spin-down during their non-accreting lifetime) they will naturally want to assume a more (less) oblate shape, but the stiffness of their crusts will prevent them from changing shape... until it breaks. If the crust breaks un-evenly in order to reconfigure to the new shape, then this can provide a source of ellipticity. It has very recently been suggested that this uneven breaking is indeed occurring and leading to ellipticities of ~10-8, and in fact this is limiting the maximum speed at which MSPs are observed to rotate at.

So, what does this mean for gravitational waves and our prospects of observing them from MSPs. Given the assumption that all MSPs have ellipticities of 10-9, we worked out what their gravitational wave amplitude at Earth would be. We then worked out how strong we expected any signals from them to be in a variety of gravitational wave detectors over an observation period of a year: the Advanced LIGO and Advanced Virgo detectors operating at their design sensitivity (see, e.g. figure 1 of this paper); a pair of upgraded Advanced LIGO detectors called A+; a potential future "third generation" gravitational wave detector called the Einstein Telescope; and, another  potential future detector called Cosmic Explorer. The signal-to-noise ratios3 (SNR) for the different detector networks are shown in the figure below, where the filled histogram represents the assumption that all pulsars have ellipticities of 10-9 (the more optimistic un-filled histogram represents the SNRs you would expect if the pulsars were emitting such that all their slow down were due to gravitational wave emission). We see that for the advanced detectors, and their upgrades to A+, the SNRs would be rather small (< 5), which would most likely be rather tricky to be able to confidently detect. However, third generation detectors would be able to see many sources with SNRs greater than 10, which should be easily detectable with high confidence. If we're feeling more optimistic, and go with ellipticities of 10-9 being a minimum, but with some distribution above that, then we might expect SNRs to be somewhere between the filled and un-filled histograms in the figure, in which case the advanced detectors may be well placed to see signals from MSPs in the not to distant future!


So, if there really is a minimum ellipticity then things a looking fairly optimistic for a gravitational wave observations. Of course, our hypothesis would be quite easy to disprove if some MSPs are observed to be in the lower left-hand gap in the P-Pdot diagram.

1Non-millisecond, aka "young", pulsars are a different matter!
2The paper mainly talks about braking indices (see, e.g. Section 1.3 of this webpage) rather than gradients in the P-Pdot diagram. The gradient in the (base-10 logarithm) P-Pdot diagram is related to the braking index, n, via -(n-2). For magnetic dipole radiation n=3, while for gravitational radiation from a star with an ellipticity you would have n=5.
3These signal-to-noise ratios are the average you would expect when not knowing the orientations that the pulsars have with respect to the detectors. If the pulsars were optimally oriented these values would be about 1.7 times higher. They also assume a particular moment of inertia and distance to the stars (for the spin-down-based limits), which are themselves uncertain by factors of about two.

Tuesday, January 31, 2017

How high are those mountains now?

A few years ago I wrote about searches undertaken (by myself, others in the LIGO Scientific Collaboration, the Virgo Collaboration, and a selection of pulsar astronomers) looking for gravitational waves from a selection of pulsars using data from the initial science runs of the LIGO and Virgo detectors. Earlier this week we released the latest search for gravitational waves from 200 pulsars using data from LIGO's first observing run (O1) during the advanced detector era (after a complete upgrade of the detectors). The main result is that no compelling evidence for gravitational waves was found from any of the pulsars, but the results still provided some interesting highlights. A summary of the paper, and main results, can be found here, and I reproduce it below:

The recent detection of gravitational waves from merging binary black hole systems (GW150914 and GW151226) has now opened up the exciting new field of gravitational wave astronomy. However, the signals from these black hole binaries were transient and were only observable in our LIGO detectors for the last second, or less, before they merged. The sources were also more than a billion light years' distant, and so lay far beyond our own Milky Way Galaxy. The signals displayed a characteristic "chirping" form, in which the frequency, or pitch, and amplitude increased over the short duration of the signal. But, we are also searching for signals that could look quite different from that chirp pattern - for example sources that are constantly "on" (i.e. continuous) and emitting at an almost fixed (or monochromatic) frequency, like a pure audio tone.

One possible source of gravitational waves that could be much more local and would provide a continuous signal is a rapidly-rotating compact object known as a neutron star. We have known these objects exist ever since they were first observed as pulsars. They are the collapsed cores of massive stars that have run out of fuel and undergone a supernova explosion. With a mass of slightly more than the Sun's (∼2.8×1030 kg) packed into a sphere of radius ∼10 km, neutron stars are about 40 trillion times denser than lead. A teaspoon of neutron star material would weigh about 10 million metric tons, roughly equivalent to the weight of a small mountain on Earth. Neutron stars are also spinning very rapidly, and for some their surfaces are rotating at up to ∼10% of the speed of light. So, these stars are very extreme objects! Our best understanding is that their pulsed emission comes from beams of radiation emanating from the magnetic poles of the stars acting like a lighthouse. If the magnetic and rotation axes of the neutron star are not aligned then a pulse is observed as the radiation beam sweeps across the Earth once per rotation.

To generate gravitational waves a pulsar must have some non-symmetric distortion in its shape - e.g. a "mountain" - that is not along its rotation axis. There are various ideas about how such a distorted star could form. The distortion could be "frozen" into the crust or core of the star after it was born in the supernova, or formed from material falling onto the star, or be produced and maintained through extremely large internal magnetic fields. However, due to the huge gravitational pull at the star's surface, the material forming the "mountain" needs to be really strong so as not to be flattened out. A mountain on Earth made of jello, for example, could not get very big before collapsing under its own weight, but one made of solid rock can become as large as, or larger than, Everest. For a pulsar with a crust made up of "normal" neutron star material (highly distorted atomic nuclei, free electrons and neutrons) the maximum deformation that could be sustained before collapsing is about 10 cm, so not very high for a "mountain" (scaling up the relative dimensions this would be equivalent to a ∼50 m hill on Earth). If the star was made up from more exotic materials, e.g. if it were a solid quark star, then it could possibly sustain a "mountain" up to ∼10 m in height. The "mountain" size can also be expressed in terms of the star's ellipticity (ε), which is a rough measure of the distortion's size as a fraction of the star's radius.

Making a few reasonable assumptions we can estimate the maximum amplitude of gravitational waves being emitted by most pulsars. To do this we use the law of conservation of energy. Pulsars are seen to slow down (known as spin-down) over time. This spin-down takes a very long time, and even the most rapidly spinning-down objects only decrease in frequency by less than a hundredth of a Hertz (or equivalently, increase their periods by less than ten microseconds) over a year. But, given the huge moment of inertia that the stars possess, this still represents a very large loss in rotational energy, corresponding to a power of ∼1031 Watts, or well over ten thousand times the Sun's luminosity. If we assume that all of this energy is being lost by emission of gravitational waves we can calculate the amplitude with which we would observe those waves at the Earth. This is called the "spin-down limit". When our searches are sensitive enough to reach below this limit we start probing interesting new territory, where gravitational wave signals could be detectable. We do, however, know that the spin-down limit is a naive upper limit, in that a large part of the spin-down can also be attributed to other mechanisms, such as magnetic dipole radiation.

Just as with the "lighthouse" model for their electromagnetic emission, gravitational wave signals are expected at a frequency related to the rotation rate, typically at twice this value. There are just over 430 known pulsars (see the Australia Telescope National Facility pulsar catalogue) spinning fast enough for their gravitational wave emission to be in the sensitive frequency band of the current Advanced LIGO detectors (∼20 to 2000 Hz). Our previous searches in LIGO and Virgo gravitational wave observatory data looked for gravitational wave signals from 195 pulsars, with the spin-down limit being surpassed for two of them (the Crab and Vela pulsars).

In this new analysis we have searched for a total of 200 of these pulsars using data from the first observing run of the Advanced LIGO detectors. To help reach the best sensitivity we have used information about these pulsars obtained through radio and gamma-ray observations; these have provided very precise knowledge of the pulsars' positions, rotational frequencies, and how their frequencies change over time. This information has allowed us to accurately track any potential gravitational wave signal in our data over the whole length of the three-month science run (a search method called "coherent integration").

From these searches we were not able to detect evidence for gravitational radiation from any of the pulsars. But, we have produced the most sensitive upper limits yet, and for eight pulsars our observations have produced limits (using three largely independent statistical methods) on the gravitational wave amplitude that are below the spin-down limits. Two of the pulsars that have long been of interest for gravitational wave searches, due to their large spin-downs, are the Crab and Vela pulsars. These are the pulsars for which we have now surpassed the spin-down limit by a factor of 20 and 10 respectively. From this we can say that, respectively, less than ∼0.2% and 1% of their spin-down energy loss is due to gravitational radiation. We can also limit the ellipticity (roughly speaking the relative deformation, or "mountain" size, compared to the star's total size) of the stars, and say that there are no "mountains" on the Crab pulsar greater than ∼10 cm in height, and none on Vela greater than ∼50 cm. Among the other pulsars, we found 32 more that are within a factor of ten of the spin-down limit. From the gravitational wave data alone we can limit the "mountain" size for some of these to less than ∼0.1 mm, although the spin-down limit is more stringent for those pulsars.

Future science runs of the Advanced LIGO and Advanced Virgo gravitational wave detectors will provide even greater sensitivity for known pulsar searches, and, at the very least, should allow us to surpass the spin-down limits for ten-or-more sources. Searches are also underway for continuous signals that are not associated with any currently known pulsars. These have to search the whole sky and a broad range of frequencies and spin-down values, making the searches computationally intensive, but opening up the possibility of observing previously unknown objects.

Thursday, June 16, 2016

Yes, we've found another one

After the excitement of February when we (the LIGO Scientific Collaboration & Virgo Collaboration, or LVC for short) announced the first direct detection of a gravitational wave signal a lot of people having understandably been asking "Well, did you see any more?". The analysis performed for the announcement of the first detection (of the source called GW150914) used just over a month of data from the start of a longer observing run (that we called O1), which ran from the 12th September 2015 until the 19th January 2016. So we did have more data "in the can". And, as it happened that additional data did indeed provide us with another highly significant detection*. This new signal was observed in what was the early hours of Boxing Day in the UK, although it was still Christmas Day in the US when it hit the two LIGO detectors - so we can call it a late Christmas present. As we work in Coordinated Universal Time (UTC), which follows Greenwich Mean Time (GMT), the signal has been given the title GW151226 (i.e. it arrived on 26th December 2015), but internally has generally been called "The Boxing Day Event".

So, is this signal different from the first one? Well, like GW150914 it appears to be the result of two black holes inspiralling in to each other and merging, although the two merging black holes are smaller at roughly 14 and 8 times the mass of the Sun. A nice illustration of where these black holes sit, in terms of mass and radius, compared to other known black holes is shown here. However, unlike GW150914 we can pretty definitively say that at one of the merging black holes is spinning. Another thing to note is that if you look at Figure 1 from our detection paper for GW151226 you can't really see the signal in the data time series (whereas GW150914 stuck out like a sore thumb!) and it pretty much looks like noise.  As the system was less massive (but at a similar distance to) GW150914 the amplitude of the signal was intrinsically smaller. The saviour to this was that it also lasts longer in the detector's sensitive frequency bands (see Figure 1 in this paper that discusses all the detection's together) which means that you can integrate (basically sum together) over the longer signal and still "see" it in the noise.

Given that we'd already announced the detection you may be wondering what's important about this new one. The main thing is that we are now starting to reveal a population of objects rather than a single one. From looking at the population you can start to understand the distribution of source properties and investigate how they form. Admittedly with just two (and a bit) sources you really can't say much - it would be hard to work out the distribution of everyone's height by measuring just two people, but you at least get a rough idea of the likely range. It also allows us to be sure that the first signal wasn't a fluke, and suggests that we'll see many more of these objects in our upcoming observing runs (the next one, O2, should start this autumn with slightly better sensitivity than O1, and hopefully include the Virgo detector).

We often say that these gravitational waves are opening a new astronomical window on the Universe. And they really are! Imagine that the sky had always been covered in cloud, so you'd never been able to see the Moon, planets or stars (although in this scenario assume that you had a pretty good theory that the diffuse light coming through the clouds was being emitted by very distant objects called "stars".). Then, imagine that one night there's a slight chink in the clouds and through that you see a black sky with a single shining point of light in it. Wow! Your theory about "stars" was right! As the nights go on the clouds clear even more and reveal even more stars and other astronomical objects and the wonders of the Universe (and the exciting physics they reveal(!)) open up to you. It's a slightly tortured analogy, but you can kind of see that we're just seeing the first few points of light as the clouds are just starting to clear.

Some further information/reading:


*The eagle-eyed of you may have know that amongst all the papers produced about our detection announcement there was also mention of a candidate gravitational event that was within the originally analysed dataset. We've estimated that this candidate, dubbed LVT151012 (for LIGO-Virgo Trigger), has a roughly 90% chance of being a real astrophysical signal, but we like to be far more certain than that to claim it as a definite signal.

Thursday, February 11, 2016

The wait is over

Did you feel anything odd at around 09.50am GMT on 14th September 2015 (I'll let you do the time zone conversion)? Did you notice a disturbance in the force? Did you feel a tingle down your spine? Did you have butterflies in your stomach? Or, did you just feel a little bit wibbly?

No!? Well, at around that time a gravitational wave slammed into you at the speed of light, tried to rip your component atoms apart and then pull them together again1, and then passed out the other side of you. But you didn't even notice, did you! It's not particularly surprising you didn't feel anything as the disturbance the wave produced was spectacular mainly in its minuscule effect - the waves would have attempted1 to 'wibble' you from head-to-toe by only about 0.000000000000000000001 m (see Fig. 1 for an illustration of the effects of such a wave!).

Fig. 1. The effect of a passing gravitational wave on an alligator (as found around the LIGO Livingston observatory), a tumbleweed (as found in plentiful supply around the LIGO Hanford observatory) and myself (as found in the School of Physics & Astronomy at the University of Glasgow). Note that there is no discernible effect on any of these, except maybe as slight noticeable increase in my happiness at the prospect of the last 13 years of my working life having not been futile!
However, we (humanity in general, but large teams of scientists - many within the LIGO Scientific Collaboration [LSC] and Virgo Collaboration, including me - more specifically) have managed to build instruments that did indeed feel something on this date and time - a signal which we've given the catchy name GW150914 (or "The Event" as it was known for a while within the collaboration). These instruments, in this case the two US-based LIGO observatories (now entering their advanced phase), one in Hanford, Washington and the other in Livingston, Louisiana2, both felt the waves' passing and saw a very consistent signal (see Fig. 2) - other than it looking like exactly what we'd expect from a gravitational wave source, our confidence that this was a real signal came from empirically estimating how often such a consistent and strong signal would have be seen by chance (i.e., from random [generally non-Gaussian] noise fluctuations in the detectors), which we work out as being less than once per 200,000 years. So, we're pretty sure (greater than 5.1σ in annoyingly frequentist statistical terminology), i.e. certain, that the signal was real. And, what's more, we've been able to use the pattern of wibbles the instruments felt to work out that this gravitational wave was emitted by two black holes, both tens of times more massive than the Sun, whacking into each other at about half the speed of light, to form the largest small(!) (~solar mass) black hole we know of. The amount of energy this event emitted was a whopping 5×1047 Joules (quite possibly the most luminous event we've ever observed), equivalent to three times the mass of the Sun being converted directly to energy (remember E=mc2). Or, if you're into some "fun" energy conversions this is apparently equivalent to "the number of kilocalories in 2×1028 cubic kilometres of butter (that's the volume of 14 billion Suns of pure butter!)" or "3 quadrillion times the energy required to destroy the planet Earth"!

Fig. 2. The GW150914 gravitational wave signal observed in the two LIGO detectors (this is figure 1 from the discovery paper, Abbott et al., Phys. Rev. Lett., 116, 061102 (2016)) (Credit: The LIGO Scientific Collaboration and the Virgo Collaboration)
This is the first ever detection3 we've had and it is quite a big deal, both scientifically (there's a whole load of awesome astrophysics that has been done using the signal, and it opens up a whole new area of astronomy) and personally. Some people have been in the gravitational wave detection game for almost five decades, but we've all had to wait patiently without seeing any definite sign of them in our detectors4 (this has occasionally been to the amusement of colleagues in other areas of physics and astronomy). As a member of the LSC myself since starting studying for my PhD at the University of Glasgow in the autumn of 2002 (when what is now known as the "initial" LIGO detectors had just started taking data) I've only been waiting 13 years, but that's still my entire working life. For some additional context here's what I wrote in my thesis acknowledgements section back in 2005:
I was attracted to the field of gravitational wave research due to the promise that we would be entering exciting times with several large scale projects bringing in unprecedented amounts of new data. Given this the discovery of gravitational waves would be just round the corner, opening up gravitational wave astronomy for real. Little did I know that this has been exactly what’s been said for around 30 years! Despite this I do actually believe that I’ve entered the field at a prime time and finding gravitational waves is just around, if not the first corner, then the next one.
So, ten years later, we have now turned that "next corner" and gravitational wave astronomy is finally with us! Many more detections should now be forthcoming in our future observation runs, hopefully including other exciting sources as well as more merging pairs of black holes.

Finally, the other rather cool, and timely, thing is that the signal arrived a century after Einstein published his General Theory of Relativity from which the prediction of gravitational waves arises. Einstein's prediction of gravitational waves and also Schwarzschild's solution to Einstein's equations from which predictions of black holes would arise, were both published a century ago in 1916.

The paper describing the detection and analysis of GW150914 has been peer reviewed and is now published (Abbott et al, Phys. Rev. Lett., 116, 061102, 2016) and further papers detailing the detectors, analyses and science results can be found here. Also, summaries (at a less detailed level) of the main science we've obtained from the signal can be found here. The data containing the signal and some example codes showing how to view, and hear(!), it are available here, so you should go ahead and take a look yourself.

More information, reactions and opinions about this amazing discovery (and the time-line of how the detection happened5) can be found in blog posts by several fellow collaboration members linked below and within a special edition of the LIGO Magazine:
Also, be sure follow @ligo and @ego_virgo on twitter along with the hastags #GravitationalWaves, #EinsteinWasRight, #BinaryBlackHole and #AdvancedLIGO.

P.S. If you want to know what I was doing when GW150914 passed by it probably involved nappies (diapers for those in the US), feeding a child/cleaning bottles, or doing laundry, as I was on the final day of paternity leave following the birth of my second child (here's my first child one simulating a gravitational wave chirp). I didn't see the growing emails about the signal until early that afternoon when I thought I should try clearing out my inbox before returning to work the following day. It definitely made going back to work that bit more exciting. However, all the work to get the analyses of this event checked has slightly eaten into my main job, which is to search for gravitational waves from pulsars.

1 Don't worry, gravity is very weak. The forces keeping your various constituent bits and bobs together are far more than enough to overcome any piddling gravitational wave that passes through you. Any displacement (stretching or squeezing) is only noticeable (if you have an exquisitely sensitive gravitational wave detector at least) between freely falling objects in the same local frame, i.e. if there are effectively no other external forces acting on the objects.
2 Note that another detector called Virgo is also due to start taking data later this year, but wasn't operational at the time, and that a smaller (and unfortunately less sensitive, but still important) detector called GEO600 was operational, but not observing when the event passed by.
3 Prior to this direct detection (some quibbles over the "directness" of detection can be found here) many gravitational waves have obviously continuously been impinging on the Earth and passing through us, but this is the first time we've had the technology to catch one in the act.
4 There are very good reasons why they've not been seen until now (basically boiling down to us not having been able to build sensitive enough detectors), but that hasn't made us any less impatient.
5 The signal was first "detected" about three minutes after it arrived by online analysis software that looked for generic transient (short-duration) coherent signals, i.e. blips in the data that appeared at the same time in both detectors and looked similar. The first reasonably detailed estimates of the source parameters (i.e. that is was two black holes merging) were with us within about a day. After a couple of days we release our estimate of the location of the source in the sky and released it to selected astronomy groups to point their telescopes at. Following that full and proper detailed studies of the signal, and very careful checks on the performance of the detectors, have taken many months of painstaking work and placed great strain6 on the collaboration. But, given the general inertia that you get within a large collaboration (which we've experienced releasing results that didn't contain any signals) we've actually turned around a finished detection paper (after 14 draft iterations), and twelve(!) companion papers in remarkably short time.
6 This is a hilarious gravitational wave pun.

Disclaimer: everything on this blog is my own personal opinion and any mistakes are my own.

Monday, March 24, 2014

GWPAW 2013: Impressions from India

In the latest issue of the LIGO Magazine I have a short article on my (relatively) recent trip to India to attend the Gravitational Wave Physics & Astronomy Workshop. Below I reproduce (a partially un-edited version of [apologies to the editors for reverting some of their changes here]) the article, with added links!
Family constraints have meant I’ve been off the conference circuit for a bit, so the 3rd Gravitational Wave Physics & Astronomy Workshop (GWPAW, formerly the Gravitational Wave Data Analysis Workshop, GWDAW, which ran on 14 occasions) seemed like a good opportunity to get back into the swing of conference attendance. Plus, its location at the Inter-University Centre for Astronomy & Astrophysics (IUCAA) in Pune, India presented the chance to visit a new country. Due to the location of the meeting, many of the other non-local attendees were able to experience a bit of India, including a group that organised a tour round Mumbai (and subsequent train journey to Pune), a couple who started their trip with a holiday in the backwaters of Kerala, and others visiting family or friends. While it would have been a great opportunity for me to see India, I was unable to bookend my trip with any site-seeing, so my experience of India outside of the confines of IUCAA mainly came from my taxi ride from Mumbai to Pune. The taxi ride itself was an interesting insight into travel in India - the first half of the approximately three and a half hour ride (it’s about a 170km journey) was just in leaving Mumbai, where the roads that are about as chaotic as they come. The system seems to be to spot a gap in the traffic, even if it looks too small for the mode of transport you are in, and then squeeze into it. Astonishingly this method (accompanied by liberal application of the horn) got us through the traffic unscathed. The freeway between Mumbai and Pune is apparently one of the best roads in India, and can supposedly offer great views as you climb up into the rocky hills, but a combination of jet lag and low clouds/smog meant that I couldn’t appreciate the trip/views fully (from the plane on my flight back from Pune to Mumbai I was able to see the views I'd previously missed). 
In Pune I stayed at the very pleasant Seasons Apartment Hotel, which as the name suggests offered large apartments with a lounge and kitchenette (and free bottled water, which is a must for travellers there). Not feeling very adventurous on my arrival I just opted for dinner at the hotel, but it was definitely worthwhile as the open air rooftop bar/restaurant offered great views of the city. The hotel was just about within walking distance of IUCAA, where the meeting was held (which I had briefly considered as a travel option), but the organisers had put on a taxi service to and from the hotel every day. On travelling to IUCAA I was thankful for this as negotiating the roads, many of which lacked pavements, may have proved daunting. IUCAA itself is situated on the Pune University campus, but is fairly self-contained with its own “housing colony” for guests, students and postdocs to stay. During the meeting we didn’t have to go far between talks in the Chandrasekhar auditorium, coffee breaks (which consisted of strong black tea really) breaks and meals. 
As well as our taxi service the organisers provided breakfast, lunch and dinner within IUCAA under a large marquee. The food was great, although you may have been hard-pressed if you didn’t like curry - not a problem for me though! Some of the dishes were pretty spicy, but I suspect they were they were probably still toned down from their usual standard heat levels. We also had freshly made roti cooked in a tandoor oven by the side of the marquee. 
Kathak dance recital
On the first evening we had entertainment put on in the form of a Kathak Dance Recital in the meeting auditorium. The singing and musical accompaniment was mesmerising. Afterwards Sathya presented the dancers and musicians with houseplants, which I can only assume is the standard thank-you gift.
And what about the science? The meeting was weighted towards compact binary coalescences (CBC) and electromagnetic follow-up, but that’s not surprising given that these are the most likely sources of the first advanced detector observations. In fact it was good to have a GWPAW where many of talks were about things that could be done in the near future, rather than having to look ahead decades, further cementing the idea that gravitational wave detections are on the horizon! A couple of standout talks were Parameswaran Ajith's overview of the status and prospects for modelling CBC waveforms and Jocelyn Read’s talk on the potential for measuring neutron star equations of state with advanced detectors. Most sessions had lively discussions following the talks, with one particular participant always ready to provide some vigorous questioning. 
The breaks and poster sessions in the grounds of the auditorium (which amongst other things contained a giant sundial and a set of swings connected as a coupled harmonic oscillator) were always buzzing with conversation, which for me yielded a potential future collaboration with an IUCAA postdoc. There were many interesting posters, but I particularly liked a couple: one was Chris Messenger’s describing a method to extract redshift information from neutron star mergers by observing modes of a potentially short-lived post-merger hyper-massive neutron star; and another was Shaon Ghosh’s on electromagnetic follow-up of CBC signals. During the meeting my own poster was upgraded to a talk (due to passport related issues for one of the invited speakers causing him to miss the meeting), so I had to quickly put together my own slides. 
The meeting turned out to be incredibly productive and fascinating, as well as welcoming and well-organised. The organisers and IUCAA staff were really friendly and helpful. It was a great chance for many Indian students and postdocs to attend the meeting and share their work, and for people from the LVC to interact with them. This was particularly useful because the distance means many collaborators in the USA and Europe got to discuss topics in person, and allowed us to develop these relationships in the run-up to LIGO India. This will be good for bringing through new local people into the field in the run up to LIGO India. There was a great deal of enthusiasm from the IUCAA director Ajit Kembhavi to keep up the efforts with the suggestion that IUCAA and other Indian institutions host summer school-type events in the future. The next GWPAW to look forward to will be in Osaka, Japan in June 2015, closely followed by Amaldi in South Korea. 
It’s a shame I didn’t get to experience more of the country, but I did I get to discover a taste for the Indian Coca-Cola equivalent, “Thums-Up”, while discussing exciting science halfway around the world.