[論文レビュー] MHD instabilities in accretion mounds on neutron star binaries

Dipanjan Mukherjee, D. Bhattacharya+1 more|arXiv (Cornell University)|Apr 26, 2013

High-pressure geophysics and materials被引用数 3

ひとこと要約

本研究では、プラズマコード「PLUTO」を用いた3次元MHDシミュレーションとGrad-Shafranov方程式の数値解法を用いて、中性子星連星系における降着マウンドの磁流力学的不安定性を調査した。圧力駆動型不安定性が、マウンド質量が約5×10⁻¹³ M☉の閾値を超えると発生し、指状構造を形成して長期間にわたる磁場埋没をダイアマグネティックスクリーニングによって制限することが判明した。

ABSTRACT

We have numerically solved the Grad-Shafranov equation for axisymmetric static MHD equilibria of matter confined at the polar cap of neutron stars. From the equilibrium solutions we explore the stability of the accretion mounds using the PLUTO MHD code. We find that pressure driven modes disrupt the equilibria beyond a threshold mound mass, forming dynamic structures, as matter spreads over the neutron star surface. Our results show that local variation of magnetic field will significantly affect the shape and nature of the cyclotron features observed in the spectra of High Mass X-ray Binaries.

研究の動機と目的

• 密度の高いフェルミプラズマの現実的な状態方程式を用いて、中性子星表面近くの高密度プラズマを正確にモデル化するため、相対論的フェルミガス状態方程式（p ∝ ρ⁵/³）を組み込んだ軸対称Grad-Shafranov方程式を用いて、降着マウンドの磁気的平衡状態をモデル化すること。
• 異なるマウンド質量および磁場曲率の下で、特に圧力駆動モードを含むMHD不安定性に対するこれらの平衡状態の安定性を調査すること。
• MHD不安定性が高質量X線連星系における観測されたサイクロトロン共鳴散乱特徴（CRSF）および中性子星磁場の長期的進化に与える影響を評価すること。
• MHD不安定性が支配的になるマウンド質量およびプラズマベータ（β）の閾値を特定し、磁束線の歪みおよびダイアマグネティックスクリーニングの制限要因を解明すること。

提案手法

• 相対論的フェルミガス状態方程式（p ∝ ρ⁵/³）を組み込んだ軸対称Grad-Shafranov方程式を解き、中性子星表面近くの高密度プラズマを正確にモデル化した磁気的平衡状態を求める。
• PLUTO MHDコードを用いて、摂動を加えた平衡解の3次元非軸対称シミュレーションを実施し、MHD不安定性の時間的成長を調査する。
• マウンド高さを磁束関数Z₀(ψ)の関数として定義し、GS解から得られる磁力線および圧力勾配を用いて、現実的なマウンド幾何形状をモデル化する。
• 時間経過に伴う密度および磁場強度の変化を追跡し、マウンドの径方向端縁に指状不安定性が発生するタイミングを特定する。
• 異なるマウンド質量におけるプラズマベータ（β）を計算し、理論的不安定性閾値（例：Litwin et al. 2001）と比較する。
• マウンド質量を約6.8×10⁻¹³ M☉から約10⁻¹² M☉の範囲で変化させたシミュレーションを実施し、不安定性の臨界質量閾値を同定する。

実験結果

リサーチクエスチョン

• RQ1 中性子星の降着マウンドの磁気的平衡状態が、どのマウンド質量でMHD不安定性によって最初に破壊されるか。
• RQ2 降着物質によって引き起こされる局所的磁場歪みが、高質量X線連星系におけるサイクロトロン共鳴散乱特徴（CRSF）の形状および変動性にどのように影響するか。
• RQ3 MHD不安定性が、中性子星における長期間にわたる磁場埋没をダイアマグネティックスクリーニングによってどの程度制限するか。
• RQ4 圧力駆動型MHD不安定性の成長 timescale は、マウンド質量および磁場曲率にどのように依存するか。
• RQ5 不安定性閾値におけるプラズマベータ（β）は、理論的予測（例：β ≈ 260）とどの程度一致するか。

主な発見

• MHD不安定性は、マウンド質量が約5×10⁻¹³ M☉を超えると、降着マウンドの平衡状態を破壊する。これは、約45 mのマウンド高さに対応する。
• 安定な45 mマウンドにおける最大プラズマベータは約293であり、これはLitwinら（2001）が予測した理論的不安定性閾値β ≈ 260に近く、一致している。
• より大きなマウンド（約10⁻¹² M☉）では、MHD不安定性が急速に成長し、数アレクサンドラ時間（約2.8×10⁻³ s）以内に径方向端縁に指状チャネルが形成される。
• より小さなマウンド（例：50 m、約6.8×10⁻¹³ M☉）では不安定性の成長が遅く、大きなマウンドに比べて約10倍の時間スケールを示す。
• 顕著な磁場歪みは、マウンド上空約1 kmまで拡大し、500 mの高さでも磁場強度がデュポール場から10%以上ずれる。
• MHD不安定性は、大規模なダイアマグネティックスクリーニング電流の形成を著しく制限し、中性子星磁場の長期的埋没を制約する。