2011/08/09

SPEX cie codeで計算に含まれる元素

SPEXのプラズマモデルはZ=1-30までの元素について組成比をパラメータとして持っていますが、実際のところCr, Mnなどは組成をどれだけ増やしてもHe-like Kα輝線とかが出てくることはありません。なぜならパラメータがあるだけで計算してないから!組成や電離度をちゃんと計算している元素がどれかは、ascdumpコマンドを叩けばわかります。

SPEX> com cie
 You have defined    1 component.
SPEX> calc
SPEX> ascd term 1 1 abund
 number of layer lines :        4076           0
 Atomic   Element  Abundance      Abundance     Average charge
 number   name     (solar units)  (absolute)
 -------------------------------------------------------------
  1        H       1.000          1.000             1.000
  2        He      1.000         9.7051E-02         2.000
  6        C       1.000         2.7733E-04         5.999
  7        N       1.000         8.1658E-05         6.996
  8        O       1.000         6.0534E-04         7.988
 10        Ne      1.000         1.2677E-04         9.923
 11        Na      1.000         2.2233E-06         10.82
 12        Mg      1.000         3.9719E-05         11.60
 13        Al      1.000         3.2584E-06         12.25
 14        Si      1.000         3.8548E-05         12.81
 16        S       1.000         1.6218E-05         14.21
 18        Ar      1.000         3.5727E-06         15.99
 20        Ca      1.000         2.3281E-06         17.87
 26        Fe      1.000         3.2659E-05         20.60
 28        Ni      1.000         1.8880E-06         22.38
てことで、N, Na, Alをのぞくodd Zの元素 (含 Mn)、Be, Ti, Crは計算してくれませんのでお気をつけを。

ちなみにlast columnのAverage chargeは楽しいですね。rtとか変えればもちろん値が変わっていきます。

0 件のコメント:

コメントを投稿