[論文レビュー] Construction of Highly Accurate Models of Rotating Neutron Stars - Comparison of Three Different Numerical Schemes

Tetsuo Nozawa, Nikolaos Stergioulas+2 more|CERN Bulletin|Apr 21, 1998

Geophysics and Gravity Measurements被引用数 4

ひとこと要約

本論文は、一般相対性理論における急速回転中性子星の高精度モデルを構築するための3つの数値コードを比較している。同一の状態方程式、補間手法、積分領域を用いて、大多数のモデルで0.1%未塔の一致を達成しており、これは、多相転移のないN > 0.5および現実的な状態方程式に対する信頼性を裏付け、このようなモデルのための標準的な物理的パラメータセットを確立している。

ABSTRACT

We conduct a direct comparison of three different representative numerical codes for constructing models of rapidly rotating neutron stars in general relativity. Our aim is to evaluate the accuracy of the codes and to investigate how the accuracy is affected by the choice of interpolation, domain of integration and equation of state. In all three codes, the same physical parameters, equations of state and interpolation method are used. We construct 25 selected models for polytropic equations of state and 22 models with realistic neutron star matter equations of state. The three codes agree well with each other (typical agreement is better than 0.1% to 0.01%) for most models, except for the extreme assumption of uniform density stars. We conclude that the codes can be used for the construction of highly accurate initial data configurations for polytropes of index N > 0.5 (which typically correspond to realistic neutron stars), when the domain of integration includes all space and for realistic equations with no phase transitions. With the exception of the uniform density case, the obtained values of physical parameters for the models considered in this paper can be regarded as ``standard'' and we display them in detail for all models.

研究の動機と目的

• 急速回転中性子星の一般相対性理論におけるモデル構築において、3つの異なる数値コードの精度を評価すること。
• 補間手法、積分領域、状態方程式が数値的精度に与える影響を調査すること。
• さまざまな物理的パラメータの範囲において、コードが一貫した結果を出力するかを特定すること。
• 数値相対性理論シミュレーションの初期データとして使用可能な、高精度な中性子星モデルのベンチマークセットを確立すること。
• 特に一様密度星のような極端な場合において、コードが失敗する条件を同定すること。

提案手法

• 本研究では、軸対称かつ定常な回転中性子星のためのアインシュタイン場方程式を解くために、それぞれ異なる数値スキームを用いた3つの異なる数値コードを採用している。
• すべてのコードが同一の物理的パラメータ、状態方程式（ポリトロープおよび現実的）および補間技術を用いており、一貫性を確保している。
• 境界効果を最小限に抑えるために、積分領域を全空間にまで拡張している。
• 反復緩和法を用いて、相対論的流体力学および重力場方程式を同時に解いている。
• 中心密度および回転速度の範囲をカバーする25のポリトロープ状態方程式（N > 0.5）および22の現実的状態方程式について、モデルを構築している。
• 収束性および精度は、質量、半径、回転周波数、コンパクトネスといった主要な物理的パラメータを、すべての3つのコード間で比較することで評価している。

実験結果

リサーチクエスチョン

• RQ1 急速回転中性子星のモデルを構築する際、3つの異なる数値コードの精度はどのように比較されるか？
• RQ2 補間手法および積分領域の選択が、モデルの数値的精度にどの程度影響を与えるか？
• RQ3 N > 0.5のポリトロープ状態方程式と現実的状態方程式の両方において、コードの性能はどのように異なるか？
• RQ4 一様密度星のような極端な状況において、コードが発散する条件は何か？
• RQ5 これらのコードが生成するモデルは、数値相対性理論シミュレーションの初期データとして標準的かつ信頼できるとみなせるか？

主な発見

• ポリトロープ指数N > 0.5および現実的状態方程式を用いた大多数のモデルについて、3つのコードは0.1%未塔の一致を示している。
• 多相転移のない現実的状態方程式を用いたモデルでは、すべてのコード間で0.01%未塔の一致が達成されている。
• 唯一顕著な不一致が観察されたのは、一様密度星という極端な場合であり、その差は1%を超えていた。
• 積分領域が全空間をカバーし、状態方程式が滑らかで相転移がない場合には、コードが一貫した高精度の結果を出力している。
• すべてのモデルの物理的パラメータ（質量、半径、回転周波数、コンパクトネス）は信頼性が高く、将来のシミュレーションのための標準と見なせる。
• 本研究では、47のモデル（25のポリトロープ、22の現実的状態方程式）からなる高精度なベンチマークデータセットを確立し、数値相対性理論における初期データとしての使用に適している。