質量光度関係のソースを表示

'''質量光度関係'''<ref name="kotobank"/><ref name="astrodic"/> ({{lang-en-short|mass–luminosity relation}}<ref name="astrodic"/>) は、[[恒星]]の[[質量]]と[[光度 (天文学)|光度]]を結び付ける式である。[[:en:Jakob Karl Ernst Halm|Jakob Karl Ernst Halm]] によって初めて記述された<ref>{{cite journal |last=Kuiper |first=G.P. |date=1938 |title=The Empirical Mass-Luminosity Relationship |journal=Astrophysical Journal |volume=88 |pages=472–506 |bibcode = 1938ApJ....88..472K |doi = 10.1086/143999 }}</ref>。恒星の質量を<math>M</math>、光度を <math>L</math> とすると、質量光度関係は以下の式で表される。

:<math>\frac{L}{L_{\odot}} = \left(\frac{M}{M_{\odot}}\right)^a</math>

ここで、<math>{L_{\odot}}</math> と <math>{M_{\odot}}</math> はそれぞれ[[太陽光度]]と[[太陽質量]]であり、a は 1〜6 の値を取る<ref name="evolutionofstars">{{cite book |last=Salaris |first=Maurizio |author2=Santi Cassisi |title=Evolution of stars and stellar populations |publisher= John Wiley & Sons |date=2005 |pages=138–140 |isbn=978-0-470-09220-0 |url=https://books.google.com/?id=r1dNzr8viRYC&lpg=PA138&dq=Mass-Luminosity%20relation&pg=PA138#v=onepage&q=&f=false}}</ref>。[[主系列星]]に対しては一般的に a = 3.5 が用いられる<ref>{{cite web |url=http://hyperphysics.phy-astr.gsu.edu/hbase/Astro/herrus.html#c3 |title=Mass-luminosity relationship |publisher=Hyperphysics |accessdate=2009-08-23}}</ref>。この関係式と一般的な a = 3.5 という値は 2〜55太陽質量を持つ主系列星には当てはまるが、[[赤色巨星]]や[[白色矮星]]には当てはまらない。また恒星の光度が[[エディントン光度]]に近づくにつれ、a = 1 に近づく。

== 概要 ==
指数の a の値は恒星質量の範囲によって異なる値をとり、それぞれの範囲に対して以下の式でよく近似できる<ref name="evolutionofstars"/><ref name="advanced">{{cite book |last= Duric |first= Nebojsa |title=Advanced astrophysics |publisher= Cambridge University Press |date=2004 |pages=19 |isbn=978-0-521-52571-8 |url=https://books.google.com/?id=-ljdYMmI0EIC&lpg=PA19&dq=Mass-luminosity%20relation&pg=PA19#v=onepage&q=&f=false}}</ref><ref>{{Cite web |url=http://jila.colorado.edu/~pja/astr3730/lecture18.pdf |title=The Eddington Limit (Lecture 18) |last= |first= |date=2003 |website=jila.colorado.edu |archive-url= |archive-date= |dead-url= |accessdate= 2019-01-22}}</ref>。

:<math>\frac{L}{L_{\odot}} \approx 0.23\left(\frac{M}{M_{\odot}}\right)^{2.3}   \qquad (M < 0.43M_{\odot})</math>
:<math>\frac{L}{L_{\odot}} = \left(\frac{M}{M_{\odot}}\right)^4   \qquad\qquad      (0.43M_{\odot} < M < 2M_{\odot})</math>
:<math>\frac{L}{L_{\odot}} \approx 1.4\left(\frac{M}{M_{\odot}}\right)^{3.5}   \qquad (2M_{\odot} < M < 20M_{\odot})</math>
:<math>\frac{L}{L_{\odot}} \approx 32000 \frac{M}{M_{\odot}}   \qquad \qquad (M > 55M_{\odot})</math>

0.43太陽質量よりも軽い恒星では内部での[[エネルギー]]の輸送は[[対流]]のみが担っており、そのため質量光度関係が大きく変わる。55太陽質量よりも重い恒星では質量と光度の関係は平坦なものとなり <math>L\propto M</math> の比例関係に近づくが<ref name="evolutionofstars"/>、このような大質量星は不安定であり強い[[恒星風]]によって急速に質量を失うため、この状態は長くは継続しない。この変化は大質量星での[[放射圧|輻射圧]]の増加によるものである<ref name="evolutionofstars"/>。これらの関係式は、[[年周視差]]や[[宇宙の距離梯子|その他の方法]]によって距離が判明している[[連星]]の恒星質量を測定することによって、経験的に導出されている。十分な数の恒星のデータを集めて質量と光度を[[両対数グラフ]]上にプロットすることで線を描くことができ、その傾きを測定することで適切な a の値を決定することが出来る。

[[K型主系列星]]に対して有効な別の表式もあり、これは指数 a の不連続性を回避するために与えられたものである<ref>{{cite journal |last=Cuntz |first=M. |author2=Wang, Zh. | date=2018 |title=The Mass–Luminosity Relation for a Refined Set of Late-K/M Dwarfs |journal=Research Notes of the American Astronomical Society |volume=2 |issue=1 |pages=19 |bibcode = 2018RNAAS...2a..19C |doi = 10.3847/2515-5172/aaaa67}}</ref>。その関係式は、

:<math>\frac{L}{L_{\odot}}  \approx \left(\frac{M}{M_{\odot}}\right)^{a(M)}   \qquad\qquad      (0.20M_{\odot} < M < 0.85M_{\odot})</math>
:<math>a(M) = -141.7 \, M^4 + 232.4 \, M^3 - 129.1 \, M^2 + 33.29 \, M + 0.215</math>

である。ただしここでは、<math>M</math> は太陽質量 <math>M_{\odot}</math> を単位とした値を用いる。この関係式は、年周視差が判明しており、[[有効温度]]と光度を正確に測定するために恒星半径が[[干渉法|干渉計]]によって測定されている晩期K型星と[[赤色矮星|M型主系列星]]の中分散[[スペクトル]]を使用したデータに基づいて導出されたものである<ref>{{cite journal |last=Mann |first=A. W. | author2=Gaidos, E. |author3=Ansdell, M. |date=2013 |title=Spectro-thermometry of M Dwarfs and Their Candidate Planets: Too Hot, Too Cool, or Just Right? |journal=Astrophysical Journal |volume=779 |issue=2 |pages=188 |bibcode = 2013ApJ...779..188M |doi = 10.1088/0004-637X/779/2/188 |arxiv=1311.0003}}</ref>。これらの恒星は[[ケプラー (探査機)|ケプラー]]での候補天体を較正する際にも用いられている。この関係式は <math>M=0.43M_{\odot}</math> での指数の不連続性を回避する以外に、<math>M=0.85M_{\odot}</math> の a = 4.0 の値もカバーしている。

質量光度関係は、通常の[[年周視差]]の測定には遠すぎる[[連星]]の距離測定に使えるため重要である。この技術は{{仮リンク|力学視差|en|dynamical parallax}}<ref>{{Cite web|和書| url = https://astro-dic.jp/dynamical-parallax/ | title = 天文学辞典 » 力学視差 | author = | authorlink = | coauthors = | date = | format = | work = 天文学辞典 | publisher = 日本天文学会 | pages = | language = | archiveurl = | archivedate = | quote = | accessdate = 2019-03-26}}</ref>と呼ばれている<ref name="double">{{cite book |last=Mullaney |first=James |title=Double and multiple stars and how to observe them |publisher= Springer |date=2005 |isbn=978-1-85233-751-3 |url=https://books.google.com/?id=OM3hozlKt9AC&lpg=PA27&dq=Mass-Luminosity%20relation%20distance%20binary&pg=PA27#v=onepage&q=&f=false}}</ref>。この技術では、連星系にある2つの恒星の質量が、一般にはまず太陽程度の質量であると推定される。次に[[天体力学]]の[[ケプラーの法則]]を用いて連星間の距離を計算する。この距離が分かれば、連星間の距離が天球上で描く弧の長さを介して連星系までの距離を暫定的に推定することが出来る。この測定と連星の両方の[[等級 (天文)|見かけの等級]]から光度を求めることができ、さらに質量光度関係を用いることでそれぞれの恒星の質量を求めることが出来る。こうして導出された質量を用いて連星間の距離を再計算し、これまでの計算を再び繰り返す。一連の計算を何度も反復することで、最高で 5% の精度でパラメータを決定することができる<ref name="double"/>。恒星の寿命がおおむね <math>M/L</math> に比例することに注目して、恒星の寿命を決定する場合にも質量光度関係が用いられる場合がある。ただし大きな質量を持つ恒星の場合は、質量光度関係から推定した値よりも短い寿命になることが知られている。より洗練された計算では、時間の経過とともに恒星が質量を失う効果が含まれる。

== 導出 ==
理論的に厳密な質量光度関係の導出には、エネルギー生成方程式と、恒星内部の熱力学的モデルの構築が必要となる。しかしいくつかの基礎的な物理と単純化した仮定を課すことで、基本的な関係式である <math>L \propto M^{3}</math> を導出することが出来る<ref name="physicsofstars">{{cite book |last=Phillips |first=A.C. |title=The Physics of Stars |publisher=John Wiley & Sons |date=1999 |isbn=978-0-471-98798-7}}</ref>。この導出は[[1924年]]に天文学者の[[アーサー・エディントン]]によって初めて行われた<ref>{{cite book |last=Lecchini |first=Stefano |title=How Dwarfs Became Giants. The Discovery of the Mass-Luminosity Relation |isbn=978-3-9522882-6-9 |publisher=Bern Studies in the History and Philosophy of Science |url=http://philoscience.unibe.ch/bern-studies/how-dwarfs-became-giants-the-discovery-of-the-mass-luminosity-relation |year=2007}}{{dead link |date= 2018-01 |bot=InternetArchiveBot |fix-attempted=yes }}</ref>。この導出の際、恒星は[[理想気体]]で近似的にモデル化出来ると仮定されたが、これは当時としては斬新で、幾分か過激なアイデアですらあった。以下に示すのは、同じ原理に基づいているがより現代的なアプローチによる導出である。

恒星の光度 (単位時間あたりに放射される[[エネルギー]]) を決める重要な要素は、恒星全体でのエネルギー散逸率である。[[対流]]が存在しない場合、散逸は主に[[光子]]の[[拡散]]によって起きる。対流が無視できる領域である[[放射層]]におけるある半径 r での表面積にわたって[[フィックの法則]]を積分することで、外向きに流れる総エネルギー流を計算することが出来る。[[エネルギー保存の法則]]より、このエネルギーのフラックスは恒星の光度と等しくなる。

::<math>L = -4\pi\,r^2 D\frac{\partial u}{\partial r}</math>

''D'' は光子の拡散係数、''u'' はエネルギー密度である。

ここで、恒星は全体が対流しているわけではなく、エネルギーを生成する過程 ([[核融合反応]]) は放射層よりも内部にある恒星の[[核 (天体)|核]]で起きていることを仮定している。これら2つの仮定は[[赤色巨星]]では正しくなく、そのため赤色巨星は通常の質量光度関係には従わない。低質量の恒星も恒星全体が[[対流層]]となっており、この法則には従わない。

恒星が[[黒体]]であると仮定すると、エネルギー密度は[[シュテファン＝ボルツマンの法則]]によって温度と結び付けられる<ref name="astrophys"/>。

::<math>u = \frac{4}{c} \sigma_{\rm B} T^{4} </math>

ここで

::<math>\sigma_{\rm B} = \frac{2\pi^{5}k_{\rm B}^{4}}{15c^{2}h^{3}}=\frac{\pi^{2}k_{\rm B}^{4}}{60\hbar^{3}c^{2}}</math>

は[[シュテファン＝ボルツマンの法則#シュテファン=ボルツマン定数|シュテファン=ボルツマン定数]]である。また ''c'' は[[光速]]、<math>k_{\rm B}</math> は[[ボルツマン定数]]、<math>h</math> と <math>\hbar</math> はそれぞれ[[プランク定数]]と[[ディラック定数|換算プランク定数]]である。

気体中の拡散の初等理論では、拡散係数 ''D'' は次の関係で近似される。

::<math> D = \frac{1}{3}c \lambda </math>

ここで &lambda; は光子の[[平均自由行程]]である。

恒星の核や、温度が恒星の核と同程度になる領域では、物質は完全に[[電離]]している。そのため光子は主に[[電子]]と衝突し、平均自由行程 &lambda; は以下の関係式を満たす。

::<math> \lambda = \frac{1}{n_{e}\sigma_{e\cdot\gamma}}</math>

ここで <math>n_e</math> は電子の[[数密度]]、

::<math>\sigma_{e\cdot\gamma} = \frac{8\pi}{3}\left(\frac{\alpha\hbar c}{m_{e} c^{2}}\right)^{2}</math>

は電子-光子散乱の[[散乱断面積]]であり、[[トムソン散乱|トムソン断面積]]に等しい。また &alpha; は[[微細構造定数]]、''m''<sub>e</sub> は[[電子の静止質量|電子の質量]]である。

恒星の平均電子密度は恒星の質量 ''M'' と半径 ''R'' と関係しており、以下のように書き表される。

::<math>\langle n_{e} \rangle = \frac{M}{m_{e} 4\pi R^{3}/3}</math>

最後に、[[ビリアル定理]]より総[[運動エネルギー]]は[[位置エネルギー|重力ポテンシャルエネルギー]] ''E''<sub>G</sub> の2分の1に等しいという関係を用い、また原子核の平均質量を ''m''<sub>n</sub> とすると、原子核1個あたりの平均運動エネルギーは以下の関係式を満たす。

::<math>\frac{3}{2}k_{\rm B} T = \frac{1}{2}E_{\rm G}\frac{m_{\rm n}}{M} = C\frac{3}{10}\frac{G M m_{\rm n}}{R}</math>

ここで、温度 ''T'' は恒星全体で平均した値である。また ''C'' は恒星の構造によって変化する 1 程度の大きさを持つ係数であり、[[ポリトロープ]]を用いて恒星を近似することで推定することが出来る。この関係式は、放射層において[[放射圧]]がガス圧を上回るような大質量星においては成り立たず、温度と質量、半径の関係は異なったものになる。これに関しては後述する。

全ての式をまとめ、またある因数で ''r'' と ''R'' が等しくなるように ''r'' を取り、さらに ''r'' での ''n''<sub>''e''</sub> をある因数で恒星全体での平均値と置き換える近似を行う。太陽の場合、2つを合わせた因数は 1/15 となり、最終的に以下の関係を得る。

:<math>L = -4\pi r^{2} D\frac{\partial u}{\partial r} \approx 4\pi R^{2} D\frac{u}{R}</math>
::<math>L \approx  \frac{1}{15}\frac{64\pi^{2}}{9}\frac {\sigma_{\rm B}}{\sigma_{e\cdot\gamma}}\frac{R^{4} T^{4}}{M}</math>
::<math>\,\, \approx  \frac{1}{15}\frac{2\pi^{3}}{9 \cdot 5^{5}}\frac {G^{4} m_{e}^{2} m_{\rm n}^{5}}{\alpha^{2}\hbar^{5}}M^{3} = 4\cdot {10^{26}} \left(\frac{M}{M_{\odot}}\right)^{3}\,\left({\rm W}\right)</math>

ここで付加した因数は実際には ''M'' に依存しているため、より詳細な計算から導いた光度の質量に対する依存性は <math>M^{3.5}</math> となる<ref name="astrophys"/>。

== 小質量および大質量での区別 ==
{{Main|エディントン光度}}
[[放射圧]]を使って上記の結果を導出することによって、恒星の質量が小さい場合と大きい場合を区別した質量光度関係を導くことが出来る。ここでは、不透明度 &kappa; を用い、恒星の内部温度 ''T''<sub>I</sub> を直接考慮することで導出を簡単化できる。より正確に取り扱う場合は、[[放射層]]の平均温度が使われる。

まずは放射圧 <math>P_{\rm rad}</math> と光度の関係を記述する所から始める。放射圧の勾配は放射の吸収による[[運動量]]の輸送と等しくなることから、以下の関係式が成り立つ。

:<math>\frac{dP_{\rm rad}}{dr} = -\frac {\kappa\rho}{c}\frac{L}{4\pi r^2}</math>

ここで ''c'' は[[光速]]である。また <math>1/\kappa\rho = l</math> は光子の平均自由行程である。

放射圧と温度には以下の関係がある。

:<math>P_{\rm rad}=\frac{4\sigma}{3c}T_{\rm I}^{4}</math>

従って

:<math>T_{\rm I}^{3}\frac{dT_{\rm I}}{dr} = -\frac{3 \kappa \rho}{16 \sigma}\frac{L}{4 \pi r^{2}}</math>

となり、この式から以下の比例関係が導かれる。

:<math>L \propto T_{\rm I}^{4} \frac{R}{\rho} \propto T_{\rm I}^{4} \frac{R^{4}}{M}</math>

[[放射層]]では、[[重力]]はガス自身の[[圧力]] ([[理想気体]]の圧力で近似される) と放射圧の両方と釣り合う。質量が小さい恒星においては後者は無視できるため、先述のように

:<math>T_{\rm I} \propto \frac{M}{R}</math>

と書き表される。より正確には、積分は半径が 0 から ''R'' まで行うため、左辺は表面温度 ''T''<sub>E</sub> を用いて <math>T_{\rm I}-T_{\rm E}</math> となるが、表面温度は内部の温度 ''T''<sub>I</sub> と比較して無視できる。

この関係を用いると、質量光度関係

:<math>L \propto M^{3}</math>

を導出できる。

恒星の質量が十分に大きい場合、放射層では放射圧がガス圧を上回るようになる<ref name="astrophys"/>。ここまで用いた理想気体のガス圧の代わりに放射圧を代入することで、

:<math>T_{\rm I}^{4} \propto \frac{M^{2}}{R^{4}}</math>

が得られ、従って大質量星に対する質量光度関係

:<math>L \propto M</math>

が導かれる<ref name="astrophys"/>。

== 核と表面の温度 ==
第一近似として、恒星が表面積 <math>4\pi R^{2}</math> を持つ[[黒体]]の放射体であるとする。従って[[シュテファン＝ボルツマンの法則]]から光度と表面温度 ''T''<sub>S</sub> が結び付けられる。

:<math>L = 4\pi R^{2}\sigma_{\rm B} T_{\rm S}^{4}</math>

この光度は恒星内部で単位時間あたりに生成される総エネルギーと等しい。このエネルギーは[[核融合反応]]によって生成されているため、通常は恒星の核で発生しており ([[赤色巨星]]の場合はそれが成り立たない)、核の温度は単位体積あたりの核融合反応の反応率によって光度と結び付いている。

:<math>L = \frac{dE}{dt} \approx  \epsilon \frac{4\pi}{3} R^{3} n_{A} n_{B} \frac{4\sqrt{2}}{\sqrt{3 m_{R}}} \sqrt{E_{0}(T)} \frac{S(E_{0}(T))}{kT} e^{-\frac{3E_{0}(T)}{kT}}</math>

ここで &epsilon; は[[陽子-陽子連鎖反応]]もしくは[[CNOサイクル]]によって生み出される総エネルギーである。また ''E''<sub>0</sub> はガモフピークのエネルギーであり、{{仮リンク|ガモフ因子|en|Gamow factor}} ''E''<sub>G</sub> に以下のように依存する<ref name="astrophys"/>。

:<math>E_{0} = \left(\sqrt{E_{\rm G}}k_{\rm B}T\right)^{2/3}</math>

''S(E)'' は「天体物理学的S因子」と呼ばれる物理量であり、''S(E)/E'' が反応断面積になる<ref name="astrophys"/>。また、''n'' は数密度である。<math>m_{R} = m_{A} m_{B}/(m_{A} + m_{B})</math> は粒子衝突の[[換算質量]]であり、''A'' と ''B'' は制限反応に関与する2つの核種を表す。例えば陽子-陽子連鎖反応の場合は ''A'' と ''B'' はどちらも[[陽子]]であり、CNOサイクルの場合は ''A'' は陽子、''B'' は {{nuclide|nitrogen|14}} などとなる。

恒星の半径 ''R'' 自身も温度と質量の関数であるため、この方程式を解くことで核の温度を求めることが出来る。

== 出典 ==
{{reflist|refs=
<ref name="kotobank">{{Cite web|和書| url = https://kotobank.jp/word/è³ªé‡å…‰åº¦é–¢ä¿‚-521799 | title = 質量光度関係(シツリョウコウドカンケイ)とは - コトバンク | publisher = コトバンク | accessdate = 2019-03-26 }}</ref>

<ref name="astrodic">{{Cite web|和書| url = https://astro-dic.jp/mass-luminosity-relation/ | title = 天文学辞典 » 質量-光度関係（星の） | author = | authorlink = | coauthors = | date = | format = | work = 天文学辞典 | publisher = 日本天文学会 | pages = | language = | archiveurl = | archivedate = | quote = | accessdate = 2019-03-26}}</ref>

<ref name="astrophys">{{Cite book |和書 |last= |first= |author= 高原文郎 |authorlink= |coauthors= |translator= |year= 1999 |title= 宇宙物理学 |publisher= 朝倉書店 |page= |id= |isbn= 978-4-254-13076-8 |quote= }}</ref>
}}

== 関連項目 ==
* [[主系列星]]
* [[ヘルツシュプルング・ラッセル図]]
* [[恒星内元素合成]]
* [[陽子-陽子連鎖反応]]
* [[CNOサイクル]]

{{恒星}}

{{DEFAULTSORT:しつりようこうとかんけい}}
[[Category:天文学に関する記事]]
[[Category:恒星物理学]]