## Zel’dovich Pancake Day!

Posted in The Universe and Stuff with tags , , , on February 16, 2021 by telescoper

Today it’s Shrove Tuesday but unfortunately I forgot to buy shroves yesterday so will have to make do with pancakes instead, but not the usual kind. I’ve blogged before about the Zel’dovich Approximation (published in Zeldovich, Ya.B. 1970, A&A, 5, 84) but there’s no harm in describing this classic again. Here’s the first page of the original paper:

In a nutshell, this daringly simple approximation considers the evolution of particles in an expanding Universe from an early near-uniform state into the non-linear regime as a sort of ballistic, or kinematic, process. Imagine the matter particles are initial placed on a uniform grid, where they are labelled by Lagrangian coordinates $vec{q}$. Their (Eulerian) positions at some later time $t$ are taken to be

$vec{r}(vec(q),t) = a(t) vec{x}(vec{q},t) = a(t) left[ vec{q} + b(t) vec{s}(vec{q},t) right].$

Here the $vec{x}$ coordinates are comoving, i.e. scaled with the expansion of the Universe using the scale factor $a(t)$. The displacement $vec{s}(vec{q},t)$ between initial and final positions in comoving coordinates is taken to have the form

$vec{s}(vec{q},t)= vec{nabla} Phi_0 (vec{q})$

where $Phi_0$ is a kind of velocity potential (which is also in linear Newtonian theory proportional to the gravitational potential).If we’ve got the theory right then the gravitational potential field defined over the initial positions is a Gaussian random field. The function $b(t)$ is the growing mode of density perturbations in the linear theory of gravitational instability.

This all means that the particles just get a small initial kick from the uniform Lagrangian grid and their subsequent motion carries on in the same direction. The approximation predicts the formation of caustics in the final density field when particles from two or more different initial locations arrive at the same final location, a condition known as shell-crossing. The caustics are identified with the main elements we find in large-scale structure. Because the initial collapse is usually along one direction the dominant structures are known as pancakes (or, as Zel’dovich himself might have called them, blini…).

Here’s a picture of a simulation showing these structures from the classic paper of Davis, Efstathiou, Frenk & White (1985):

Despite its simplicity this approximation is known to perform extremely well at reproducing the morphology of the cosmic web, although it breaks down after shell-crossing has occurred. In reality, bound structures are formed whereas the Zel’dovich approximation simply predicts that particles sail straight through the caustic which consequently evaporates.

## The Zel’dovich Lens

Posted in The Universe and Stuff with tags , , , , on June 30, 2014 by telescoper

Back to the grind after an enjoyable week in Estonia I find myself with little time to blog, so here’s a cute graphic by way of  a postscript to the IAU Symposium on The Zel’dovich Universe. I’ve heard many times about this way of visualizing the Zel’dovich Approximation (published in Zeldovich, Ya.B. 1970, A&A, 5, 84) but this is by far the best graphical realization I have seen. Here’s the first page of the original paper:

In a nutshell, this daringly simple approximation considers the evolution of particles in an expanding Universe from an early near-uniform state into the non-linear regime as a sort of ballistic, or kinematic, process. Imagine the matter particles are initial placed on a uniform grid, where they are labelled by Lagrangian coordinates $\vec{q}$. Their (Eulerian) positions at some later time $t$ are taken to be

$\vec{r}(\vec(q),t) = a(t) \vec{x}(\vec{q},t) = a(t) \left[ \vec{q} + b(t) \vec{s}(\vec{q},t) \right].$

Here the $\vec{x}$ coordinates are comoving, i.e. scaled with the expansion of the Universe using the scale factor $a(t)$. The displacement $\vec{s}(\vec{q},t)$ between initial and final positions in comoving coordinates is taken to have the form

$\vec{s}(\vec{q},t)= \vec{\nabla} \Phi_0 (\vec{q})$

where $\Phi_0$ is a kind of velocity potential (which is also in linear Newtonian theory proportional to the gravitational potential).If we’ve got the theory right then the gravitational potential field defined over the initial positions is a Gaussian random field. The function $b(t)$ is the growing mode of density perturbations in the linear theory of gravitational instability.

This all means that the particles just get a small initial kick from the uniform Lagrangian grid and their subsequent motion carries on in the same direction. The approximation predicts the formation of caustics  in the final density field when particles from two or more different initial locations arrive at the same final location, a condition known as shell-crossing. The caustics are identified with the walls and filaments we find in large-scale structure.

Despite its simplicity this approximation is known to perform extremely well at reproducing the morphology of the cosmic web, although it breaks down after shell-crossing has occurred. In reality, bound structures are formed whereas the Zel’dovich approximation simply predicts that particles sail straight through the caustic which consequently evaporates.

Anyway the mapping described above can also be given an interpretation in terms of optics. Imagine a uniform illumination field (the initial particle distribution) incident upon a non-uniform surface (e.g. the surface of the water in a swimming pool). Time evolution is represented by greater depths within the pool.  The light pattern observed on the bottom of the pool (the final distribution) displays caustics with a very similar morphology to the Cosmic Web, except in two dimensions, obviously.

Here is a very short  but very nice video by Johan Hidding showing how this works:

In this context, the Zel’dovich approximation corresponds to the limit of geometrical optics. More accurate approximations can presumably be developed using analogies with physical optics, but this programme has only just begun.

## The Zel’dovich Universe – Day 4 Summary

Posted in History, The Universe and Stuff with tags , , , , , , , on June 27, 2014 by telescoper

And on the fourth day of this meeting about “The Zel’dovich Universe”  we were back to a full schedule (9am until 7.30pm) concentrating on further studies of the Cosmic Web. We started off with a discussion of the properties of large-scale structure at high redshift. As someone who’s old enough to remember the days when “high redshift” meant about z~0.1 the idea that we can now map the galaxy distribution at redshifts z~2. There are other measures of structure on these huge scales, such as the Lyman alpha forest, and we heard a bit about some of them too.

The second session was about “reconstructing” the Cosmic Web, although a more correct word have been “deconstructing”. The point about this session is that cosmology is basically a backwards subject. In other branches of experimental science we set the initial conditions for a system and then examine how it evolves. In cosmology we have to infer the initial conditions of the Universe from what we observe around us now. In other words, cosmology is an inverse problem on a grand scale.  In the context of the cosmic web, we want to infer the pattern of initial density and velocity fluctuations that gave rise to the present network of clusters, filaments and voids. Several talks about this emphasized how proper Bayesian methods have led to enormous progress in this field over the last few years.

All this progress has been accompanied by huge improvements in graphical visualisation techniques. Thirty years ago the state of the art in this field was represented by simple contour plots, such as this (usually called the Cosmic Chicken):

You can see how crude this representation is by comparing it with a similar plot from the modern era of precision cosmology:

Even better examples are provided by the following snapshot:

It’s nice to see a better, though still imperfect,  version of the chicken at the top right, though I found the graphic at the bottom right rather implausible; it must be difficult to skate at all with those things around your legs.

Here’s another picture I liked, despite the lack of chickens:

Incidentally, it’s the back of Alar Toomre‘s head you can see on the far right in this picture.

The afternoon was largely devoted to discussions of how the properties of individual galaxies are influenced by their local environment within the Cosmic Web. I usually think of galaxies as test particles (i.e. point masses) but they are interesting in their own right (to some people anyway). However, the World Cup intervened during the evening session and I skipped a couple of talks to watch Germany beat the USA in their final group match.

That’s all for now. Tonight we’ll have the conference dinner, which is apparently being held in the “House of Blackheads” on “Pikk Street”. Sounds like an interesting spot!

## The Zel’dovich Universe – Day 3 Summary

Posted in History, The Universe and Stuff with tags , , , , , , on June 26, 2014 by telescoper

Day Three of this meeting about “The Zel’dovich Universe” was slightly shorter than the previous two, in that it finished just after 17.00 rather than the usual 19.00 or later. That meant that we got out in time to settle down for a beer in time the World Cup football. I watched an excellent game between Nigeria and Argentina, which ended 3-2 to Argentina but could have been 7-7. I’ll use that as an excuse for writing a slightly shorter summary.

Anyway we began with a session on the Primordial Universe and Primordial Signatures led off by Alexei Starobinsky (although there is some controversy whether his name should end -y or -i). Starobinsky outlined the theory of cosmological perturbations from inflation with an emphasis on how it relates to some of Zel’dovich’s ideas on the subject. There was then a talk from Bruce Partridge about some of the results from Planck. I’ve mentioned already that this isn’t a typical cosmology conference, and this talk provided another unusual aspect in that there’s hardly been any discussion of the BICEP2 results here. When asked about at the end of his talk, Bruce replied (very sensibly) that we should all just be patient.

Next session after coffee was about cosmic voids, kicked off by Rien van de Weygaert with a talk entitled “Much Ado About Nothing”, which reminded me of the following quote from the play of the same name:

“He hath indeed better bettered expectation than you must expect of me to tell you how”

The existence of voids in the galaxy distribution is not unexpected given the presence of clusters and superclusters, but they are interesting in their own right as they display particular dynamical evolution and have important consequences on observations. In 1984, Vincent Icke proved the so-called “Bubble Theorem” which showed that an isolated underdensity tends to evolve to a spherical shape.Most cosmologists, including myself, therefore expected big voids to be round, which turns out to be wrong; the interaction of the perimeter of the void with its surroundings always plays an important role in determining the geometry. Another thing that sprang into my mind was a classic paper by Simon White (1979) with the abstract:

We derive and display relations which can be used to express many quantitative measures of clustering in terms of the hierarchy of correlation functions. The convergence rate and asymptotic behaviour of the integral series which usually result is explored as far as possible using the observed low-order galaxy correlation functions. On scales less than the expected nearest neighbour distance most clustering measures are influenced only by the lowest order correlation functions. On all larger scales their behaviour, in general, depends significantly on correlations of high order and cannot be approximated using the low-order functions. Bhavsar’s observed relation between density enhancement and the fraction of galaxies included in clusters is modelled and is shown to be only weakly dependent on high-order correlations over most of its range. The probability that a randomly placed region of given volume be empty is discussed as a particularly simple and appealing example of a statistic which is strongly influenced by correlations of all orders, and it is shown that this probability may obey a scaling law which will allow a test of the small-scale form of high-order correlations.

The emphasis is mine. It’s fascinating and somewhat paradoxical that we can learn a lot about the statistics of where the galaxies are fom the regions where galaxies are not.

Another thing worth mentioning was Paul Sutter’s discussion of a project on cosmic voids which is a fine example of open science. Check out the CosmicVoids website where you will find void catalogues, identification algorithms and a host of other stuff all freely available to anyone who wants to use them. This is the way forward.

After lunch we had a session on Cosmic Flows, with a variety of talks about using galaxy peculiar velocities to understand the dynamics of large-scale structure. This field was booming about twenty years ago but which has been to some extent been overtaken by other cosmological probes that offer greater precision; the biggest difficulty has been getting a sufficient number of sufficiently accurate direct (redshift-independent) distance measurements to do good statistics. It remains a difficult but important field, because it’s important to test our models with as many independent methods as possible.

I’ll end with a word about the first speaker of this session, the Gruber prize winner Marc Davis. He suffered a stroke a few years ago which has left him partly paralysed (down his right side). He has battled back from this with great courage, and even turned it to his advantage during his talk when he complained about how faint the laser pointer was and used his walking stick instead.

## The Zel’dovich Universe – Day 2 Summary

Posted in History, The Universe and Stuff with tags , , , on June 25, 2014 by telescoper

Day Two of this enjoyable meeting involved more talks about the cosmic web of large-scale structure of the Universe. I’m not going to attempt to summarize the whole day, but will just mention a couple of things that made me reflect a bit. Unfortunately that means I won’t be able to do more than merely mention some of the other fascinating things that came up, as phase-space flip-flops and one-dimensional Origami.

One was a very nice review by John Peacock in which he showed that a version of Moore’s law applies to galaxy redshift surveys; since the first measurement of the redshift of an extragalactic object by Slipher in 1912, the number of redshifts has doubled every 2-3 years ago. This exponential growth has been driven by improvements in technology, from photographic plates to electronic detectors and from single-object spectroscopy to multiplex technology and so on. At this rate by 2050 or so we should have redshifts for most galaxies in the observable Universe. Progress in cosmography has been remarkable indeed.

The term “Cosmic Web” may be a bit of a misnomer in fact, as a consensus may be emerging that in some sense it is more like a honeycomb. Thanks to a miracle of 3D printing, here is an example of what the large-scale structure of the Universe seems to look like:

One of the issues that emerged from the mix of theoretical and observational talks concerned the scale of cosmic homogeneity. Our standard cosmological model is based on the Cosmological Principle, which asserts that the Universe is, in a broad-brush sense, homogeneous (is the same in every place) and isotropic (looks the same in all directions). But the question that has troubled cosmologists for many years is what is meant by large scales? How broad does the broad brush have to be? A couple of presentations discussed the possibly worrying evidence for the presence of a local void, a large underdensity on scale of about 200 MPc which may influence our interpretation of cosmological results.

I blogged some time ago about that the idea that the Universe might have structure on all scales, as would be the case if it were described in terms of a fractal set characterized by a fractal dimension $D$. In a fractal set, the mean number of neighbours of a given galaxy within a spherical volume of radius $R$ is proportional to $R^D$. If galaxies are distributed uniformly (homogeneously) then $D = 3$, as the number of neighbours simply depends on the volume of the sphere, i.e. as $R^3$, and the average number-density of galaxies. A value of $D < 3$ indicates that the galaxies do not fill space in a homogeneous fashion: $D = 1$, for example, would indicate that galaxies were distributed in roughly linear structures (filaments); the mass of material distributed along a filament enclosed within a sphere grows linear with the radius of the sphere, i.e. as $R^1$, not as its volume; galaxies distributed in sheets would have $D=2$, and so on.

We know that $D \simeq 1.2$ on small scales (in cosmological terms, still several Megaparsecs), but the evidence for a turnover to $D=3$ has not been so strong, at least not until recently. It’s just just that measuring $D$ from a survey is actually rather tricky, but also that when we cosmologists adopt the Cosmological Principle we apply it not to the distribution of galaxies in space, but to space itself. We assume that space is homogeneous so that its geometry can be described by the Friedmann-Lemaitre-Robertson-Walker metric.

According to Einstein’s theory of general relativity, clumps in the matter distribution would cause distortions in the metric which are roughly related to fluctuations in the Newtonian gravitational potential $\delta\Phi$ by $\delta\Phi/c^2 \sim \left(\lambda/ct \right)^{2} \left(\delta \rho/\rho\right)$, give or take a factor of a few, so that a large fluctuation in the density of matter wouldn’t necessarily cause a large fluctuation of the metric unless it were on a scale $\lambda$ reasonably large relative to the cosmological horizon $\sim ct$. Galaxies correspond to a large $\delta \rho/\rho \sim 10^6$ but don’t violate the Cosmological Principle because they are too small in scale $\lambda$ to perturb the background metric significantly.

The discussion of a fractal universe is one I’m overdue to return to. In my previous post I left the story as it stood about 15 years ago, and there have been numerous developments since then, not all of them consistent with each other. I will do a full “Part 2” to that post eventually, but in the mean time I’ll just comment that current large surveys, such as those derived from the Sloan Digital Sky Survey, do seem to be consistent with a Universe that possesses the property of large-scale homogeneity. If that conclusion survives the next generation of even larger galaxy redshift surveys then it will come as an immense relief to cosmologists.

The reason for that is that the equations of general relativity are very hard to solve in cases where there isn’t a lot of symmetry; there are just too many equations to solve for a general solution to be obtained. If the cosmological principle applies, however, the equations simplify enormously (both in number and form) and we can get results we can work with on the back of an envelope. Small fluctuations about the smooth background solution can be handled (approximately but robustly) using a technique called perturbation theory. If the fluctuations are large, however, these methods don’t work. What we need to do instead is construct exact inhomogeneous model, and that is very very hard. It’s of course a different question as to why the Universe is so smooth on large scales, but as a working cosmologist the real importance of it being that way is that it makes our job so much easier than it would otherwise be.

PS. If anyone reading this either at the conference or elsewhere has any questions or issues they would like me to raise during the summary talk on Saturday please don’t hesitate to leave a comment below or via Twitter using the hashtag #IAU308.

## The Zel’dovich Universe – Day 1 Summary

Posted in Biographical, History, The Universe and Stuff with tags , on June 24, 2014 by telescoper

I’m up possibly bright but definitely early to get ready for day two of IAU Symposium No. 308 The Zel’dovich Universe. The weather was a bit iffy yesterday, with showers throughout the day, but that didn’t matter much in practice as I was indoors most of the day attending the talks. I have to deliver the conference summary on Saturday afternoon so I feel I should make an effort to attend as much as I can in order to help me pretend that I didn’t write my concluding talk in advance of the conference.

Day One began with some reflections on the work and personality of the great Zel’dovich by two of his former students, Sergei Shandarin and Varun Sahni, both of whom I’ve worked with in the past.
Zel’dovich (left) was born on March 8th 1914. To us cosmologists Zel’dovich is best known for his work on the large-scale structure of the Universe, but he only started to work on that subject relatively late in his career during the 1960s. He in fact began his life in research as a physical chemist and arguably his greatest contribution to science was that he developed the first completely physically based theory of flame propagation (together with Frank-Kamenetskii). No doubt he also used insights gained from this work, together with his studies of detonation and shock waves, in the Soviet nuclear bomb programme in which he was a central figure, and which no doubt led to the chestful of medals he’s wearing in the photograph. In fact he was awarded the title of  Hero of Socialist Labour no less than three times.

My own connection with Zel’dovich is primarily through his scientific descendants, principally his former student Sergei Shandarin, who has a faculty position at the University of Kansas, but his work has had a very strong influence on my scientific career. For example, I visited Kansas back in 1992 and worked on a project with Sergei and Adrian Melott which led to a paper published in 1993, the abstract of which makes it clear the debt it owed to the work of Ze’dovich.

The accuracy of various analytic approximations for following the evolution of cosmological density fluctuations into the nonlinear regime is investigated. The Zel’dovich approximation is found to be consistently the best approximation scheme. It is extremely accurate for power spectra characterized by n = -1 or less; when the approximation is ‘enhanced’ by truncating highly nonlinear Fourier modes the approximation is excellent even for n = +1. The performance of linear theory is less spectrum-dependent, but this approximation is less accurate than the Zel’dovich one for all cases because of the failure to treat dynamics. The lognormal approximation generally provides a very poor fit to the spatial pattern.

The Zel’dovich Approximation referred to in this abstract is based on an extremely simple idea but which, as we showed in the above paper, turns out to be extremely accurate at reproducing the morphology of the “cosmic web” of large-scale structure.

Zel’dovich passed away in 1987. I was a graduate student at that time and had never had the opportunity to meet him. If I had done so I’m sure I would have found him fascinating and intimidating in equal measure, as I admired his work enormously as did everyone I knew in the field of cosmology. Anyway, a couple of years after his death a review paper written by himself and Sergei Shandarin was published, along with the note:

The Russian version of this review was finished in the summer of 1987. By the tragic death of Ya. B.Zeldovich on December 2, 1987, about four-fifths of the paper had been translated into English. Professor Zeldovich would have been 75 years old on March 8, 1989 and was vivid and creative until his last day. The theory of the structure of the universe was one of his favorite subjects, to which he made many note-worthy contributions over the last 20 years.

As one does if one is vain I looked down the reference list to see if any of my papers were cited. I’d only published one paper before Zel’dovich died so my hopes weren’t high. As it happens, though, my very first paper (Coles 1986) was there in the list. That’s still the proudest moment of my life!

We then went into a Dick Bond Special, with a talk entitled: From Superweb Simplicity to Complex Intermittency in the Cosmic Web. The following pic will give you a flavour:

It’s all very straightforward, really. Um…

The rest of the day consisted of a number of talks about the Cosmic Web of large-scale structure using techniques inspired by the work of Zel’dovich, particularly the Zel’dovich approximation which I’ve mentioned already. There were many fascinating talks but I had to single out Johan Hidding of Groningen for the best use of graphics. Here’s a video of his from Youtube as an example:

Well, I must get going for the start of Day Two. The first session starts at 9am (7am UK time) and the day ends at 19.30. Conferences like this are hard work!

PS. If anyone reading this either at the conference or elsewhere has any questions or issues they would like me to raise during the summary talk on Saturday please don’t hesitate to leave a comment below or via Twitter using the hashtag #IAU308.

## One Hundred Years of Zel’dovich

Posted in The Universe and Stuff with tags , , , , on March 12, 2014 by telescoper

Lovely weather today, but it’s also been an extremely busy day with meetings and teachings. I did realize yesterday however that I had forgotten to mark a very important centenary at the weekend. If I hadn’t been such a slacker that I took last Saturday off work I would probably have been reminded…

The great Russian physicist Yakov Borisovich Zel’dovich (left) was born on March 8th 1914, so had he lived he would have been 100 years old last Saturday. To us cosmologists Zel’dovich  is best known for his work on the large-scale structure of the Universe, but he only started to work on that subject relatively late in his career during the 1960s.  He in fact began his life in research as a physical chemist and arguably his greatest contribution to science was that he developed the first completely physically based theory of flame propagation (together with Frank-Kamenetskii). No doubt he also used insights gained from this work, together with his studies of detonation and shock waves, in the Soviet nuclear bomb programme in which he was a central figure, and which no doubt led to the chestful of medals he’s wearing in the photograph.

My own connection with Zel’dovich is primarily through his scientific descendants, principally his former student Sergei Shandarin, who has a faculty position at the University of Kansas. For example, I visited Kansas back in 1992 and worked on a project with Sergei and Adrian Melott which led to a paper published in 1993, the abstract of which makes it clear the debt it owed to the work of Ze’dovich.

The accuracy of various analytic approximations for following the evolution of cosmological density fluctuations into the nonlinear regime is investigated. The Zel’dovich approximation is found to be consistently the best approximation scheme. It is extremely accurate for power spectra characterized by n = -1 or less; when the approximation is ‘enhanced’ by truncating highly nonlinear Fourier modes the approximation is excellent even for n = +1. The performance of linear theory is less spectrum-dependent, but this approximation is less accurate than the Zel’dovich one for all cases because of the failure to treat dynamics. The lognormal approximation generally provides a very poor fit to the spatial pattern.

The Zel’dovich Approximation referred to in this abstract is based on an extremely simple idea but which, as we showed in the above paper, turns out to be extremely accurate at reproducing the morphology of the “cosmic web” of large-scale structure.

Zel’dovich passed away in 1987. I was a graduate student at that time and had never had the opportunity to meet him. If I had done so I’m sure I would have found him fascinating and intimidating in equal measure, as I admired his work enormously as did everyone I knew in the field of cosmology.  Anyway, a couple of years after his death a review paper written by himself and Sergei Shandarin was published, along with the note:

The Russian version of this review was finished in the summer of 1987. By the tragic death of Ya. B.Zeldovich on December 2, 1987, about four-fifths of the paper had been translated into English. Professor Zeldovich would have been 75 years old on March 8, 1989 and was vivid and creative until his last day. The theory of the structure of the universe was one of his favorite subjects, to which he made many note-worthy contributions over the last 20 years.

As one does if one is vain I looked down the reference list to see if any of my papers were cited. I’d only published one paper before Zel’dovich died so my hopes weren’t high. As it happens, though, my very first paper (Coles 1986) was there in the list. That’s still the proudest moment of my life!

Anyway, this post gives me the opportunity to advertise that there is a special meeting called The Zel’dovich Universe coming up this summer in Tallinn, Estonia. It looks a really interesting conference and I really hope I can find the time to fit it into my schedule. I’ve never been to Estonia…