Sobol, Oleksandr O.Oleksandr O.Sobolvon Eckardstein, RichardRichardvon EckardsteinSchmitz, KaiKaiSchmitzKoch, Elias R.Elias R.KochThiele, UweUweThieleGurevich, Svetlana V.Svetlana V.Gurevich2026-06-242026-06-242026-05-25Sobol, O. O., von Eckardstein, R., Schmitz, K., Koch, E. R., Thiele, U., Gurevich, S. V. (2026). Numerical data for the paper "Nonlinear physics of axion inflation" [Data set]. https://doi.org/10.5281/zenodo.2037273510.5281/zenodo.20372735https://doi.org/10.5281/zenodo.20372735https://ir.library.knu.ua/handle/15071834/24958Data 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]ukdataset