next up previous
: 光学的に薄い熱制動放射 (試験公開) : スペクトル解析 : スペクトル解析

光学的に薄い非熱ジャイロシンクロトロン放射 (試験公開)

注意:改良についてのコメント歓迎。

放射領域物理量から放射を求めることができる。
Dulk (1985) の方法とRamaty (1969)の方法とが用意されている。 ここでは電子の分布関数が以下のようなものであると 仮定している。

\begin{displaymath}
{{d NV(E)} \over {d E}} = K \left({{E}\over{E_{0}}}\right)^{-\delta}
\end{displaymath}

$E$は電子のエネルギー(keV)・$NV(E)$はエネルギー$E$以上の電子数(particles)。 この分布はパラメータ$\delta$$E_{0}$$K$とによって記述されるが、$K$のかわりに $NV(E_0)=K/(\delta -1)/E_0^{\delta-1}$を与える方が物理的な意味が明確になるので そちらを実際は用いる。以下のプロシジャではDulk (1985)にならい、 $E_0 = 10\ {\rm keV}$で固定している。

(a) 物理量から放射を求める

放射領域の、電子分布べき指数$\delta$・磁場強度(G) ・視線と磁場とのなす角度(度)、さらに (エネルギー10keV以上の)非熱電子数(個)を 入力すると、放射束密度(SFU)・円偏波率が以下のように得られる。 Dulk (1985)の方法では
IDL> dulk_gysy,delta,bb,theta,nv,freq,fi,rc $<$CR$>$
さらに放射源の大きさ(立体角であたえる。単位はsterad)をあたえると 光学的厚さが求まる。
IDL> dulk_gysy,delta,bb,theta,nv,freq,fi,rc,omega,tau $<$CR$>$
同様にRamaty (1969)の方法では
IDL> ramaty_gysy,delta,bb,theta,nv,freq,fi,rc $<$CR$>$
IDL> ramaty_gysy,delta,bb,theta,nv,freq,fi,rc,omega,tau $<$CR$>$


(b) 放射から物理量を求める

Dulk (1985) の近似モデルにもとづいて、放射から放射領域物理量 を求めることができる。
$\alpha$をもとめたあと(節3.7参照)、 非熱電子の分布関数のべき指数がただちに求められる。
IDL> norp_alpha,freq,fi,mvd,mvdfit,alpha_tk,alpha_tn,freqpk,fluxpk $<$CR$>$
IDL> norp_alpha2delta,alpha_tn,delta $<$CR$>$



Takaaki YOKOYAMA 平成12年10月25日