Trispectrum during Inflation

After a lot of work, “Large trispectrum in two-field slow-roll inflation” was released on the arXiv yesterday as arXiv:1203.6844. In this article Joe Elliston, Laila Alabidi, David Mulryne, Reza Tavakol and I look at the generation of higher order statistics during inflation in the early universe.

In the early universe the curvature perturbations, which later are seen as temperature fluctuations in the Cosmic Microwave Background (CMB), are initially thought to be Gaussian, but can become skewed during inflation depending on the physics of their evolution. In the last few years a lot of work has been done to both find evidence of this non-Gaussianity, and to construct physical models in which it is generated in the early universe.

In the past most of the focus has been on the 3-point function or bispectrum, and discussion of non-gaussianity has boiled down to finding bounds on the parameter f_{\mathrm{NL}}. In terms of the CMB the bispectrum in essence considers whether the temperature of three points on the sky is correlated. The WMAP satellite has not seen any definitive evidence of a non-zero value for f_{\mathrm{NL}} but the Planck satellite should be able to detect a signal if it is moderately large. In this work we look beyond the bispectrum to the 4-point function or trispectrum.  For the trispectrum the correlation we attempt to measure is between four different points on the sky.

In this work we have tried to find models which generate a large value for the trispectrum during inflation. We have found some new expressions for the parameters f_{\mathrm{NL}}\tau_{\mathrm{NL}} and g_{\mathrm{NL}}. The last two of these parametrise two different contributions to the trispectrum. The bottom line is that it is quite difficult to find conditions where the trispectrum can be large, at least under the assumptions we made of sum- and product-separable potentials using the \delta N formalism.

In the course of searching for models which give large values to these parameters we plotted the coefficient functions which need to be large as heatmaps, following Byrnes et al (arXiv:0807.1101). In order to generate these heatmaps I relied on the combination of Python, Numpy and Matplotlib, which I have used before on Pyflation. The script I used to generate the heatmap figures in the paper is now available as a repository on Bitbucket.

PS Someone really needs to work on the Wikipedia Non-Gaussianity page!


Two new papers on the arXiv

Last Tuesday was a bit hectic for me as I tried to coordinate the last minute changes needed to put two papers on the arXiv servers for the next day. The two articles which are now available are numbered 1111.6919 and 1111.6940:

  • Calculating Non-adiabatic Pressure Perturbations during Multi-field Inflation
    Ian Huston, Adam J. Christopherson
    Abstract: Isocurvature perturbations naturally occur in models of inflation consisting of more than one scalar field. In this paper we calculate the spectrum of isocurvature perturbations generated at the end of inflation for three different inflationary models consisting of two canonical scalar fields. The amount of non-adiabatic pressure present at the end of inflation can have observational consequences through the generation of vorticity and subsequently the sourcing of B-mode polarisation. We compare two different definitions of isocurvature perturbations and show how these quantities evolve in different ways during inflation. Our results are calculated using the open source Pyflation numerical package which is available to download.
  • Quantifying the behaviour of curvature perturbations during inflation
    Ellie Nalson, Adam J. Christopherson, Ian Huston, Karim A. Malik
    Abstract: How much does the curvature perturbation change after it leaves the horizon, and when should one evaluate the power spectrum? To answer these questions we study single field inflation models numerically, and compare the evolution of different curvature perturbations from horizon crossing to the end of inflation. In particular we calculate the number of efolds it takes for the curvature perturbation at a given wavenumber to settle down to within a given fraction of their value at the end of inflation. We find that e.g. in chaotic inflation, the amplitude of the comoving and the curvature perturbation on uniform density hypersurfaces differ by up to 180 % at horizon crossing assuming the same amplitude at the end of inflation, and that it takes approximately 3 efolds for the curvature perturbation to be within 1 % of its value at the end of inflation.

As part of the work on the first of these papers, my Pyflation code has been updated to version 0.2 and now allows multiple fields systems to be used at first order.

There seemed to be a glut of papers released last Wednesday, with the post-Thanksgiving queue being worked through and people aiming to apply for jobs in December. Apologies to those irked by the overwhelming astro-ph listing that day for my part in it!


New paper and Pyflation software package

PyflationMy latest paper has just hit the arXiv and is now available. The paper builds on the numerical work I previously completed on cosmological perturbations beyond linear order. The new results do not assume slow-roll in the calculation of the source term for the second order equations of motion and so allow a much greater range of potentials to be analysed. The paper is called “Second Order Perturbations During Inflation Beyond Slow-roll” and already has a record on SPIRES and Inspire Beta.

Accompanying the paper is the release of the software package used to create the results. Pyflation is a Python package which calculates the first and second order perturbation results including the source term required at each time step. It requires Numpy, Scipy and PyTables and a fairly recent version of Python. The full calculation for a large range of Fourier wavemodes takes a long time on a single CPU, but there is support for splitting the calculation into many separate jobs which can be queued on a multi-core or cluster based setup. Alternatively there is also support for calculating the second order results for a single wavemode, which significantly reduces the execution time.

The Pyflation website contains all the details and download links, installation notes and links to the relevant papers. We’ve released the package under a very liberal open source license (modified BSD license) but as citations are increasingly the oil that greases the wheels of academia we ask that anyone who uses results they have obtained with Pyflation to please cite one or more of the related papers.


Checklist for arXiv submission

I am currently finishing up a paper that is about to be submitted to the arXiv and I thought I would go through the list of things I normally do just before sending the work off. It’s mainly a common sense list but in the rush to get something out it is good to have a list to work from to make sure you don’t miss anything (misspelling collaborators names is not helpful!).

  • Check bib style is correct.
    If the journal you are planning to submit to has a particular house style for the bibliography it is probably worth using it in the arXiv submission.
  • Check bibliography text is correct.
    Even though I think BiBTeX is the great, and much easier than preparing the bibliography by hand (which some people still do), there can be instances when a stray space or mistakenly capitalised letter appear. If you get your BiBTeX entries from SPIRES or Inspire this is mostly taken care of, but there is always the odd paper or book entry that you have typed in by hand. Check the final output in the paper, not just the .bib file contents. Continue reading…

Inspire Beta: New Interface to SPIRES database

INSPIRE BetaA colleague mentioned today that the front page of the venerable SPIRES database of High Energy Physics papers is now promoting the new Inspire interface which was announced a few years ago.

The website for the current beta phase of the project is

It is unclear to me whether “Beta” is part of the site name, as suggested by the URL and the text on the INSPIRE page, or whether this is just the beta phase of the INSPIRE project as the SPIRES homepage seems to imply. It would certainly be an odd decision to use a different URL for the beta phase and force everyone to change bookmarks, references in blogs, literature etc., once the beta phase is over.

These are quick first impressions because I haven’t had much time to use the new service. First off it is fast. Very fast. So fast that when searching for my name INSPIRE claims that the “Search took 0.00 seconds”. It feels almost instantaneous. This might be because of a light load before the hordes using SPIRES are switched over. It is certainly an improvement on the interminable and often futile stretches of time needed with the SPIRES engine.

Continue reading…