例１） 球対称な星の進化
状態方程式は基本的に良く判っている
圧力＝輻射圧＋ガス圧＋縮退圧
核反応も主要なものは良く判っている
対流によるエネルギー輸送＝？だったが太陽ニュートリノでOK?
輻射の輸送、輻射の物質による散乱吸収再放出等も主要な部分は良く判っている
結論として基礎方程式と数値解法が確立している、
その結果著者にあまりよらない結果が得られていて教科書に書かれている

宇宙物理学と状態方程式

京都大学大学院理学研究科
中村卓史
２００５年１２月２６日理論懇シンポ

例2） 回転していたり質量放出する星の進化
状態方程式は例１）と同じ角運動量の輸送の理論の基礎方程式や質量放出の基礎方程式がはっきりしていないので教科書に書くレベルには至っていない。
1234

例３） 宇宙圧力＝輻射圧＋ダークエネルギーの負の圧力
ダークエネルギーの理論はない！！
観測が先に進んでいるが、独立な検証が必要！！
種々の方法の適用と提案
１） TypeIa 超新星、
２） WMAP,
３） Baryon Acoustic Oscilation
４） GRBの光度ーピークエネルギー関係
5) 合体1年前の連星中性子星中性子星の発見
ベルが最初に見た電波パルサー

5P=1.3373011秒

61• 重力が遠心力以上でないと、天体はバラバラになってしまう。
• 未知の天体の質量をＭ、半径をＲとしてみよう• 原子核の密度• 星を作っている普通の物質
陽子：中性子：電子＝１：１：１• 普通の密度では電子は原子核のはるか外にあるが、このような高密度では電子と陽子が接するので電子捕獲が進んで中性子過剰物質となる。

観測時間（秒）
現在約１７００個 周期は１．４ｍｓから 数秒
78910dP/dt
B
Age

BinaryMillisecond Pulsar

P
11

Binary122

Spiral structure

ＴＯＶ方程式• 一般相対論的平衡解(Tolman OppenheimerVolkov 方程式)f >11314

中性子星の最大質量2.11.44• ＴＯＶ方程式が正しい• 物質の音速は光速以下である。（因果律）• 核密度以下での状態方程式を知っている。
• 中性子星の最大質量は重力平衡解はない

どの状態方程式が正しいのか全く不明

15

ブラックホール候補X線星

１）ブラックホールの存在証明に状態方程式の詳細は不要
２）観測された中性子星の質量より最大質量が小さい状態方程式はダメである。
1618個のブラックホール候補天体

太陽質量の16倍

X線星

普通の星

太陽質量の6倍太陽質量の9倍太陽質量の9倍1718

3

中性子星の質量の観測

中性子星の質量が４桁の精度で決まった連星パルサーＰＳＲ１９１３＋１６（１９９３年度ノーベル物理学賞）

19

連星の軌道パラメター

20• 約１５年間のＰＳＲ１９１３＋１６からの電波パルスの到着時間のデータ解析から• 軌道長半径 ～ ７０万ｋｍ（ほぼ太陽の半径）• 公転周期 ～ ８時間• 楕円 （離心率 ～ ０．６）
• 電波パルサー ＝ 正確な時計• 例 ＰＳＲ１９３７＋２１のパルス周期
１．５５７８０６４６８１９７９４ミリ秒• 連星パルサー ＰＳＲ１９１３＋１６は自然が与えて21くれた 「一般相対論の実験場」

伴星の質量

22

主星の質量
23

244

Shapiro Time delay
Orbital phase

光の曲がりは０．０２％の精度で確認されている。（Ｗｉｌｌ ２００１）
シャピロ時間遅れは０．１％の精度で確認されている。

近日点離角

Binary Pulsar Casers25

両中性子星の質量が正確に決まった
•

太陽質量

•

太陽質量

• それ以外の大発見として、公転周期が１年間に ７６．９４マイクロ秒短くなっていることが分かった。
• これは、何を意味するのか？⇒重力波の存在の間接的な証明• ケプラーの法則によると公転周期の２乗は軌道長半径の３乗に比例するので、両中性子星が近づいていることになる。

26• ＰＳＲ１９１３＋１６は連星系なので加速度運動をしている。⇒重力波が放射されている。
⇒重力波は正のエネルギーを持っているので連星系はエネルギーを失う。⇒連星は近づ
く⇒周期は短くなる• アインシュタインの一般相対性理論による理論計算によると１年間に７６．１５マイクロ秒公転周期が短くなるはずである。
• 先に得られた観測値との誤差はわずかに０．
１％程度である。
しかし、重力波が直接観測されたわけではない。

2728

2930

発見された連星中性子星のリスト

2004

J1906+0746

144

40.085

2.6

0.144

300

5

中性子星と白色矮星の連星質量降着時期があるので中性子星の質量が増大する可能性あり。

3132

軌道周期の減少とSin i

Keplerの第３法則から近日点移動
Mch以下３つの変数に１つの関係のみ
これ以外に、１） 近日点移動 ２）シャピロタイムディレー
３）重力波放出による公転周期の減少の２つが必要
Sin I
進化理論

Projected semi major axis

33

34

35

36

6

37383940

cos iは精度はあまりよくならないが、周期の変化は観測時間の2.5乗で良くなる。

基礎方程式はQCDの理論
陽子や中性子間に働いて原子核を結合させている強い力はQCD(Quantum Chromo Dynamics)で記述されることが確立している。
Quark F=6

グルオン場 a=1-8ゲルマン行列

SU(3)の構造定数41Quark u(up), d(down), s(strange), c(charm), b(bottom), t(top)427このあたりが核密度 4443

２つの問題

状態方程式を重力波から得る

１） 格子の数が少ない。⇒数値物理の初期には必ずそうなる問題なし。

中性子星の状態方程式

２） 温度が低いときの計算法が確立していない⇒数値物理の初期には必ずそうなる問題なし。１発のアイデアで解決する場合を数々見てきた。

G=c=1の単位系

例１） 数値相対論 初期は２８X28、計算の不安定性はConformal TT Decompositionで解決
例２） Newton N体問題
Multi-time step法、Tree法、Grape等の専用計算機
4546

最後の３分間連星中性子星の軌道の変化• チャープ信号

現在

ブラックホール

重力波の振幅３億年後

約２５００回転

47

•
•

２０Ｈｚ←

３分間

高精度な理論予測が必要
ｋＨｚ→

48

8

連星中性子星の合体のシミュレーション（密度）ＴＯＶ方程式(1939)• 一般相対論的平衡解(Tolman Oppenheimer
Volkov 方程式)

49

50

中性子星

逆問題

白色矮星

中心密度と状態方程式を与えると、中性子星の半径と質量はTOV方程式を中心から解くことによって決まる。
中性子星の半径が質量の関数として得られたら状態方程式は求まるか？
YES． Lindblom 1992

5152

エンタルピー ｈ
TOV方程式We can extend the equation of the state up to53

54

9問題点 たくさんの質量と半径のデータがいるが、それは可能か？
M,Rが1個しか得られない場合
では状態方程式が判っているとすると仮定

55565758

如何にして中性子星の質量と半径を決定するか？
質量は連星中性子星とブラックホールのチャープシグナルから決定できる？
１） 重力波のpower spectrum にform factorとして半径の逆数に比例した振動数の間隔でピークが表れる
Saijo and Nakamura 2001

状態方程式に依存しないが、振動数は高い物の大きさを知るには何かを当てるか破壊する。
5960

10構造を持つ場合構造を持たない点粒子

6162

２） ISCO（Inner Most Stable Orbit) での重力波の振動数が決定可能とする。 つまり、ISCO以内では急速に重力波の振幅が下がると考える。

ISCOは連星の各質量と状態方程式に依存する。
⇒半径、質量と同様のプログラムが可能？

63

一般相対論でのケプラー問題

64

中性子星の状態方程式は全く不明？？？
１）QCDにもとずいた計算が完成するのを待つか？

円軌道の角運動量依存性

２）重力波を使って質量と半径の関係を決めて逆問題を解くか？

重力波放出で角運動量が減少⇒円軌道の半径が減少⇒重力波の振動数増大

３）両者が同時進行するのがベスト！！

最小円軌道の後はほぼ自由落下65

6611