GRB余辉模型
afterglowp的Python项目详细描述
半解析GRB余辉模型
一个python模块来计算grb余辉光曲线和光谱。方法详情见Ryan等人,2019年。建立在范艾尔滕和麦克法登2010年和范艾尔滕2018年的基础上。此代码正在开发中。
归因
如果您在出版物中使用此代码,请按名称引用该软件包,并引用"Ryan,G.,van Eerten,H.,Troja,E.,and Piro,L.(in prep)"。描述代码的出版物将在不久的将来出版。
功能
余辉计算相对论性冲击波正向冲击产生的同步辐射。它包括:
- 通过恒定密度介质的全反相对论激波演化。
- 在冲击面等观测者时间片上的动态积分。
- 喷射扩散的近似处方。
- 任意视角。
- 角结构射流,即e(θ)
- 球速分层流,即e(u)
- 反向喷射。
它对以下方面的支持有限(这些应该被认为是实验性的):
- 初始能量注入
- 逆comption谱
- 空间分辨率强度图
- 早期滑行阶段
它还不包括:
- 外部风介质,即n∮r-2
- 同步辐射自吸收
- 反向冲击发射
Afterglowpy已校准为BoxFit代码(van eerten、van der horst和macfadyen 2011,可在Afterglow库和p在轴上和轴外都能产生类似的顶帽喷气式飞机的光照曲线(使用相同参数时在50%以内)。其喷气模型默认不包括初始滑行阶段,这可能会影响早期观测的预测。
安装/建筑
余辉可通过pip
:
$ pip install afterglowpy
如果您正在处理此repo的本地副本并希望从源安装,则可以从项目的顶级目录运行以下命令。
$ pip install -e .
使用
此界面将在不久的将来更新为更合理的界面
在python代码中,使用 感兴趣的主要功能是 对于喷射状余辉( 对于茧状余辉( 关键字参数是:import afterglowpy as grb>导入库
电子中的热能分数grb.fluxdensity(t,nu,jettype,spectype,*pars,**kwargs)
。有关一个简单的示例,请参见tests/plotlc.py
。jettype
可以是-1(顶帽)、0(高斯)、1(幂律w/核)、2(高斯w/核)、3(cocoon)或4(平滑幂律)。spectype
可以是0(全局冷却时间,无逆康普顿)或1(全局冷却时间,逆康普顿)。jettypes
-2、-1、0、1、2和4)pars
有14个位置参数:0 thetav
以弧度表示的视角
1e0
轴上各向同性等效能量,单位为erg2 thetac
以弧度表示的喷射核心的半宽(特定于喷射类型)3喷嘴的"机翼"截断角,单位为弧度
4b
幂律结构的幂,θ-b5l0
能量注入基准亮度,单位为erg/s,通常为0。6q
能量注入的时间幂律指数,通常为0。7 ts
能量注入基准时间标度,单位为秒,通常为0。8n0
ISM的数量密度,单位为厘米-39p
电子分布幂律指数(p>;2)12xiu n
加速电子的分数13 d_l
光度距离,单位为厘米jettype
3)pars
有14个位置参数:0 umax
初始最大流出速度4-velocity1umin
最小流出4-速度2e i
速度分布中的基准能量,e(>;u)=eiu-k3k
能量速度分布幂律指数4mej
太阳质量中"umax"处的物质质量5l0
能量注入基准亮度,单位为erg/s,通常为0。6q
能量注入的时间幂律指数,通常为0。7 ts
能量注入基准时间标度,单位为秒,通常为0。8n0
ISM的数量密度,单位为厘米-39p
电子分布幂律指数(p>;2)12xiu n
加速电子的分数13 d_l
光度距离,单位为厘米z
redshift(默认为0)tres
冲击演化方案的时间分辨率,每十年的时间采样点数latres
结构化喷流的纬度分辨率,每个θc的壳数
rtol
通量积分的目标相对公差spread
boolean(默认为true),是否允许喷射扩散。推荐PyPI第三方库