Archive for Cosmology

Poisson (d’Avril) Point Processes

Posted in Uncategorized with tags , , , on April 2, 2019 by telescoper

I was very unimpressed by yesterday’s batch of April Fool jokes. Some of them were just too obvious:

I’m glad I didn’t try to do one.

Anyway, I noticed that an old post of mine was getting some traffic and when I investigated I found that some of the links to pictures were dead. So I’ve decided to refresh it and post again.


I’ve got a thing about randomness. For a start I don’t like the word, because it covers such a multitude of sins. People talk about there being randomness in nature when what they really mean is that they don’t know how to predict outcomes perfectly. That’s not quite the same thing as things being inherently unpredictable; statements about the nature of reality are ontological, whereas I think randomness is only a useful concept in an epistemological sense. It describes our lack of knowledge: just because we don’t know how to predict doesn’t mean that it can’t be predicted.

Nevertheless there are useful mathematical definitions of randomness and it is also (somtimes) useful to make mathematical models that display random behaviour in a well-defined sense, especially in situations where one has to take into account the effects of noise.

I thought it would be fun to illustrate one such model. In a point process, the random element is a “dot” that occurs at some location in time or space. Such processes occur in wide range of contexts: arrivals of buses at a bus stop, photons in a detector, darts on a dartboard, and so on.

Let us suppose that we think of such a process happening in time, although what follows can straightforwardly be generalised to things happening over an area (such a dartboard) or within some higher-dimensional region. It is also possible to invest the points with some other attributes; processes like this are sometimes called marked point processes, but I won’t discuss them here.

The “most” random way of constructing a simple point process is to assume that each event happens independently of every other event, and that there is a constant probability per unit time of an event happening. This type of process is called a Poisson process, after the French mathematician Siméon-Denis Poisson, who was born in 1781. He was one of the most creative and original physicists of all time: besides fundamental work on electrostatics and the theory of magnetism for which he is famous, he also built greatly upon Laplace’s work in probability theory. His principal result was to derive a formula giving the number of random events if the probability of each one is very low. The Poisson distribution, as it is now known and which I will come to shortly, is related to this original calculation; it was subsequently shown that this distribution amounts to a limiting of the binomial distribution. Just to add to the connections between probability theory and astronomy, it is worth mentioning that in 1833 Poisson wrote an important paper on the motion of the Moon.

In a finite interval of duration T the mean (or expected) number of events for a Poisson process will obviously just be proportional to the product of the rate per unit time and T itself; call this product λ.

The full distribution is then of the form:

This gives the probability that a finite interval contains exactly x events. It can be neatly derived from the binomial distribution by dividing the interval into a very large number of very tiny pieces, each one of which becomes a Bernoulli trial. The probability of success (i.e. of an event occurring) in each trial is extremely small, but the number of trials becomes extremely large in such a way that the mean number of successes is l. In this limit the binomial distribution takes the form of the above expression. The variance of this distribution is interesting: it is alsol.  This means that the typical fluctuations within the interval are of order the square root of l on a mean level of l, so the fractional variation is of the famous “one over root n” form that is a useful estimate of the expected variation in point processes.  Indeed, it’s a useful rule-of-thumb for estimating likely fluctuation levels in a host of statistical situations.

If football were a Poisson process with a mean number of goals per game of, say, 2 then would expect must games to have 2 plus or minus 1.4 (the square root of 2)  goals, i.e. between about 0.6 and 3.4. That is actually not far from what is observed and the distribution of goals per game in football matches is actually quite close to a Poisson distribution.

This idea can be straightforwardly extended to higher dimensional processes. If points are scattered over an area with a constant probability per unit area then the mean number in a finite area will also be some number l and the same formula applies.

As a matter of fact I first learned about the Poisson distribution when I was at school, doing A-level mathematics (which in those days actually included some mathematics). The example used by the teacher to illustrate this particular bit of probability theory was a two-dimensional one from biology. The skin of a fish was divided into little squares of equal area, and the number of parasites found in each square was counted. A histogram of these numbers accurately follows the Poisson form. For years I laboured under the delusion that it was given this name because it was something to do with fish, but then I never was very quick on the uptake.

This is all very well, but point processes are not always of this Poisson form. Points can be clustered, so that having one point at a given position increases the conditional probability of having others nearby. For example, galaxies like those shown in the nice picture are distributed throughout space in a clustered pattern that is very far from the Poisson form. But it’s very difficult to tell from just looking at the picture. What is needed is a rigorous statistical analysis.


The statistical description of clustered point patterns is a fascinating subject, because it makes contact with the way in which our eyes and brain perceive pattern. I’ve spent a large part of my research career trying to figure out efficient ways of quantifying pattern in an objective way and I can tell you it’s not easy, especially when the data are prone to systematic errors and glitches. I can only touch on the subject here, but to see what I am talking about look at the two patterns below:


You will have to take my word for it that one of these is a realization of a two-dimensional Poisson point process and the other contains correlations between the points. One therefore has a real pattern to it, and one is a realization of a completely unstructured random process.

I show this example in popular talks and get the audience to vote on which one is the random one. The vast majority usually think that the top  is the one that is random and the bottom one is the one with structure to it. It is not hard to see why. The top pattern is very smooth (what one would naively expect for a constant probability of finding a point at any position in the two-dimensional space) , whereas the bottom one seems to offer a profusion of linear, filamentary features and densely concentrated clusters.

In fact, it’s the bottom  picture that was generated by a Poisson process using a  Monte Carlo random number generator. All the structure that is visually apparent is imposed by our own sensory apparatus, which has evolved to be so good at discerning patterns that it finds them when they’re not even there!

The top  process is also generated by a Monte Carlo technique, but the algorithm is more complicated. In this case the presence of a point at some location suppresses the probability of having other points in the vicinity. Each event has a zone of avoidance around it; the points are therefore anticorrelated. The result of this is that the pattern is much smoother than a truly random process should be. In fact, this simulation has nothing to do with galaxy clustering really. The algorithm used to generate it was meant to mimic the behaviour of glow-worms which tend to eat each other if they get  too close. That’s why they spread themselves out in space more uniformly than in the random pattern.

Incidentally, I got both pictures from Stephen Jay Gould’s collection of essays Bully for Brontosaurus and used them, with appropriate credit and copyright permission, in my own book From Cosmos to Chaos. I forgot to say this in earlier versions of this post.

The tendency to find things that are not there is quite well known to astronomers. The constellations which we all recognize so easily are not physical associations of stars, but are just chance alignments on the sky of things at vastly different distances in space. That is not to say that they are random, but the pattern they form is not caused by direct correlations between the stars. Galaxies form real three-dimensional physical associations through their direct gravitational effect on one another.

People are actually pretty hopeless at understanding what “really” random processes look like, probably because the word random is used so often in very imprecise ways and they don’t know what it means in a specific context like this.  The point about random processes, even simpler ones like repeated tossing of a coin, is that coincidences happen much more frequently than one might suppose.

I suppose there is an evolutionary reason why our brains like to impose order on things in a general way. More specifically scientists often use perceived patterns in order to construct hypotheses. However these hypotheses must be tested objectively and often the initial impressions turn out to be figments of the imagination, like the canals on Mars.

Now, I think I’ll complain to wordpress about the widget that links pages to a “random blog post”. I’m sure it’s not really random….




Machine Learning in the Physical Sciences

Posted in The Universe and Stuff with tags , , , , , on March 29, 2019 by telescoper

If, like me, you feel a bit left behind by goings-on in the field of Machine Learning and how it impacts on physics then there’s now a very comprehensive review by Carleo et al on the arXiv.

Here is a picture from the paper, which I have included so that this post has a picture in it:

The abstract reads:

Machine learning encompasses a broad range of algorithms and modeling tools used for a vast array of data processing tasks, which has entered most scientific disciplines in recent years. We review in a selective way the recent research on the interface between machine learning and physical sciences.This includes conceptual developments in machine learning (ML) motivated by physical insights, applications of machine learning techniques to several domains in physics, and cross-fertilization between the two fields. After giving basic notion of machine learning methods and principles, we describe examples of how statistical physics is used to understand methods in ML. We then move to describe applications of ML methods in particle physics and cosmology, quantum many body physics, quantum computing, and chemical and material physics. We also highlight research and development into novel computing architectures aimed at accelerating ML. In each of the sections we describe recent successes as well as domain-specific methodology and challenges.

The next step after Machine Learning will of course be Machine Teaching…

BICEP2: Is the Signal Cosmological?

Posted in Astrohype, The Universe and Stuff with tags , , on March 28, 2019 by telescoper

An article in Physics Today just reminded me just now that I have missed the fifth anniversary of the BICEP2 announcement of `the detection of primordial gravitational waves’. I know I’m a week but I thought I’d reblog the post I wrote on March 19th 2014.You will see that I was sceptical…

..and it subsequently turned out that I was right to be so.

In the Dark

I have a short gap in my schedule today so I thought I would use it to post a short note about the BICEP2 results announced to great excitement on Monday.

There has been a great deal of coverage in the popular media about a “Spectacular Cosmic Discovery” and this is mirrored by excitement at a more technical level about the theoretical implications of the BICEP2 results. Having taken a bit of time out last night to go through the discovery paper, I think I should say that I think all this excitement is very premature. In that respect I agree with the result of my straw poll.

First of all let me make it clear that the BICEP2 experiment is absolutely superb. It was designed and built by top-class scientists and has clearly functioned brilliantly to improve its sensitivity so much that it has gone so…

View original post 1,015 more words

Fine-tuning in Cosmology

Posted in The Universe and Stuff with tags , , , , , on March 25, 2019 by telescoper

I forgot to post a link to a paper by Fred Adams that appeared on the arXiv last month on the topic of the fine-tuning of the Universe which I had bookmarked for a blog a while ago.

My heart always sinks when the arXiv informs me that the abstract of a paper is `abridged’ so here’s the full version from the PDF you can download for yourself here. Please be aware, though, that it’s a lengthy paper running to over two hundred pages:

My own view on this topic is that it is indeed remarkable that the Universe is finely-tuned to exactly the extent required to allow authors to write such long papers about the fine-tuning of the Universe…


The Most Ancient Heavens

Posted in Art, Biographical, Poetry, The Universe and Stuff with tags , , , , , , , , on March 21, 2019 by telescoper

So here I am, in that London, getting ready for the start of a two-day conference at the Royal Astronomical Society on cosmology, large-scale structure, and weak gravitational lensing, to celebrate the work of Professor Alan Heavens, on (or near) the occasion of his 60th birthday. Yes, it is a great name for an astronomer.

I was honoured to be invited to give a talk at this meeting, though my immediate reaction when I was told about was `But he can’t be sixty! He’s only a few years older than me…oh.’ I gather I’m supposed to say something funny after the conference dinner tomorrow night too.

Courtesy of alphabetical order it looks like I’m top of the bill!

Anyway, I’ve known Alan since I was a research student, i.e. over thirty years, and we’re co-authors on 13 papers (all of them since 2011).

Anyway, I’m looking forward to the HeavensFest not only for the scientific programme (which looks excellent) but also for the purpose of celebrating an old friend and colleague.

Just to clear up a couple of artistic points.

First, the title of the meeting, The Most Ancient Heavens, is taken from Ode to Duty by William Wordsworth.

Second, the image on the conference programme shown above is a pastiche of The Creation of Alan Adam which is part of the ceiling of the Sistine Chapel, waswhich painted by Michelangelo di Lodovico Buonarroti Simoni, known to his friends as Michelangelo. Apparently he worked flat out painting this enormous fresco. It was agony but the ecstasy kept him going. I’ve often wondered (a) who did the floor of the Sistine Chapel and (b) how could Michelangelo create such great art when it was so clearly extremely cold? Anyway, I think that is a picture of Alan at high redshift on the far right, next to the man with beard who at least had the good sense to wear a nightie to spare his embarrassment.

Anyway, that’s all for now. I must be going. Time for a stroll down to Piccadilly.

Update: you can find a bunch of pictures of this conference here.

New Publication at the Open Journal of Astrophysics!

Posted in Open Access, The Universe and Stuff with tags , , , , , , , , , on March 20, 2019 by telescoper

It’s nice to be able to announce that the Open Journal of Astrophysics has just published another paper. Here it is!

It’s by Darsh Kodwani, David Alonso and Pedro Ferreira from a combination of Oxford University and Cardiff University.

You can find the accepted version on the arXiv here. This version was accepted after modifications requested by the referee and editor.

This is another one for the `Cosmology and Nongalactic Astrophysics’ folder. We would be happy to get more submissions from other areas of astrophysics. Hint! Hint!

P.S. A few people have asked why the Open Journal of Astrophysics is not listed in the Directory of Open Access Journals. The answer to that is simple: to qualify for listing a journal must publish a minimum of five papers in a year. Since OJA underwent a failure long hiatus after publishing its first batch of papers we don’t yet qualify. However, so far in 2019 we have published four papers and have several others in the pipeline. We will reach the qualifying level soon and when we do I will put in the application!

A Bayesian Look at Cosmic Anomalies

Posted in Cosmic Anomalies with tags , , , on March 3, 2019 by telescoper

I’ve posted a few times on this blog about Cosmic Anomalies, by which I mean apparent departures from the predictions of the standard cosmological model. From time to time I also talk about this subject at seminars and conferences.

There’s an interesting new paper on this topic on the arXiv now by Shaikh et al., with the following abstract:

You can click on the image to make it larger. You can also find the PDF version of the full paper here.

I find this Bayesian analysis of two of the apparent anomalies (low amplitude in the power spectrum at large angular scales and hemispherical power asymmetry) may be different manifestations of the same underlying phenomenon, which would make them easier to account for without invoking new physics. Rather than being two independent statistical flukes these measurements might both be the result of one, which would be more likely to occur in the standard model. This analysis however suggests that this might not be the case after all, and these are two different things after all. This presupposes, however, that the model chosen to describe the asymmetries is appropriate. Anyway, this paper is well worth a read if you’re into Bayesian model testing (which you should be)…

This also gives me the excuse to post the following poll, which has been running for several years (even longer than Brexit):