You can get the data of the effective degrees of freedom here.

The data format is:

Temperature [GeV], g_{s}, g_{\rho}

This data file is directly available in the Program Micromegas v_5.

In order to use the present estimation, simply put

loadHeffGeff("path/to/eos2020.dat");

in the main program.

Though we esimate g_{s}(T_0) = 3.931, the Micromegas seems to adopt g_{s}(0) = 3.910 when converting Y_{DM} to Omega h^2.

Therefore multiplying factor 1.005 by the Micromegas output with EOS2020.dat reproduces well our estimation.

There are several uncertainties of the estimation of the degrees of freedom.

We also provide the data for the upper error and lower error.

See also EOS2018 for the application to the gravitational wave spectrum.

We also apply it to simplified WIMP dark matter models.

We estimate the DM annihilation cross section to realize the current observation \[\Omega h^2 = 0.120 \pm 0.001.\]

S-wave

real scalar, complex scalar or Majorana fermion, vector, Dirac fermion

P-wave

real scalar, complex scalar or Majorana fermion, vector, Dirac fermion

The data format is:

DM Mass [GeV], sigma v (a_s or b_p) [cm^3/s /10^{-26}], error from DM abudance mesurement, error from EOS

If you find any problems or request additional data, please contact me (satoshi.shirai_at_ipmu.jp)