Archive for Power spectrum

Two New Publications at the Open Journal of Astrophysics

Posted in OJAp Papers, Open Access, The Universe and Stuff with tags , , , , , , , , , , , on February 1, 2024 by telescoper

As the first month of 2024 is now over, I thought I’d post an update relating to the  Open Journal of Astrophysics.  Since the last update we have published two papers, taking  the count in Volume 7 (2024) up to 9 (the total for January) and the total published by OJAp up to 124. We will have others soon, but I will be travelling for the first few days of February so the next update will be in a week or so.

Using our sophisticated forecasting algorithm, based on the first month of 2024 as input, I predict that we will publish around 9×12=108 papers in 2024, more than double last year’s total of 50.

Both the current papers discussed here are in the folder marked Cosmology and NonGalactic Astrophysics, our most popular category.

Anyway, the first paper of the most recent pair – published on January 30th – is “Capse.jl: efficient and auto-differentiable CMB power spectra emulation”, by Marco Bonici (INAF Milano, Italy & Waterloo, Canada), Federico Bianchini (Stanford, USA) and Jaime Ruiz-Zapatero (Oxford, UK). This paper presents an emulator for rapid and accurate prediction of Cosmic Microwave Background (CMB) temperature, polarization, and lensing angular power spectra, that works much faster than traditional methods. The code is written in Julia, in which language we are seeing an increasing number of submissions to OJAp.

Here is a screen grab of the overlay which includes the abstract:

 

You can click on the image of the overlay to make it larger should you wish to do so. You can find the officially accepted version of the paper on the arXiv here.

The second paper was published yesterday (31st January 2024) and has the title “Cosmological Inflation in N-Dimensional Gaussian Random Fields with Algorithmic Data Compression” which is a  study of inflationary models with Gaussian random potentials for multiple scalar fields, tracking the evolutionary trajectories numerically. The authors are Connor Painter and Emory Bunn, both the Physics Department at the University of Richmond, Virginia (USA). Ted Bunn (as he is usually known) is a longstanding member of the Editorial Board of the Open Journal of Astrophysics (and was thereby excluded from any involvement in the editorial process for this paper).

Here is a screen grab of the overlay which includes the abstract:

 

 

 

You can click on the image of the overlay to make it larger should you wish to do so. You can find the officially accepted version of the paper on the arXiv here.

Here Endeth the Update.

Leave a comment »

With the Cosmic Web in Mind..

Posted in Astronomy Lookalikes, The Universe and Stuff with tags , , , , , on November 23, 2019 by telescoper

Some time ago I posted one of my Astronomy Look-alikes about the remarkable similarity between the structure of the human brain and that revealed by computer simulations of the large-scale structure of the Universe:

I wonder whether this means that the Cosmic Web is really just all in the mind?

Anyway I just came across an article by Franco Vazza and Alberto Fenetti that takes the comparison between brain cells (among other things) and the Cosmic Web a bit further, including a look at the corresponding power spectra:

The main point to take from this picture is that many naturally occurring patterns have approximately power-law power spectra, at least over a limited range of scales. However, as I have pointed out before on this blog, the power spectrum on its own does not really quantify pattern in any meaningful way. Here for example are two patterns with exactly the same power spectrum:

The point is that the power spectrum does not contain any information about the phase correlations of the Fourier modes, which are important in quantifying localised features. For further discussion of this issue, see here.

2 Comments »

The Power Spectrum and the Cosmic Web

Posted in Bad Statistics, The Universe and Stuff with tags , , , , , , on June 24, 2014 by telescoper

One of the things that makes this conference different from most cosmology meetings is that it is focussing on the large-scale structure of the Universe in itself as a topic rather a source of statistical information about, e.g. cosmological parameters. This means that we’ve been hearing about a set of statistical methods that is somewhat different from those usually used in the field (which are primarily based on second-order quantities).

One of the challenges cosmologists face is how to quantify the patterns we see in galaxy redshift surveys. In the relatively recent past the small size of the available data sets meant that only relatively crude descriptors could be used; anything sophisticated would be rendered useless by noise. For that reason, statistical analysis of galaxy clustering tended to be limited to the measurement of autocorrelation functions, usually constructed in Fourier space in the form of power spectra; you can find a nice review here.

Because it is so robust and contains a great deal of important information, the power spectrum has become ubiquitous in cosmology. But I think it’s important to realise its limitations.

Take a look at these two N-body computer simulations of large-scale structure:

The one on the left is a proper simulation of the “cosmic web” which is at least qualitatively realistic, in that in contains filaments, clusters and voids pretty much like what is observed in galaxy surveys.

To make the picture on the right I first  took the Fourier transform of the original  simulation. This approach follows the best advice I ever got from my thesis supervisor: “if you can’t think of anything else to do, try Fourier-transforming everything.”

Anyway each Fourier mode is complex and can therefore be characterized by an amplitude and a phase (the modulus and argument of the complex quantity). What I did next was to randomly reshuffle all the phases while leaving the amplitudes alone. I then performed the inverse Fourier transform to construct the image shown on the right.

What this procedure does is to produce a new image which has exactly the same power spectrum as the first. You might be surprised by how little the pattern on the right resembles that on the left, given that they share this property; the distribution on the right is much fuzzier. In fact, the sharply delineated features  are produced by mode-mode correlations and are therefore not well described by the power spectrum, which involves only the amplitude of each separate mode. In effect, the power spectrum is insensitive to the part of the Fourier description of the pattern that is responsible for delineating the cosmic web.

If you’re confused by this, consider the Fourier transforms of (a) white noise and (b) a Dirac delta-function. Both produce flat power-spectra, but they look very different in real space because in (b) all the Fourier modes are correlated in such away that they are in phase at the one location where the pattern is not zero; everywhere else they interfere destructively. In (a) the phases are distributed randomly.

The moral of this is that there is much more to the pattern of galaxy clustering than meets the power spectrum…

26 Comments »

SPT and the CMB

Posted in The Universe and Stuff with tags , , , , , , , , , on November 30, 2012 by telescoper

I’ve been remiss in not yet passing on news  from the South Pole Telescope, which has recently produced a number of breakthrough scientific results, including:  improved cosmological constraints from the SPT-SZ cluster survey (preprint here); a new catalogue of 224 SZ-selected cluster candidates from the first 720 square-degrees of the survey (preprint here); the first measurement of galaxy bias from the gravitational lensing of the CMB (preprint here); the first CMB-based constraint on the evolution of the ionized fraction during the epoch of reionization (preprint here); the most-significant detection of non-Gaussianity induced from the gravitational lensing of the CMB (preprint here); and the most precise measurement of the CMB damping tail and improved constraints on models of Inflation (preprint here).

Here’s the graph that drew my eye (from this paper). It shows the (angular) power spectrum of the cosmic microwave for very high (angular) frequency spherical harmonics; the resolution of SPT allows it to probe finer details of the spectrum that WMAP (also shown, at lower l).

Slide1

This is an amazing graph, especially for oldies like me who remember being so impressed by the emergence of the first “acoustic peak” at around l=200 way back in the days of Boomerang and Maxima and gobsmacked by WMAP’s revelation of the second and third. Now there are at least six acoustic peaks, although of progressively lower amplitude. The attenuation of the CMB fluctuations at high frequencies is the result of diffusion damping – similar to the way high-frequency sound waves are attenuated when they pass through a diffusive medium (e.g. a gas).  The phenomenon in this case is usually called Silk Damping, as it was first worked out back in the 1960s by Joe Damping Silk.

Anyway, there’ll be a lot more CMB news early (?) next year from Planck which will demonstrate yet again that cosmic microwave background physics has certainly come a long way from pigeon shit…

2 Comments »

Power versus Pattern

Posted in Bad Statistics, The Universe and Stuff with tags , , , , , on June 15, 2012 by telescoper

One of the challenges we cosmologists face is how to quantify the patterns we see in galaxy redshift surveys. In the relatively recent past the small size of the available data sets meant that only relatively crude descriptors could be used; anything sophisticated would be rendered useless by noise. For that reason, statistical analysis of galaxy clustering tended to be limited to the measurement of autocorrelation functions, usually constructed in Fourier space in the form of power spectra; you can find a nice review here.

Because it is so robust and contains a great deal of important information, the power spectrum has become ubiquitous in cosmology. But I think it’s important to realise its limitations.

Take a look at these two N-body computer simulations of large-scale structure:

The one on the left is a proper simulation of the “cosmic web” which is at least qualitatively realistic, in that in contains filaments, clusters and voids pretty much like what is observed in galaxy surveys.

To make the picture on the right I first  took the Fourier transform of the original  simulation. This approach follows the best advice I ever got from my thesis supervisor: “if you can’t think of anything else to do, try Fourier-transforming everything.”

Anyway each Fourier mode is complex and can therefore be characterized by an amplitude and a phase (the modulus and argument of the complex quantity). What I did next was to randomly reshuffle all the phases while leaving the amplitudes alone. I then performed the inverse Fourier transform to construct the image shown on the right.

What this procedure does is to produce a new image which has exactly the same power spectrum as the first. You might be surprised by how little the pattern on the right resembles that on the left, given that they share this property; the distribution on the right is much fuzzier. In fact, the sharply delineated features  are produced by mode-mode correlations and are therefore not well described by the power spectrum, which involves only the amplitude of each separate mode.

If you’re confused by this, consider the Fourier transforms of (a) white noise and (b) a Dirac delta-function. Both produce flat power-spectra, but they look very different in real space because in (b) all the Fourier modes are correlated in such away that they are in phase at the one location where the pattern is not zero; everywhere else they interfere destructively. In (a) the phases are distributed randomly.

The moral of this is that there is much more to the pattern of galaxy clustering than meets the power spectrum…

8 Comments »