Dark Matter Halos in the context of FRBs [v1.1]



In [1]:

    
# import
from importlib import reload
import numpy as np

from astropy import units
from astropy.coordinates import SkyCoord

from frb import halos

Modified NFW

Init



In [2]:

    
mnfw = halos.ModifiedNFW(alpha=2, y0=2)



In [3]:

    
# Mass in the halo
mnfw.M_b.to('M_sun') * mnfw.f_hot









    Out[3]:





$1.9018718 \times 10^{11} \; \mathrm{M_{\odot}}$



In [4]:

    
mnfw.r200









    Out[4]:





$240.54004 \; \mathrm{kpc}$

Check the Mass



In [5]:

    
Mb_check = mnfw.mass_r(mnfw.r200)
Mb_check









    



/home/xavier/Projects/FRB/frb/halos.py:561: RuntimeWarning: divide by zero encountered in reciprocal
  rho = self.rho0_b * self.f_hot / y**(1-self.alpha) / (self.y0+y)**(2+self.alpha)






    Out[5]:





$1.9030119 \times 10^{11} \; \mathrm{M_{\odot}}$

DM from Rperp



In [6]:

    
Ne = mnfw.Ne_Rperp(100*units.kpc)
Ne









    Out[6]:





$63.075835 \; \mathrm{\frac{pc}{cm^{3}}}$

Milky Way model



In [7]:

    
Galaxy = halos.MilkyWay()



In [8]:

    
Galaxy.M_halo.to('M_sun')









    Out[8]:





$1.5135612 \times 10^{12} \; \mathrm{M_{\odot}}$



In [9]:

    
Galaxy.r200









    Out[9]:





$236.87581 \; \mathrm{kpc}$

Halo DM (from Sun)

This requires NE2001



In [10]:

    
from ne2001.density import NEobject



In [11]:

    
# Zero out inner 10kpc of the halo where we assume the ISM domiantes
Galaxy.zero_inner_ne = 10. # kpc



In [12]:

    
# NE2001 parameters
params = dict(F=1., e_density=1.)



In [13]:

    
model_ne = NEobject(Galaxy.ne, **params)

Pick a sightline



In [14]:

    
l,b = 0., 90.  # deg in Galactic coordinates



In [15]:

    
DM = model_ne.DM(l,b, Galaxy.r200.value)
DM









    Out[15]:





$65.12174 \; \mathrm{\frac{pc}{cm^{3}}}$

M31



In [16]:

    
#reload(halos)
m31 = halos.M31()

DM through M31 from the Sun



In [17]:

    
coord = SkyCoord('J004244.3+413009', unit=(units.hourangle, units.deg))



In [18]:

    
DM = m31.DM_from_Galactic(coord)
DM









    Out[18]:





$133.63955 \; \mathrm{\frac{pc}{cm^{3}}}$

ICM



In [19]:

    
reload(halos)
icm = halos.ICM()



In [20]:

    
icm.M_b.to('M_sun') * icm.f_hot









    Out[20]:





$5.6 \times 10^{13} \; \mathrm{M_{\odot}}$

Test



In [21]:

    
icm.r200









    Out[21]:





$1637.4852 \; \mathrm{kpc}$



In [22]:

    
icm.mass_r(icm.r200)









    Out[22]:





$5.6 \times 10^{13} \; \mathrm{M_{\odot}}$

DM



In [23]:

    
icm.Ne_Rperp(1*units.Mpc)









    Out[23]:





$187.08701 \; \mathrm{\frac{pc}{cm^{3}}}$



In [ ]: