Programs and data for Figure 9 for "Instability of global de Sitter space to particle creation," Physical Review D89, 104039 (2014)

Abstract

Each fortran program is for different values of m and k and solves the mode equation for the Bunch-Davies state for that value of m and k. For this figure m = 10. The starting values for the modes which come from the Mathematica notebook were copied and pasted into the Fortran programs. Note that the Fortran programs all use the Numerical Recipies routine odebs and its subroutines to numerically solve the mode equation. This can easily be replaced by a call to any differential equation solver.