Numerical Calculation of Second Order Perturbations

My new paper, written with Karim Malik, has just been released into the wild. We show that it is possible to numerically simulate second order perturbations for a single scalar field with a canonical action. I’ve been working on this for a long time and learned a lot about the mechanics of inflation in the process. I’ve also churned out quite a lot of python code, and learned some rudimentary parallel programming. We hope to be able to release the code in some form in the future, probably after publication. In the meantime the paper describes the numerical process in some detail.

Now that we have finished the preprint version I can concentrate solely on writing up my PhD thesis. It’s going to be a long few months!