For our numerical study we implement the model in LanHEP 3.2 [1]} and perform the calculation of \(\Omega _{\rm DM} h^2\)  using MicrOMEGAs 5.0.6 [2]}. Moreover, we perform an independent calculation in Mathematica.
