# ON THE HEAVY RELIC NEUTRINO - GALACTIC GAMMA HALO CONNECTION

## 1 Introduction

The recent observations of EGRET telescope show a diffuse $`\gamma `$ \- ray emission in the halo of our galaxy . Actually models trying to explain a gamma halo range between   
a) a high galactic latitude distribution of high energy cosmic ray sources (fast running pulsar ”Geminga” like) ,   
b) collisions of high energy protons (tens of GeV) in molecular clouds (mostly $`H_2`$), whose formation is favored in galactic halos ,   
c) a cold dark matter scenario with neutralino, the lightest supersymmetric particle in the minimal SUSY model, whose annihilations in heavy fermions ($`c\overline{c},b\overline{b},t\overline{t}`$) and bosons ($`W^+W^{},ZZ,gg`$) could lead to gamma secondaries emission ,   
d) Inverse Compton Scattering (ICS) of interstellar photons off cosmic ray electrons with a spectra harder than previous predictions .   
A heavy Dirac neutrino of a fourth generation in a Cold Dark Matter (CDM) model, at masses $`m_N>M_Z/2`$, may offer an elegant solution to the gamma ray signal detected, though it can not solve all the dark matter problem in the halo.   
Neutrino annihilations at high galactic latitude could produce   
1) $`\gamma `$ rays by ICS of relativistic electrons (primaries or as secondary decay products of heavier leptons in the annihilation chains $`N\overline{N}l^+l^{}e^+e^{}`$) onto thermal photons ($`IR`$, $`optical`$) near and above the galactic plane,   
2) direct gammas by neutral secondary pions decay.   
The estimated flux we derived here is roughly close to EGRET results in the hypothesis of a smooth and homogeneous galactic halo.   

## 2 The heavy neutrino model

LEP I has fixed sever constrains on the number of ”light” ($`m_\nu M_Z/2`$) neutrino families from Z width data. However there is no experimental prohibition on the existence of an additional heavy neutral lepton with mass $`M_N>M_Z/2`$. Such a heavy stable neutrino belonging to a fourth fermion family was introduced nearly 20 years ago as a CDM candidate ,. An apparently simple fourth generation model is not so easy to build. Anyway models predicting a fourth heavy stable neutrino can be found in ,,. In an expanding universe scenario heavy neutrinos decouple from a thermal equilibrium condition when global temperature drops below their rest mass energy ($`T<m_N`$) and weak interactions become too slow to keep neutrinos in equilibrium with the cosmological fluid. From this moment the only change in neutrino density is due to cosmic expansion. The relic abundance is given by 

$$n_N\frac{2\times 10^{18}}{g_{}^{1/2}M_pm_N(\overline{\sigma \beta })_f}\left[40+\mathrm{ln}\left(\frac{g_s}{g_{}^{1/2}}M_pm_N(\overline{\sigma \beta })_f\right)\right]n_\gamma (T)$$
(1)   
where $`g^{}=N_{bos}+\frac{7}{8}N_{ferm}`$ is the number of effective degrees of freedom at temperature T, $`g_s`$ is the number of particle spin states $`M_p`$ is the proton mass, $`(\overline{\sigma \beta })_f`$ is the thermally averaged annihilation cross section at freeze out, $`n_\gamma =0.24T^3`$ is the cosmic photon number density.   
Annihilations of heavy neutrinos in the universe happen through two main channels   
channel 1) $`N\overline{N}f\overline{f}`$ if $`M_Z/2<m_N<m_W`$   
where the cross section decreases as $`4m_N^2/(4m_N^2M_Z^2)^2`$ for growing $`m_N`$   
channel 2) $`N\overline{N}W^+W^{}`$ if $`m_N>m_W`$   
with a cross section growing like $`m_N^2`$ for increasing $`m_N`$ .   
As $`\rho _N\sigma ^1`$, neutrino relic density exhibits a maximum (with $`\rho _{max}/\rho _c10^2h^2`$) near $`m_NM_W`$ and then starts to decrease as $`m_N^2`$ in the mass range $`m_N>m_W`$, without reaching the critical value $`\mathrm{\Omega }=1`$, at least in the mass range where Standard Model may be applied ($`m_N<1TeV`$).   
Clustering during galactic structure formation determine an increase in neutrino density ,, that in the central part of the galaxy could be as large as 5 $`÷`$ 7 orders of magnitude (the exact value depends on the other CDM, HDM densities and masses) .   
A spherical halo around our galaxy made of heavy neutrinos is the model proposed in order to explain $`\gamma `$ emission observed at high galactic latitude. In a galactic halo, neutrinos with a higher density distribution could annihilate again leading to a flux of ordinary particles beyond the galactic plane, potential sources of high energy radiation.   
Constrains on neutrino mass come either from cosmological data (not too high pollution of $`e^+e^{}`$ cosmic rays is observed in the range $`M_Z<m_N<300GeV`$ ) or DAMA detector, where recent signals could be attributed to heavy neutrino in the mass window $`45GeV<m_N<50GeV`$.

## 3 Neutrino annihilation products as gamma ray source

### 3.1 Relativistic electron pairs: ICS on the galactic interstellar radiation field(ISRF)

Heavy neutrinos could directly annihilate in relativistic electron pairs (either prompt ones or born through secondary decay processes of heavier particles $`\mu ,\tau `$, as in channel 1 way) . Channel 2 leads to electron pairs through leptonic decay of W ($`N\overline{N}W^+W^{}l^+l^{}`$). Electron pairs may be generated even by W, Z hadronic decay through charged pions and neutrons production.   
Electrons and positrons are trapped by galactic magnetic field, and propagating through the Galaxy loose either ”memory” of their ”place of birth” as well as energy for bremsstrahlung, synchrotron or ICS. These processes determine a broadening of different electron ”lines” ($`N\overline{N}l^+l^{}e^+e^{}`$), so that $`e^{}(e^+)`$ spectra (even considering electrons and positrons that come from hadron decays) are at final stages described by the consequent approximated power law $`J=KE^\alpha `$ (where K is a normalization constant). Numerical simulation of $`N\overline{N}`$ annihilation performed with the package PYTHIA 5.7  with suitable modifications to include a fourth generation of fermions, show that such $`N\overline{N}`$ relic electron fluxes are considerably lower than observed neighbor galactic background (in the range of masses $`45GeV<m_N<M_Z`$). Such a cosmic ray input can not be used to confirm or refute heavy neutrino presence in galactic halo.   
ICS of ”soft” background photons (isotropic CBR or anisotropic infrared and optical interstellar radiation field) off relativistic electrons created out of the galactic plane is a possible source of radiation. Energies order of magnitude is fixed by the ICS characteristic relation $`E_\gamma =4/3ϵ_{ph}(E_e/m_ec^2)^2`$, where $`ϵ_{ph}`$ is the target photon energy.   
We excluded here collisions with microwave photons which would require too large neutrino masses ($`m_N>1TeV`$). No clear theory is still available for such a heavy particle. ICS on IR and optical photons needs respectively $`E_e50GeV`$ and $`E_e10GeV`$, and is more efficient in gamma ray production. The Galaxy is a disk-like radiative source of radius $``$ 15 $`kpc`$, so the interstellar radiation field has a vertical extent of several kpc. We assumed this radiation to be represented by the obvious scaling law

$$n_{ph}(r)=\frac{n_{ph}(0)}{1+r^2/a_\gamma ^2}.$$
(2)   
where $`r`$ is the distance from the galactic plane, and $`a_\gamma =10kpc`$ is the characteristic length of interstellar radiation distribution in the Galaxy. Photon density could be considered roughly constant in a region of radius $`a_\gamma `$ .   
An electron distribution ($`KE^\alpha `$) interacting by ICS with photons at energy $`ϵ_{ph}`$ and density $`n_{ph}`$ generate radiation whose intensity is 

$$J_\gamma (E_\gamma )=\frac{2}{3}Ka_\gamma n_{ph}\sigma _T\left(\frac{\overline{ϵ}_{ph}}{(mc^2)^2}\right)^{(\alpha 1)/2}E_\gamma ^{(\alpha +1)/2}$$
(3)   
where $`n_{ph}`$ is the target photon background density, $`\overline{ϵ}_{ph}`$ its average energy and $`\sigma _T`$ is the Thomson cross section.   
Gamma intensity has been calculated for $`m_N=45,\mathrm{\hspace{0.17em}50},\mathrm{\hspace{0.17em}100},\mathrm{\hspace{0.17em}300}GeV`$. The largest flux has been obtained for ICS on optical photons.   
For $`m_N=50GeV`$ the calculated flux is

$$\frac{dN_\gamma }{dSdtd\mathrm{\Omega }dE_\gamma }210^7A(\psi )\left(\frac{E_\gamma }{GeV}\right)^{1.55}\left(\frac{a_\gamma }{10kpc}\right)cm^2s^1sr^1GeV^1$$
(4)   
and for $`m_N=100GeV`$ one finds

$$\frac{dN_\gamma }{dSdtd\mathrm{\Omega }dE_\gamma }310^7A(\psi )\left(\frac{E_\gamma }{GeV}\right)^{1.5}\left(\frac{a_\gamma }{10kpc}\right)cm^2s^1sr^1GeV^1.$$
(5)   
$`A(\psi )`$ is the adimensional integral of interstellar photon density along the line of sight L, defined by the angular coordinate $`\psi `$ (angle between L and the direction of the galactic centre). $`A(\psi )`$ is of few unities and corresponds to

$$A(\psi )=\frac{1}{a_\gamma }_{lineofsight}\frac{dr(\psi )}{(1+r(\psi )^2/a_\gamma ^2)}$$
(6)   
Gamma intensity due to infrared background is less abundant than optical photons as a consequence of the spectral power law $`E^\alpha `$.   
Assuming an average $`N\overline{N}`$ clustering $`\rho _N^{gal}/\rho _N^{cosm}=10^6`$, the flux obtained for two values of $`m_N`$   
1)$`\mathrm{\Phi }_\gamma (E>1GeV)410^7A(\psi )\left(\frac{a_\gamma }{10kpc}\right)cm^2s^1sr^1`$, ($`m_N50GeV`$),   
2)$`\mathrm{\Phi }_\gamma (E>1GeV)610^7A(\psi )\left(\frac{a_\gamma }{10kpc}\right)cm^2s^1sr^1`$ ($`m_N100GeV`$),   
is comparable with EGRET observations:   
$`\mathrm{\Phi }_\gamma (E>1GeV)810^7cm^2s^1sr^1`$.   
Additional tests of this model could be obtained with the nearly detectable signal of hundreds KeV radiation in the halo with flux $`J_\gamma 10^2cm^2s^1sr^1`$ at peak energy $`E_\gamma 300keV`$ as well as a flux $`J_X0.3cm^2s^1sr^1`$ at $`E_X3KeV`$ due to ICS of tens of $`GeV`$ and $`GeV`$ electrons with CBR. An addtional parassite radio background arises at high galactic latitudes due to synchrotron losses of the same electrons at $`E_e10GeV`$, with typical density flux

$$J_{sync}510^4\left(\frac{B}{1\mu G}\right)\left(\frac{\gamma }{210^4}\right)^2\left(\frac{U_{rad}}{0.2eVcm^3}\right)^1Jy.$$
(7)   
and average frequency

$$\nu =\gamma ^2\left(\frac{eB}{2\pi m_e}\right)1GHz\left(\frac{B}{1\mu G}\right)\left(\frac{\gamma }{210^4}\right)^2$$
(8)   
where we used as characteristic scale for the magnetic field B = 1 $`\mu G`$,and as optical background energy density $`U_{rad}=0.2eVcm^3`$.   

### 3.2 Annihilations in gamma photons

Gamma radiation could also be produced in neutrino annihilations due to neutral pions secondaries in Z, W hadronic decay. This kind of emission does not need to introduce any kind of radiative background and is directly related to neutrino distribution in the halo.   
Photon flux is described by the following expression:

$$J_\gamma =\frac{1}{4\pi m_N^2}\underset{i}{}\sigma _iv\frac{dN^i}{dE}_{lineofsight}\rho ^2(r)𝑑r(\psi )$$
(9)   
where $`\psi `$ is the angle between the line of sight and the galactic center, $`\rho (r)`$ is heavy neutrino density as a function of galactocentric radius, and $`_i\sigma _iv\frac{dN^i}{dE}`$ counts all possible final photon channels ($`\frac{dN^i}{dE}`$) which could contribute to gamma photons emission. The integral of neutrino density along the line of sight L depends on the halo model chosen for dark matter distribution, which is generally described as

$$\rho (r)\frac{1}{(\frac{r}{a})^\gamma [1+(\frac{r}{a})^\alpha ]^{(\beta \gamma )/\alpha }}$$
(10)   
The simplest density profile is described by an isothermal sphere with $`\alpha =2,\beta =0,\gamma =0`$) and length scale $`a10kpc`$.   
In the spherical model the square density integral leads to an adimensional intensity $`I(\psi )`$

$$I(\psi )=\frac{1}{(1+(r/a)^2)^2}𝑑r(\psi )/a;$$

this intensity has a characteristic behaviour which is maximum in the direction of the galactic center ($`\psi =0`$), and then decreases for $`0<\psi <\pi `$, but it doesn’t vary more than a factor ten with the angular coordinate.   
At high latitudes $`I(\psi )`$ is generally of order unity.   
Models with a singular behaviour towards the galactic center or which postulate a clumpy distribution of dark matter could contribute to enhance the total gamma flux in the halo, but we shall neglect them here.   
Monte Carlo simulations of neutrino annihilations  have been also used to compare EGRET flux, showing that it is possible to extrapolate a power law for gamma spectrum.   
An approximated integral flux for $`m_N=50GeV`$ and $`a10kpc`$ is roughly

$$\mathrm{\Phi }_\gamma >610^7I(\psi )cm^2s^1sr^1,$$
(11)   
while for $`m_N=100GeV`$ at $`a10kpc`$

$$\mathrm{\Phi }_\gamma >410^7I(\psi )cm^2s^1sr^1,$$
(12)   
In conclusion there are at least two independent processes able to solve the puzzle of a GeV gamma halo by the role of a cosmic relic DM made of fourth generation heavy neutrinos. The neutrino masse are constrained into a narrow energy window ($`45GeV<m_N<60GeV`$), in order to combine at once the DAMA data and the other underground detector. This reality will be soon confirmed or excluded by LEP II search for $`e^{}e^+N\overline{N}\gamma `$ events  in this energy band.