| 研究概要 |
これまでに申請者は,非加法的Tsallisエントロピーを用いて重力多体系の平衡状態の解析を行い,エントロピーの変分の極値から平衡状態を求め,圧力P,密度pに対する状態方程式がポリトロープ関係Pχp^<1+1/n,n=1/q-1+3/2に従うポリトロープ状態が存在することを示した.
ところで,重力多体系には,collision time t_cという粒子間の2体重力散乱による緩和の時間スケールが存在する.従って,上で求めたポリトロープ状態は,実際にはcollision time t_cで時間変化すると予想されていた.現在進行中の断熱壁に閉じこめたN体数値シミュレーションの結果では,系は,まずあるポリトロープ状態に向かって緩和し,その後はポリトロープ状態の系列を等温状態に向かって時間発展して行くことが示唆されている.
そこで今年度は,Tsallisエントロピーの極値で与えられるポリトロープ状態が重力多体系の準定常状態である,という理解を確立すること目的として研究を行った.具体的には,(1)より一般的な初期状態からの系統的な重力多体系のN体数値シミュレーションを実行し,ポリトロープ状態が準定常状態であることを確認下.さらに,まだ非常に特別な例しか得ていないが,ポリトロープ状態に漸近しない初期状態の系列を発見した.
(2)"比熱が負"という重力系の本質が具現されているcore-collapseという現象にも本研究の視点・手法を適用することを計画し,そのための数値的手法の準備を行った.ポリトロープ状態から始まるCollapseを解析するため,(1)一般化された変分原理によるFokker-Planck方程式の解法,(2)連星形成を取り扱えるN体計算のコード開発,を行った.
|
