Mass distribution of neutron stars in binary systems

技術的サマリー：連星系における中性子星の質量分布

問題提起
観測結果は、連星系に存在する既知の中性子星（NS）の質量分布が、約 1.4 M$_\odot$ と $\sim$1.8 M$_\odot$ 付近にピークを持つ二峰性のパターンを示すことを示している。1.4 M$_\odot$ 付近のピークは、通常のコア崩壊型超新星や共通エンベロープ進化に起因することが多いが、$\sim$1.8 M$_\odot$ 付近のピークの起源については依然として議論が続いている。中性子星の質量成長を降着によって説明しようとする以前の理論モデルは、エディントン限界によって制限されていたか、あるいは超エディントン降着領域を十分に考慮していなかった。さらに、質量移動の効率、プロペラー効果、磁場進化に関する不確実性が、降着する中性子星の最終的な質量分布の正確な予測を妨げてきた。本研究は、超エディントン降着モデルが、連星系における中性子星の観測された二峰性質量分布を再現できるかどうかを調査することを目的としている。

手法
著者は、2 つの異なるコードを組み合わせたハイブリッドな集団合成アプローチを採用した。

1. 連星進化（BSE）： $10^7$ 個の原始連星系の形成と初期進化をシミュレートするために使用された。このコードは、初期質量関数、質量比、離心率、および出生時のキック速度を取り入れ、主系列星の伴星を持つ中性子星連星の集団を生成した。系は、安定した質量移動が可能なものを選別し、不安定な共通エンベロープ進化に陥るものや極端な質量移動率（$>10^4 \dot{M}_{Edd}$）を持つものを除外した。
2. 恒星天体物理学実験用モジュール（MESA）： BSE の出力から選ばれた中性子星連星は、MESA（バージョン r24.08.1）を用いた詳細な進化計算に付された。この段階は、安定したローシュ限界溢出条件における降着過程に焦点を当てた。

主要な物理モデルと仮定：

• 超エディントン降着： エディントン限界によって制限されていた以前のモデルとは異なり、本研究は Lipunov (1987) が記述する超エディントン降着の物理を取り入れている。このモデルは、アルヴェーン半径（$R_a$）と球状化半径（$R_s$）の関係に基づき、臨界未満の領域と臨界超過の領域を区別する。
• 降着効率： 中性子星によって降着される移動質量の割合（$f_{mt,a}$）は、降着状態によって決定される。臨界超過の領域では、過剰な物質は中性子星近傍の放射圧によって吹き飛ばされ、特定の軌道角運動量を運び去る。保持される質量については、降着効率を 0.5 と仮定する。
• 磁場進化： 中性子星の磁場は、降着質量（$\Delta M$）の関数として、指数法則に従って減衰し、初期の $5 \times 10^{12}$ G から最小値 $5 \times 10^8$ G まで減少する。
• スピン進化： このモデルは、中性子星のスピン周期を 3 つの段階（放射：降着なし、プロペラー：物質の放出、降着：スピンアップ）で追跡する。これらの段階間の遷移は、光円筒、同期半径、および磁気圏半径の相対的な位置に依存する。
• 角運動量の損失： シミュレーションは、重力波放射、磁気ブレーキング、および放出された質量と連星周囲円盤の形成による角運動量の損失を考慮する。

主要な結果

1. 質量成長メカニズム： 本研究は、中性子星の質量成長が連星の軌道周期と質量供給星の質量に大きく依存することを示している。超エディントンモデルにおいて、中性子星は降着平衡状態に長期間留まることができ、これにより顕著な質量蓄積（1.4 M$_\odot$ の出生質量から $\sim$1.0 M$_\odot$ までの成長）が可能となる。これは、中性子星がしばしば質量成長を制限する長期のプロペラー段階に入る Li ら (2021) のモデル（エディントン限界）とは対照的である。
2. 二峰性の再現： 集団合成の結果は、観測された二峰性分布を成功裡に再現した。
   • $\sim$1.8 M$_\odot$ のピーク： このピークは、主に質量供給星の質量が $M_d < \sim$1.6 M$_\odot$ で軌道周期が $P_{orb} < 20$ 日である中性子星連星から生じる。これらの系では、質量移動率が長期間にわたって効率的な準平衡の亜エディントン降着を可能にする。
   • 1.4 M$_\odot$ のピーク： このピークは、共通エンベロープ進化を経る中性子星系（しばしばより質量の大きな供給星や特定の軌道構成を伴う）に起因し、質量成長は最小限である。
   • 高質量テール： 中間質量の供給星（$1.6 \le M_d \le 3.0$ M$_\odot$）と短い軌道周期を持つ系は、2.1 M$_\odot$ を超えるのに十分な物質を降着させることができる。
3. スピン周期分布： このモデルは、低質量 X 線連星（LMXB）において、中性子星のスピン周期が 1 ms から 15 ms の間に集中することを予測する。結果は一部の観測データよりもわずかに大きいが、著者らはこれを、降着円盤が中性子星表面に近すぎないようにする最小磁場 $5 \times 10^8$ G の仮定に起因すると帰している。

意義と主張
本論文は、中性子星の質量分布における高質量ピークを説明するために、超エディントン降着の物理を取り込むことが不可欠であると主張している。質量移動率がエディントン限界を超えていても中性子星が効率的に降着することを可能にすることで、このモデルは、理論的に予測される低い質量成長と、観測される質量の大きな中性子星の集団との間の不一致を解決する。著者らは、自らのモデルが高質量ピークのための異様な形成チャネルを必要とせず、標準的な連星進化と超エディントン降着に起因する二峰性分布に対して一貫した物理的説明を提供すると主張している。

限界と不確実性
著者らは、自らのモデルにおけるいくつかの不確実性を認めている。

• 降着効率： 仮定された効率 0.5 と、中性子星近傍での質量損失の扱いは簡略化されている。
• 磁場減衰： 減衰率と最小磁場強度は、スピン進化と降着効率に大きく影響するパラメータである。
• 状態方程式（EoS）： 慣性モーメントは $10^{45}$ g cm$^2$ に固定されており、これは下限に近い。中性子星の状態方程式の変動は慣性モーメントを変化させ、スピンアップの時間スケールやプロペラー領域への遷移に影響を与える可能性がある。
• 除外されたチャネル： 本研究は、降着誘起崩壊（AIC）によって形成された中性子星や、二峰性質量分布を持って生まれた中性子星の寄与を定量化しておらず、代わりに標準的なコア崩壊とそれに続く連星進化によって形成された中性子星に焦点を当てている。