We derive model-independent astrophysical constraints on leptophilic dark matter (DM), considering its thermal production in a supernova core and taking into account core temperature fluctuations within the framework of $q$-deformed Tsallis statistics. In an effective field theory approach, where the DM fermions interact with the Standard Model via dimension-six operators of either scalar-pseudoscalar, vector-axial vector, or tensor-axial tensor type, we obtain bounds on the effective cut-off scale $\Lambda$ from supernova cooling and free-streaming of DM from supernova core, and from thermal relic density considerations, depending on the DM mass and the $q$-deformation parameter. Using Raffelt's criterion on the energy loss rate from SN1987A, we obtain a lower bound on $\Lda \gtrsim 3$ (12) TeV corresponding to $q = 1.0~(1.1)$ and an average supernova core temperature of $T_{\rm SN}=30$ MeV. From the optical depth criterion on the free-streaming of DM fermions from the outer 10\% of the SN1987A core, the cooling bound is restricted to $\Lda \gtrsim 1$ TeV. Both cooling and free-streaming bounds are insensitive to the DM mass $m_\chi$ for $m_\chi\lesssim T_{\rm SN}$, whereas for $m_\chi\gg T_{\rm SN}$, the bounds weaken significantly due to the Boltzmann-suppression of the DM number density. We also calculate the thermal relic density of the DM particles in this setup and find that it imposes an upper bound on $\Lambda^4/m_\chi^2$, which together with the cooling/free-streaming bound significantly constrains light leptophilic DM.

The cold dark matter fermions with mass MeV scale, pair produced inside the supernova SN1987A core, can freely stream away from the supernovae and hence contributes to its energy loss rate. Similar type of DM fermions(having similar kind of coupling to the standard model photon), produced from some other sources earlier, could have contributed to the relic density of the Universe. Working in a theory with an effective dark matter-photon coupling (inversely proportional to the scale $\Lambda$) in the formalism of Tsallis statistics, we find the dark matter contribution to the relic density and obtain a upper bound on $\Lambda$ using the experimental bound on the relic density for cold non-baryonic matter i.e. $\Omega h^2 = 0.1186 \pm 0.0020 $. The upper bound obtained from the relic density is shown with the lower bound obtained from the Raffelt's criterion on the emissibity rate of the supernovae SN1987A energy loss $\dot{\varepsilon}(e^+ e^- \to \chi \overline{\chi}) \le 10^{19}~\rm{erg~g^{-1}s^{-1}}$ and the optical depth criteria on the free streaming of the dark matter fermion (produced inside the supernovae core). As the deformation parameter $q$ changes from $1.0$ (undeformed scenario) to $1.1$(deformed scenario), the relic density bound on $\Lambda$ is found to vary from $ \sim 4.9 \times 10^7 $ TeV to $1.6 \times 10^8$ TeV for a fermion dark matter($\chi$) of mass $m_\chi = 30~\rm{MeV}$, which is almost $10$ times more than the lower bound obtained from the SN1987A energy loss rate and the optical depth criteria. \noindent {{\bf Keywords}: Dark matter, Relic density, Supernova cooling, Tsallis statistics, free-streaming, } }

Light dark matter($\simeq 1-30~\rm{MeV}$) particles pair produced in electron-positron annihilation $ e^-e^+ \stackrel{\gamma}{\longrightarrow} \chi \bar{\chi}$ inside the supernova core can take away the energy released in the supernova SN1987A explosion. Working within the formalism of $q$-deformed statistics (with the average value of the supernovae core temperature(fluctuating) being $T_{SN} = 30~\rm{MeV}$) and using the Raffelt's criterion on the emissivity for any new channel $\dot{\varepsilon}(e^+ e^- \to \chi \overline{\chi}) \le 10^{19}~{erg~g^{-1}s^{-1}}$, we find that as the deformation parameter $q$ changes from $1.0$ (undeformed scenario) to $1.1$(deformed scenario), the lower bound on the scale $\Lambda$ of the dark matter effective theory varies from $3.3\times 10^6$ TeV to $3.2 \times 10^7$ TeV for a dark matter fermion of mass $m_\chi = 30~\rm{MeV}$. Using the optical depth criteria on the free streaming of the dark matter fermion, we find the lower bound on $\Lambda \sim 10^{8}~\rm{TeV}$ for $m_\chi = 30~\rm{MeV}$. In a scenerio,where the dark matter fermions are pair produced in the outermost sector of the supernova core (with radius $0.9 R_c \le r \le R_c$, $R_c (=10~\rm{km})$ being the supernova core radius or the radius of proto-neutron star), we find that the bound on $\Lambda$ ($\sim 3 \times 10^7$ TeV)obtained from SN cooling criteria (Raffelt's criteria) is comparable with the bound obtained from free streaming (optical depth criterion) for light fermion dark matter of mass $m_{\chi}=10 - 30$ MeV.