Can accreting isolated neutron stars be detected?

기술적 요약: 인구 합성을 통한 고립된 중성자별의 강착 감지 가능성

문제 제기
은하계 내 중성자별 (NS) 의 개수는 약 $10^9$개 수준으로 추정되지만, 관측된 고립된 중성자별은 $10^4$개 미만이며, 이는 주로 젊은 전파 펄사들이다. 노년기 고립 중성자별의 운명은 여전히 중요한 미해결 과제로 남아 있다. 수십 년간 노년 중성자별이 성간 매질 (ISM) 물질을 강착하여 X 선을 방출함으로써 관측 가능해질 수 있다는 가설이 제기되어 왔으나, ROSAT, Chandra, XMM-Newton 과 같은 관측 노력들은 그러한 소스의 확인된 집단을 식별하지 못했다. 이전의 인구 합성 연구들은 낮은 중성자별 속도, 높은 강착 효율, 그리고'프로펠러'단계를 통한 급속한 진화를 가정하여 낙관적인 추정치를 제시했다. 그러나 높은 태생 킥 (natal kick) 속도의 발견과 낮은 강착 효율을 시사하는 수치 시뮬레이션은 이러한 가정들에 도전했다. 본 논문은 은하계 내 고립 중성자별에 대한 상세한 인구 합성을 수행하여 이러한 불일치를 해결하는 것을 목표로 하며, 최신 ISM 모델, 자기 - 회전 진화, 그리고 불확실한 프로펠러 단계 물리학을 통합한다.

방법론
저자들은 은하계의 수명 ($13.6$ Gyr) 동안 고립된 중성자별의 진화를 시뮬레이션하는 인구 합성 접근법을 사용한다. 방법론은 다음과 같은 주요 구성 요소들을 통합한다:

1. 초기 매개변수와 운동학: 합성 중성자별은 전파 펄사와 자기성 (magnetar) 집단으로부터 유도된 분포에 기반한 초기 자전 주기 ($P_0$), 자기장 ($B_0$), 그리고 공간적 위치로 생성된다. 초기 속도는 원형 은하 속도, 잔류 전구체 속도, 그리고 이모달 맥스웰 분포에서 추출된 태생 킥 속도의 벡터 합으로 계산된다. 궤적은 현실적인 은하 중력 퍼텐셜 내에서 계산된다.
2. 성간 매질 (ISM) 모델: 두 가지 ISM 모델이 비교된다:
   • 단일 냉각 상 (분자, 중성, 온전 이온화 가스 및 먼지) 으로 구성된 단일 상 모델.
   • 냉각 상과 고온 상 (코로나/헤일로 가스) 을 구분하며, 중성자별 궤적을 따라 각 상의 국소 충전 계수를 결정하기 위해 압력 평형을 활용하는 이중 상 모델.
3. 자기 - 회전 진화: 연구는 중성자별을 방출기 (ejector), 프로펠러, 강착기 (accretor), 그리고 지회전체 (georotator) 의 네 가지 진화 단계를 거쳐 추적한다. 전이는 임계 반경 (광원통, 알프벤 반경, 슈바르츠만 반경, 동조 반경) 에 의해 지배된다.
   • 자전 감속: 오일러 방정식을 사용하여 자전 진화를 모델링한다. 특히, 연구는 프로펠러 단계 동안의 자전 감속 토크 ($K_P$) 에 대해 네 가지 서로 다른 모델 (A, B, C, D) 을 테스트하며, 이는 매우 효율적인 (모델 A) 에서 비효율적인 (모델 D) 에 이르는 광범위한 이론적 효율을 나타낸다.
   • 자기장: 두 가지 자기장 진화 시나리오가 테스트된다: 일정한 자기장 (CF) 과 자기성에서의 홀 캐스케이드 효과 및 장기적인 옴 소산을 고려한 지수적으로 감쇠하는 자기장 (ED).
4. 관측 가능성 계산: 강착기 단계에 도달하는 중성자별에 대해, $100\%$ 효율 ($\zeta=1$) 의 본디 - 호일 - 리틀턴 강착을 가정하여 X 선 광도를 계산한다. 스펙트럼은 극관 (polar caps) 에서의 흑체로 모델링되며, 중력 적색 편이와 시간 지연이 보정된다. 성간 흡수는 HI4PI 데이터와 일관된 column density 지도를 적용한다. eROSITA 망원경에 의한 감지 가능성은 계수율을 계산하여 평가하며, 감지 임계값은 $0.01$ counts s$^{-1}$로 설정된다.

주요 기여 및 결과

• 프로펠러 단계 불확실성의 우세: 주요 발견은 자전 감속 토크 모델의 선택에 의해 결정되는 프로펠러 단계의 지속 시간이 관측 가능한 강착 중성자별의 수에 영향을 미치는 가장 중요한 요소라는 것이다.
  • 효율적인 자전 감속 (모델 A 및 B): 프로펠러 단계가 짧다면 (급속한 자전 감속), 중성자별의 상당 부분 ($17-53\%$) 이 강착기 단계에 도달한다. 가장 낙관적인 시나리오 (모델 A, 일정한 자기장, 단일 상 ISM) 에서 연구는 약 $13,000$개의 잠재적 관측 가능 소스를 예측한다. 더 현실적인 시나리오 (모델 B, 이중 상 ISM, 감쇠하는 자기장) 에서는 이 수가 약 $3,000-4,000$개로 감소한다.
  • 비효율적인 자전 감속 (모델 C 및 D): 자전 감속이 비효율적이라면 (예: Davies 와 Pringle 1981 에 기반한 모델 D), 프로펠러 단계는 은하계의 나이보다 더 길어질 수 있다. 결과적으로 강착기 단계에 도달하는 중성자별의 수는 무시할 수 있을 정도로 작아져, eROSITA 에 의해 감지 가능한 소스는 사실상 제로가 된다.
• 속도 의존성: 중성자별의 ISM 에 대한 특성 속도는 두 번째로 중요한 매개변수이다. 강착률은 $\dot{M} \propto v^{-3}$로 비례한다. 따라서 관측 가능한 소스는 거의 예외 없이 저속 천체들이다 (모델 A/B 의 경우 $v \lesssim 60-70$ km s$^{-1}$, 모델 C 의 경우 $v \lesssim 40$ km s$^{-1}$). 고속 중성자별은 강착기 단계에 도달하거나 감지될 만큼 충분한 플럭스를 생성하는 경우가 드물다.
• 자기장 진화: 자기장 감쇠는 방출기와 프로펠러 단계에 있는 중성자별의 비율을 크게 변화시킨다 (상호작용 반경을 축소하고 전이를 용이하게 함). 그러나 이는 강착기의 총 수에 미치는 영향이 프로펠러 자전 감속 모델에 비해 부차적이다. 다만, 자기장 감쇠는 자기권과 극관 반경의 축소로 인해 관측 가능한 소스의 유효 온도를 ($\sim 5 \times 10^6$ K 에서 $\sim 2 \times 10^6$ K 로) 낮춘다.
• 공간적 분포: 잠재적 관측 가능 강착기는 상대적으로 태양에 가까운 곳 (약 $400-550$ pc 이내) 에 위치하며, 은하면 근처의 얇은 원반 ($z < 470$ pc) 에 제한될 것으로 예측된다.
• ISM 모델 영향: 단일 상에서 이중 상 ISM 모델로의 전환은 중성자별이 진화가 느린 저밀도 영역에서 더 많은 시간을 보내게 함으로써 예측된 강착자의 수를 modest 하게 감소시킨다 (2 배 미만의 감소).

의의 및 주장
본 논문은 지금까지 강착 고립 중성자별이 감지되지 않은 이유는 소스의 부재나 관측 민감도 부족이 아니라, 프로펠러 단계 물리학의 불확실성 때문일 가능성이 높다고 주장한다. 연구는 프로펠러 자전 감속이 충분히 빠르다면 eROSITA 가 수천 개의 이러한 소스 집단을 감지할 잠재력을 가지고 있음을 보여준다. 반대로, 자전 감속이 비효율적이라면 그러한 집단은 관측 가능한 우주에 사실상 존재하지 않을 수 있다.

저자들은 현재 프로펠러 영역에 대한 이론적 모델이 관측에 의해 제약받지 않고 있음을 강조한다. 그들은 Gaia 에 의해 최근 발견된 광대역 저질량 쌍성계에서의 중성자별 관측이 중요하다고 제안한다. 이러한 시스템은 중성자별이 성간 매질과 유사한 방식으로 항성풍과 상호작용하므로, 프로펠러 자전 감속 토크와 강착 효율을 제약하는 데 필요한 경험적 데이터를 제공할 수 있으며, 이를 통해 고립된 강착 중성자별의 감지 가능성을 명확히 할 수 있다. 논문은 프로펠러 단계와 저율 강착 효율의 불확실성을 해결하지 않는 한 관측 가능한 집단에 대한 정확한 예측은 여전히 요원하다고 겸손하게 결론지었다.