Data files containing numerical results illustrating the fundings of the paper "Nonlinear physics of axion inflation" https://arxiv.org/abs/2603.02570
%---------------------------------------------------------------------------------------------------------------------------------------%Relation to figures in the paper.%---------------------------------------------------------------------------------------------------------------------------------------%
Figures 1 and 7:- background color --> file "Re-zeta.dat"- red solid line --> file "SR-UB-boundary.dat"- red dashed line --> file "SB-UB-boundary.dat"
Figures 2 and 3:- red solid line --> file "SR-UB-boundary.dat"- red dashed line --> file "SB-UB-boundary.dat"- blue solid line --> file "delta-KG-criterion.dat"
Figure 4:- panel (a): green dashed line, red dots -->  file "bifurcation.dat"- panel (a): green solid and dashed lines --> solution \xi=\xi(b) to Eq. (3.7) for fixed v=200; function \bar{g}_0(\xi) tabulated in file "g0.dat"- panle (b) --> file "phase-space-xi2.89-b0.10853.dat"- panel (c) --> file "phase-space-xi2.977-b0.16972.dat"
Figure 5:- panels (a1) and (b1) --> file "phase-space-v3.6e7-b1.4e-6.dat"- panel (c1) --> file "Fourier-v3.6e7-b1.4e-6.dat" - panels (a2) and (b2) --> file "phase-space-v3.6e7-b1.9e-6.dat"- panel (c2) --> file "Fourier-v3.6e7-b1.9e-6.dat" - panels (a3) and (b3) --> file "phase-space-v3.6e7-b3.3e-6.dat"- panel (c3) --> file "Fourier-v3.6e7-b3.3e-6.dat" - panels (a4) and (b4) --> file "phase-space-v3.6e7-b2.4e-5.dat"- panel (c4) --> file "Fourier-v3.6e7-b2.4e-5.dat" 
Figure 6 (a):- background color --> solution of Eq. (4.3) w.r.t. \beta, where v^{new}_{thr} is tabulated in file "SR-UB-boundary.dat"- white contour lines --> Eq. (4.3) for fixed values of \beta, where v^{new}_{thr} is tabulated in file "SR-UB-boundary.dat"- yellow [cyan] solid and dashed lines --> parametric curves given by Eq. (4.4) with potential in Eq. (4.1) [Eq. (4.2)]
Figure 6 (b):- yellow solid curve --> file "m2-b20_gef.dat"- yellow dashed curve --> file "m2-b20_inf.dat"- cyan solid curve --> file "tanh-b60_gef.dat"- cyan dashed curve --> file "tanh-b60_inf.dat"
Figure 6(c):- yellow solid curve --> files "m2-b20_gef.dat" and "m2-b20_inf.dat"- cyan solid curve --> files "tanh-b60_gef.dat" and "tanh-b60_inf.dat"
%---------------------------------------------------------------------------------------------------------------------------------------%Contents of each file.%---------------------------------------------------------------------------------------------------------------------------------------%
File "g0.dat":- 1 column --> values of \xi [dimensionless]- 2 column --> values of log_10(g_0(\xi)) [dimensionless]
Files "SR-UB-boundary.dat",  "SB-UB-boundary.dat", and "delta-KG-criterion.dat":- 1 column --> values of log_10(v) [dimensionless]- 2 column --> values of log_10(b) [dimensionless]- 3 column --> values of \xi [dimensionless]
File "Re-zeta.dat":- 1 column --> values of log_10(v) [dimensionless]- 2 column --> values of log_10(b) [dimensionless]- 3 column --> values of Re(\zeta_1) [dimensionless]
File "bifurcation.dat":- 1 column --> values of b [dimensionless]- 2 column --> values of \overline(\xi) [dimensionless]- 3 column --> values of \xi_{min} [dimensionless]- 4 column --> values of \xi_{max} [dimensionless]
Files "phase-space-xi2.89-b0.10853.dat" and "phase-space-xi2.977-b0.16972.dat":- 1 column --> values of \xi [dimensionless]- 2 column --> values of d\xi/dN [dimensionless]- 3 column --> values of \rho_E/\rho_{E,0}-1 [dimensionless]
Files "phase-space-v3.6e7-b1.4e-6.dat", "phase-space-v3.6e7-b1.9e-6.dat", "phase-space-v3.6e7-b3.3e-6.dat", and "phase-space-v3.6e7-b2.4e-5.dat":- 1 column --> values of \Delta N (number of e-folds) [dimensionless]- 2 column --> values of \xi [dimensionless]- 3 column --> values of d\xi/dN [dimensionless]- 4 column --> values of \rho_E/\rho_{E,0} [dimensionless]
Files "Fourier-v3.6e7-b1.4e-6.dat", "Fourier-v3.6e7-b1.9e-6.dat", "Fourier-v3.6e7-b3.3e-6.dat", and "Fourier-v3.6e7-b2.4e-5.dat":- 1 column --> values of Fourier frequencies \omega_n (w.r.t. number of e-folds as a time variable) [dimensionless]- 2 column --> values of Fourier coefficients A_n [dimensionless]- 3 column --> values of Fourier coefficients B_n [dimensionless]- 4 column --> values of Fourier coefficients D_n [dimensionless]
Files "m2-b20_gef.dat" and "tanh-b60_gef.dat":- 1 column --> time t [reduced Planck units, M_P^-1]- 2 column --> number of e-folds N [dimensionless]- 3 column --> inflaton field \phi [reduced Planck units, M_P]- 4 column --> inflaton velocity \dot\phi [reduced Planck units, M_P^2]- 5 column --> inflaton acceleration \ddot\phi [reduced Planck units, M_P^3]- 6 column --> Hubble rate H [reduced Planck units, M_P]- 7 column --> instability variable \xi [dimensionless]- 8 column --> mean value of elecric field squared [reducel Planck units, M_P^4]- 9 column --> mean value of magnetic field squared [reducel Planck units, M_P^4]- 10 column --> mean value of Chern-Pontryagin density [reducel Planck units, M_P^4]- 11 column --> instability threshold momentum k_h [reduced Planck units, M_P]
 Files "m2-b20_inf.dat" and "tanh-b60_inf.dat":- 1 column --> time t [reduced Planck units, M_P^-1]- 2 column --> number of e-folds N [dimensionless]- 3 column --> inflaton field \phi [reduced Planck units, M_P]- 4 column --> inflaton velocity \dot\phi [reduced Planck units, M_P^2]- 5 column --> inflaton acceleration \ddot\phi [reduced Planck units, M_P^3]- 6 column --> Hubble rate H [reduced Planck units, M_P]