Mass distribution of neutron stars in binary systems

기술 요약: 쌍성계 중성자별의 질량 분포

문제 제기
관측에 따르면, 쌍성계에 있는 알려진 중성자별 (NS) 들의 질량 분포는 1.4 M$_\odot$와 $\sim$1.8 M$_\odot$ 부근에 피크를 보이는 이분모 (bimodal) 패턴을 나타냅니다. 1.4 M$_\odot$ 피크는 종종 표준 핵붕괴 초신성이나 공통-envelope 진화로 귀속되지만, $\sim$1.8 M$_\odot$ 피크의 기원은 여전히 논쟁의 대상입니다. 강착을 통한 중성자별 질량 증가를 설명하려는 이전의 이론적 모델들은 종종 에딩턴 한계에 의해 제한받거나, 초에딩턴 강착 영역을 완전히 고려하지 못했습니다. furthermore, 질량 전달 효율, 프로펠러 효과, 그리고 자기장 진화에 대한 불확실성은 강착 중성자별의 최종 질량 분포에 대한 정확한 예측을 방해해 왔습니다. 본 연구는 초에딩턴 강착 모델이 쌍성계 중성자별의 관측된 이분모 질량 분포를 성공적으로 재현할 수 있는지 조사하는 것을 목표로 합니다.

방법론
저자들은 두 가지 서로 다른 코드를 결합한 하이브리드 집단 합성 (population synthesis) 접근법을 사용했습니다:

1. 쌍성 진화 (BSE): $10^7$개의 원시 쌍성계의 형성과 초기 진화를 시뮬레이션하는 데 사용되었습니다. 이 코드는 초기 질량 함수, 질량비, 이심률, 그리고 태어날 때의 속도 (natal kick velocities) 를 포함하여 주계열성 동반자를 가진 중성자별 쌍성군을 생성했습니다. 시스템은 안정된 질량 전달이 가능한 것들을 식별하기 위해 필터링되었으며, 불안정한 공통-envelope 진화나 극단적인 질량 전달률 ($>10^4 \dot{M}_{Edd}$) 을 가진 시스템은 제외되었습니다.
2. 항성 천체물리학 실험 모듈 (MESA): BSE 출력에서 선택된 중성자별 쌍성계는 MESA(버전 r24.08.1) 를 사용하여 상세한 진화 계산을 받았습니다. 이 단계는 안정된 로슈 로브 넘침 (Roche lobe overflow) 조건 하에서의 강착 과정에 초점을 맞추었습니다.

주요 물리 모델 및 가정:

• 초에딩턴 강착: 에딩턴 한계에 제한된 이전 모델들과 달리, 본 연구는 Lipunov (1987) 가 기술한 초에딩턴 강착의 물리를 포함합니다. 이 모델은 알프벤 반경 ($R_a$) 과 구형화 반경 ($R_s$) 사이의 관계에 기반하여 아임계 (sub-critical) 영역과 초임계 (super-critical) 영역을 구분합니다.
• 강착 효율: 중성자별에 강착되는 전달 질량의 비율 ($f_{mt,a}$) 은 강착 상태에 의해 결정됩니다. 초임계 영역에서는 과잉 물질이 중성자별 근처의 복사압에 의해 날아가며, 특정 궤도 각운동량을 운반합니다. 유지되는 질량에 대한 강착 효율은 0.5 로 가정됩니다.
• 자기장 진화: 중성자별의 자기장은 강착된 질량 ($\Delta M$) 의 함수로 지수 법칙에 따라 감쇠하며, 초기 $5 \times 10^{12}$ G 에서 최소 $5 \times 10^8$ G 로 감소합니다.
• 스핀 진화: 이 모델은 중성자별의 스핀 주기를 세 단계 (무강착 상태인 전파, 물질 분출 상태인 프로펠러, 스핀 업 상태인 강착) 를 추적합니다. 이러한 단계 간의 전이는 광원통 (light cylinder), 동조 반경 (corotation radius), 그리고 자기권 반경의 상대적 위치에 따라 결정됩니다.
• 각운동량 손실: 시뮬레이션은 중력파 복사, 자기 제동, 그리고 분출된 질량과 쌍성계 외곽 원반 형성으로 인한 각운동량 손실을 고려합니다.

주요 결과

1. 질량 증가 메커니즘: 본 연구는 중성자별의 질량 증가가 쌍성 궤도 주기와 donor 별의 질량에 크게 의존함을 보여줍니다. 초에딩턴 모델에서 중성자별은 확장된 기간 동안 강착 평형 상태에 머무를 수 있어, 1.4 M$_\odot$의 탄생 질량에서 최대 $\sim$1.0 M$_\odot$까지의 상당한 질량 축적을 가능하게 합니다. 이는 Li et al. (2021) 의 모델 (에딩턴 제한) 과 대조적인데, 해당 모델에서는 중성자별이 종종 질량 증가를 제한하는 장기간의 프로펠러 단계에 진입합니다.
2. 이분모 재현: 집단 합성 결과는 관측된 이분모 분포를 성공적으로 재현합니다:
   • $\sim$1.8 M$_\odot$ 피크: 이 피크는 주로 donor 질량 $M_d < \sim$1.6 M$_\odot$이고 궤도 주기 $P_{orb} < 20$일인 중성자별 쌍성계에서 발생합니다. 이러한 시스템에서 질량 전달률은 장기간에 걸쳐 효율적인 준평형 아에딩턴 강착을 가능하게 합니다.
   • 1.4 M$_\odot$ 피크: 이 피크는 공통-envelope 진화를 겪는 중성자별 시스템 (종종 더 무거운 donor 또는 특정 궤도 구성을 포함) 에 귀속되며, 여기서 질량 증가는 미미합니다.
   • 고질량 꼬리: 중간 질량 donor ($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$로 고정되어 있으며, 이는 하한선 근처입니다. 중성자별의 EoS 변화는 관성 모멘트를 변경하여 스핀 업 시간 척도와 프로펠러 영역으로의 전이에 영향을 줄 수 있습니다.
• 제외된 채널: 본 연구는 강착 유도 붕괴 (AIC) 를 통해 형성된 중성자별이나 이분모 질량 분포로 태어난 중성자별의 기여도를 정량화하지 않으며, 대신 표준 핵붕괴와 후속 쌍성 진화를 통해 형성된 중성자별에 초점을 맞춥니다.