Optically-Thin Non-Thermal Gyrosynchrotron Emission (Preliminary)

Note: The procedures described in this subsection are under testing. Comments are welcome for improvement.

Based on approximation models, the relation between the physical variables of emitting region and the emission can be derived. The electron distribution function is assumed to be the power law.

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

where $E$ is electron energy (keV), $NV(E)$ is number of electrons (particles) that has larger energy than $E$. This distribution is described with parameters $\delta$, $E_{0}$, and $K$. But it is usual to use $NV(E_0)=K/(\delta -1)/E_0^{\delta-1}$ instead of $K$. And we fix $E_{0} = 10\ {\rm keV}$ after Dulk (1985).

(a) From Physical Variables to Emission

Inputs; power-law index $\delta$, magnetic field strength (G), angle between magnetic field and line of sight (degree), and $NV$ -- total number (particles) of non-thermal electron (of $E> 10\ {\rm keV}$): Outputs; flux density (SFU), circular polarization degree:
By Dulk's (1985) method:
IDL> dulk_gysy,delta,bb,theta,nv,freq,fi,rc $<$CR$>$
If the size of emitting source (by solid angle in unit of sterad) is given additionally, optical depth is also obtained
IDL> dulk_gysy,delta,bb,theta,nv,freq,fi,rc,omega,tau $<$CR$>$
By Ramaty's (1969) method:
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) From Emission to Physical Variables

Based on Dulk's (1985) model we may derive the physical variables from emission. After deriving $\alpha$ (see 3.7), power law index of non-thermal electron distribution function is
IDL> norp_alpha,freq,fi,mvd,mvdfit,alpha_tk,alpha_tn,freqpk,fluxpk $<$CR$>$
IDL> norp_alpha2delta,alpha_tn,delta $<$CR$>$