C:ハヤシライン C: ハヤシライン
C.1.赤色巨星の構造 C.1.1.UV平面 星の構造は次の4本の方程式で与えられる ① ② ③ ④ C: ハヤシライン
その解を表すのによく使われるのが(U,V) 平面である。U,Vの定義は、 である。星の中心では、(U,V)=(0,3)、表面では(U、V)=(0,∞) となる。 したがって、星の構造は(U,V)面では左下のようになる。 化学組成に飛びがある赤色巨星のような場合には、右下のようになることは前回話した通りである。 表面 外層 V V 境界 中心核 境界 中心 3 0 3 0 U U 組成が一様な星のUVカーブ 組成に飛びがある星のUVカーブ C: ハヤシライン
C.1.2.ポリトロープ 前に述べた4つ構造方程式は5つの未知数M,P,T,ρ,Fを含んでいて式の数が不足している。それを与えるのは状態方程式で、P=P(T,ρ)であるが、この状態方程式がTを含まず、 と与えられると、構造方程式はM,P,ρに関する次の3つで閉じてしまう。 初めの2つの式から、ポアッソン方程式 が得られる。 中心密度=ρCとして、 ρ= ρC・θN と表わすと、 この式を上のポアッソン方程式に代入し、 C: ハヤシライン
この式はポリとロープ指数Nのエムデン方程式と呼ばれる。 等温の場合はN=∞に対応するが、初めに戻ってP=K・ρをポアッソン方程式に代入すると、 が等温のエムデン方程式である。 エムデン解から物理量に戻すには、パラメターとしてρC とαを選び、その2つからKを計算して、以下の式に入れればよい。 C: ハヤシライン
前ページの式を使って、U,Vを計算すると以下のようになる。 このように、中心から伸びるポリトロープのU-V曲線は中心の密度や圧力に関係なく、ポリトロープ指数Nだけで決まる。 N>5では表面にまで達する(U=0,V=∞)解を得ることはできない。 C: ハヤシライン
UVカーブでガス塊の構造を分類してみよう。 タイプ1: 正常型。自己重力系 C.1.3.重力系の分類 UVカーブでガス塊の構造を分類してみよう。 タイプ1: 正常型。自己重力系 自己の質量が重力の源となる。UV カーブで考えると、 星の表面: U=0 V=∞ 星の中心: U=3 V=0 となる。 タイプ1重力系のUVカーブ U 1 2 3 10 20 表面 中心 V 星の表面: U=0 V=∞ 星の中心: U=3 V=0
タイプ2: 質点の周りに広がるガス。 タイプ2のUVカーブ タイプ2: 質点の周りに広がるガス。 中心に質点があり、その重力場の影響のもとで広がるガス。中心から遠くではタイプ1と似た分布であるが、中心に近づくと質点重力場が支配的となる。 表面: U=0 V=∞ 中心近くは、 なので、 U→0 タイプ2のUVカーブ U 1 2 3 10 20 表面 中心 V タイプ2で重要な性質は、 UVカーブ上の点を0と名付けるとき、点0が中心に向かうにつれ、 (R1/R0)∞ となることである。
実際の星ではどうだろうか? 表面 (1)主系列星 タイプ1 UVカーブは中心から外側へ伸びる単純なタイプ1型である。 2 実際の星ではどうだろうか? (1)主系列星 タイプ1 表面 1 UVカーブは中心から外側へ伸びる単純なタイプ1型である。 log V 0 太陽の内部構造をUVカーブで表わしてみるとタイプ1型になっているのが分かる。 中心 -1 log U
中心核の重力場が外層には質点重力場のように働いていることが分かる。 中心核表面 (2)赤色巨星=複合型 log U log V -1 -2 1 3 2 星表面 中心核+外層 という2重構造 中心核UVカーブ=タイプ1 外層UVカーブ=タイプ2 左下のU=0、V=N+1へ向かう。 中心核の重力場が外層には質点重力場のように働いていることが分かる。 中心核表面 中心核表面 中心 このような構造になっていると、なぜ星の半径は主系列に比べ大きくなるのだろうか? それは、赤色巨星の中心にあるのが半径ゼロの質点でなく、有限の半径を持つ縮退核だからである。
AE AI BE BI 赤色巨星の半径は進化と共に増加する 中心核の半径Ro ≒白色矮星R0半径でほぼ一定。 外層の半径 R1 = Ro ・( R1 / Ro ) 進化が進むと中心核の質量が増加する。その結果、外層と核の境界が例えばA点からB点に移る。外層の解はタイプ2であり、UVカーヴが内側に進行すると(R1/R0)∞に増大する。こうして、星の半径R1は大きくなっていくのである。 赤色巨星の構造 赤色巨星構造の進化 1 0.5 1.5 log V 1.5 表面(R=R1) AE AI 表面(R=R1) 1 (R1/ R0)B> (R1/ R0)A log V BE BI 中心核表面 0.5 外層下面(R=R0) log U -1 -2 -2 -1 log U
無次元量でなく、実際の密度、半径で考えると、 中心核半径はあまり変わらないので、質量の増加は中心集中度を高めることで吸収している。その結果、下図に見るように、AからBに進むと中心密度は増加するが、核と外層の境界密度、圧力は逆に低下する。低い圧力でも支えられるようにするには外側が広がって重量(質量でなく)を軽減する必要がある。 ρ A B R R1(A) R1(B)
赤色巨星構造の中間まとめ これまで見てきたように、赤色巨星の構造の特徴は (1) 中心核と外層との境界で平均分子量がジャンプする。(μC/μE)>1 (2) 中心核の質量が増加していくと、(端の密度/平均密度)0となる。(U0) (3) しかし、中心核の半径R0はあまり変化しない。 (4) (2)で実現された(端の密度/平均密度)0は、外層にとってはタイプ2型の構造でR0 0の場合に対応する。その時はR/R0 ∞である。 (5) このため、R=R0 ・(R/R0 ) ∞ そのままだと、R ∞、Te0になってしまう。 このプロセスに歯止めをかけるのが、次に話す赤色巨星大気の表面条件である。 さらに詳しい議論は以下を参照するとよい。 Sugimoto,D.,Fujimoto,M.Y. 2000, ApJ 538, 837-853. R R0 μC μE
C.2.Hayashi Line C.2.1.赤色矮星の表面条件 C.2.Hayashi Line C.2.1.赤色矮星の表面条件 前回、Hoyle,Schwarzchild 1955 が赤色巨星の構造を考える際には表面条件が重要であると述べていたことを学んだ。低温度星の表面条件の問題は当初赤色矮星から始まった。 L.Bierman 1935、Astronomische Nachrichten, 257, 269 “Konvektion im Innern der Sterne” 中心から表面まで対流平衡な星では、P=Kργー1 が成立する。 表面でP=0,ρ=0という境界条件をおくと、Kが定まらない。したがって、このよ うな星は存在しないか、少なくとも中立平衡である。 これに対し、次のような反論が加えられた。 T.G.Cowling 1938、 Monthly Notices of the royal Astronomical Society,98, 734 “The Stability of the Convective Stars” 上の議論は表面近く、少なくともτ<1では、輻射でエネルギーを宇宙空間に放射 になっていることを見逃している。表面条件を考えると次の関係式が成立する。
(1) 光球(Photoshere)条件、T(τ=2/3)=Te 表面(r=R)付近で、 (1) 光球(Photoshere)条件、T(τ=2/3)=Te 表面(r=R)付近で、 を仮定すると、 (2) L=4πσTe4R2 (3) T=Kργー1 のポリトロープなので、 (4) 赤色矮星の場合さらに質量光度関係 L=A・Mαが成立する。 結局、M、R、L,Teff、Kに対し、4つの関係式が得られる。その結果、与えられたMに対してKが定まるのである。 C: ハヤシライン
このようにして、赤色矮星では表面条件を正しく扱う必要が明らかにされた。 Osterbrock1953は対流層を特徴付けるのは層の単位質量当たりエントロピーであることに着目して一連のモデル計算を行った。 D.E.Osterbrock 1953、Astrophys.J.118, 529-546 “The Internal Structure of Red Dwarf Stars” M型連星Castor Cは平均値、M=0.60Mo, L=0.063Lo, R=0.63Roと求まっていた。しかし、輻射層を仮定してしてモデルを作ると、光度が観測値よりはるかに明るくなることが問題とされていた。Stremgren(1952)は赤色矮星に対流層が発達していると温度勾配が小さいので中心温度が低くなって光度が下がるのではないか、と提案した。 Castor A,B,Cはそれぞれが二重星で、全体で六重星と分っている。Castor C はフレア星YY Gem として有名である。 右図は赤外Jバンド、A,B星は明るすぎて飽和している。左図はX線、C星以外にA,B星もX線で明るいことに注意。 XMMーニュートン(X線) 2MASS Jバンド(5分角) C: ハヤシライン
線形近似ではφ(τ)=τ+2/3 となる。κはHー吸収の平均値。 表面 灰色大気モデルを適用する。 輻射層 線形近似ではφ(τ)=τ+2/3 となる。κはHー吸収の平均値。 灰色大気の温度勾配が対流不安定になったら、対流層開始。 まず、開始点のエントロピーを計算。 S=一定のラインが電離層を通過し、logT=6,logP=13.8でほぼ完全電離(1-x)<<1。P=KT2.5の関係(logK=-1.15)が成立する。その先は下の無次元量で計算する。 対流層 E=13.6となる。色々なE(=8-40)に対して計算しておく。 吸収係数: κ=κ0ρ0.5T-3.5 エネルギー発生: ε= ε0ρT4.5 (p-pサイクル) 無次元量の計算では中心部での dlogP/dlogT=(n+1)Cがパラメターになる。 輻射核 C: ハヤシライン
対流層の解でE≒45が(3,0)を通ることを注意しておく。 鎖線:対流層で数字はEの値。 実線:輻射核で数字は(n+1)Cの値 点線:(n+1)=2.5となる実線の 端を結んだ線。 対流層と輻射核は点線の上でつながる。 その中で、CastorCは前に求めた E=13.6が解になるべきだが、 水素量Xから決まる(n+1)Cも適当 な値になるという条件を考えて、少し 異なる値 E=19.9、X=0.7,Y=o.26 を採用した。 対流層 n+1=2.5 16 24 32 E=40 対流層の解でE≒45が(3,0)を通ることを注意しておく。 C: ハヤシライン
C.2.2.赤色巨星の表面条件 前回、Hoyle,Schwarzchild 1955 が赤色巨星の構造を考える際には表面条件が重要であると述べていたことを学んだ。この問題をシステマティックに扱ったのが1961年の林忠四郎と蓬茨霊運の論文である。 Hayashi,C.,Hoshi,R. 1961,Publ.Astron.Soc.Japan,13,442-449. “The Outer Envelope of Giant Stars with Surface Convection Zones” 彼らはOsterbrock(1953)と同様に、星の表面付近の構造を左図のように考えた。 P=光球(Photosphere) PC=表面輻射層 C=輻射フラックスが対流フラックスと等しくなる点 CD=不完全電離層の断熱線(S=一定) DE=完全電離ガスの断熱線 P=KT2.5
重力g=GM/R2 と フラックスF は一定として、構造方程式を解く。 (1) P-C 表面輻射層 重力g=GM/R2 と フラックスF は一定として、構造方程式を解く。 表面近くは輻射層である。 Te=0の解を漸近輻射構造線と呼ぶ。表面から奥に入るとTeによらずこの解に収束していく。
(2). C-D:輻射層の奥に入ると対流不安定になる。その先は エントロピー一定で表わされる対流構造線になる。 エントロピー一定で表わされる対流構造線になる。 対流構造線 輻射構造線 漸近輻射構造線
対流不安定となるC点は通常電離層の中にある。 電離度xの等高線を示す。T=1万度付近で電離が起こることが分かる。
エントロピー(S/k)の等高線。電離層のところで傾きが変化。
星の構造線を電離度x、エントロピーS/k等高線に重ねてみると logTe=3.5S/k=40、logTe=3.45S/k=35と対応している。 Teの違いによって運ぶべきエネルギーが変わり、対流開始点に差が生じる。 対流開始
前ページの図では、M=1Mo,L=1000Loの星の場合、 有効温度 logTe 3.45 3.50 対流エントロピー(S/k) 35 41 (3) D-E: 完全電離。対流層 前ページの図では、M=1Mo,L=1000Loの星の場合、 有効温度 logTe 3.45 3.50 対流エントロピー(S/k) 35 41 この先は完全電離となるので、P=K・T2.5の関係が成立する。ここで、対流の エントロピー(S/k)とポリトロピックな関係の定数Kとの関係を求めておこう。 エントロピー(S/k)は下の式で定義していた。 このままでは完全電離x=1の場合が計算できない。 そこで、サハの式を用いてこの式を書き換える。 を一番上の式に代入して、 C: ハヤシライン
次にKの方だが、単純のため水素のみのガスを考える。 これで、x=1も計算できるようになった。 次にKの方だが、単純のため水素のみのガスを考える。 このKを上の(S/k)の式に代入して、 一方、圧力P、温度Tの無次元化 を行うと P=KT2.5は ポアッソン方程式 は、 となるので、 とおいて、エムデン方程式 を得る。 C: ハヤシライン
この方程式の解の様子は、Osterbrock1955のところで見た。 E>45.53だと空洞解となってしまう。これはKが大きいすなわち(S/k)が小さいことを意味する。つまり、表面温度が低すぎて対流層のエントロピーが小さくなると、圧力Pを低温T高密度ρで作ることになる。このため、同じ圧力差ΔPを作るためにより多くの質量Mを消費し、結局中心にたどりつく前に星の質量Mを全て使い尽くしてしまうのである。 E=45.53はn=1.5のエムデン解に対応し、中心まで対流となっている解である。また、 E<45.53は中心まで行ってもまだ使いきれない質量が残っている場合で、この残った質量は中心に置かれた質点と解釈される。これはタイプ2型の解で赤色巨星に対応する。 この様なわけで、HR図上にE=45.53のラインを引くとそのラインより低温側には星が存在できないことが分かった。このラインをハヤシラインと呼ぶ。 ハヤシラインは星の質量Mとメタル量Z毎に引かれる。 C: ハヤシライン
S/k=50、中心に空洞が出来てしまう。対応する天体はない。 ここではn=1.5の代わりにn=2.7の図を示すが様子は同じである。 S/k=60、中心に芯がある、タイプ2。赤色巨星に対応する。 S/k=55、タイプ1とタイプ2の境界線。HAYASHI lineの星。
こうして表面温度が低すぎる星(膨れ上がりすぎた星)はエントロピーが低く、中心部に空洞が生じて存在不可能であることが判った。そのぎりぎりの境界線をHAYASHI LINEと呼ぶ。 1Mo星の 禁止領域