宇宙物理学と状態方程式

京都大学大学院理学研究科
中村卓史
２００５年１２月２６日理論懇シンポ1例１） 球対称な星の進化
状態方程式は基本的に良く判っている
圧力＝輻射圧＋ガス圧＋縮退圧
核反応も主要なものは良く判っている
対流によるエネルギー輸送＝？だったが太陽ニュートリノでOK?
輻射の輸送、輻射の物質による散乱吸収再放出等も主要な部分は良く判っている
結論として基礎方程式と数値解法が確立している、
その結果著者にあまりよらない結果が得られていて教科書に書かれている
例2） 回転していたり質量放出する星の進化
状態方程式は例１）と同じ角運動量の輸送の理論の基礎方程式や質量放出の基礎方程式がはっきりしていないので教科書に書くレベルには至っていない。
2例３） 宇宙
圧力＝輻射圧＋ダークエネルギーの負の圧力
ダークエネルギーの理論はない！！
観測が先に進んでいるが、独立な検証が必要！！
種々の方法の適用と提案
１） TypeIa 超新星、
２） WMAP,
３） Baryon Acoustic Oscilation
４） GRBの光度ーピークエネルギー関係
5) 合体1年前の連星中性子星

34

5

中性子星の発見
ベルが最初に見た電波パルサー

P=1.3373011秒

6観測時間（秒）現在約１７００個 周期は１．４ｍｓから 数秒
7• 重力が遠心力以上でないと、天体はバラバラになってしまう。
• 未知の天体の質量をＭ、半径をＲとしてみよう• 原子核の密度• 星を作っている普通の物質
陽子：中性子：電子＝１：１：１• 普通の密度では電子は原子核のはるか外にあるが、このような高密度では電子と陽子が接するので電子捕獲が進んで中性子過剰物質となる。
89

10

dP/dtB
Age

BinaryP
11

Millisecond Pulsar

Binary12

Spiral structure13

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

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

中性子星の最大質量• ＴＯＶ方程式が正しい• 物質の音速は光速以下である。（因果律）• 核密度以下での状態方程式を知っている。
• 中性子星の最大質量は重力平衡解はない
１）ブラックホールの存在証明に状態方程式の詳細は不要
２）観測された中性子星の質量より最大質量が小さい状態方程式はダメである。
16

ブラックホール候補X線星
太陽質量の16倍

X線星

普通の星

太陽質量の6倍太陽質量の9倍太陽質量の9倍1718個のブラックホール候補天体

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

20

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

• 約１５年間のＰＳＲ１９１３＋１６からの電波パルスの到着時間のデータ解析から

22

23

伴星の質量主星の質量24

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

25

Shapiro Time delay
Orbital phase

近日点離角

Binary Pulsar Casers26

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

太陽質量

•

太陽質量

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

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

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

2004

J1906+0746

144

40.085

2.6

0.144

300

29

30

31

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

32

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

33

軌道周期の減少とSin i
近日点移動Mch以下

Sin I
進化理論

34

Projected semi major axis

35

36

37

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

39

40

基礎方程式は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)

42

43

このあたりが核密度 44

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

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

状態方程式を重力波から得る
中性子星の状態方程式
G=c=1の単位系

46

連星中性子星の軌道の変化
現在３億年後47

最後の３分間• チャープ信号

ブラックホール

重力波の振幅

約２５００回転
•
•

２０Ｈｚ←

３分間

高精度な理論予測が必要
ｋＨｚ→48連星中性子星の合体のシミュレーション（密度）49ＴＯＶ方程式(1939)• 一般相対論的平衡解(Tolman OppenheimerVolkov 方程式)

50

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

51

中性子星
白色矮星

52

エンタルピー ｈ
TOV方程式53

We can extend the equation of the state up to54

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

55

56

57

58

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

状態方程式に依存しないが、振動数は高い

59

物の大きさを知るには何かを当てるか破壊する。
60

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

61

62

63

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

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

64一般相対論でのケプラー問題
円軌道の角運動量依存性

重力波放出で角運動量が減少⇒円軌道の半径が減少⇒重力波の振動数増大
最小円軌道の後はほぼ自由落下65中性子星の状態方程式は全く不明？？？
１）QCDにもとずいた計算が完成するのを待つか？
２）重力波を使って質量と半径の関係を決めて逆問題を解くか？
３）両者が同時進行するのがベスト！！

66