은하계 고립 구형 중성자 별의 공간·속도 분포와 시뮬레이션 결과

초록

본 논문은 은하계 내 고립된 오래된 중성자 별(NS)의 현재 위치와 속도를 추정하기 위해 최신 출생 시 공간·속도 분포와 형성률을 반영한 몬테카를로 궤도 시뮬레이션을 수행한다. 디스크와 팽대부 두 개의 출생 집단을 가정하고, 이중 피크와 단일 피크 두 종류의 초기 속도 분포를 각각 3 × 10⁶개의 NS에 적용한다. 동역학적 가열 효과는 무시할 수 있음을 확인했으며, 전체 NS 수를 10⁹개로 가정할 때 태양 근처의 NS 밀도는 2–4 × 10⁻⁴ pc⁻³, 규모 높이는 0.3–0.6 kpc 정도임을 제시한다. 시뮬레이션 결과는 라디오 트랜시언트와의 연관성을 검증하는 데 활용될 수 있다.

상세 분석

이 연구는 은하계 고립 구형 중성자 별(NS)의 현재 3차원 분포와 운동학적 특성을 정량화하기 위해 대규모 몬테카를로 궤도 시뮬레이션을 수행한 점이 가장 큰 강점이다. 저자는 두 개의 출생 집단을 설정했는데, 하나는 지난 12 Gyr 동안 은하 디스크에서 일정한 형성률로 생성된 NS이며, 다른 하나는 12 Gyr 전 팽대부에서 한 번에 형성된 집단이다. 디스크와 팽대부 각각이 전체 NS 인구의 40 %와 60 %를 차지한다는 가정은 관측된 별 형성 역사와 질량 분포를 반영한 합리적인 선택이다. 초기 속도 분포에 대해서는 관측된 펄사 데이터에 기반한 이중 피크(bimodal)와 단일 피크(unimodal) 모델을 각각 적용했으며, 각각 3 × 10⁶개의 가상 NS를 추적함으로써 통계적 불확실성을 최소화하였다.

시뮬레이션은 은하 중력 포텐셜(디스크, 팽대부, 암흑 물질 halo)를 포함하고, 동역학적 가열(spiral arm, molecular cloud 충돌 등)의 효과를 평가했지만, 결과적으로 가열이 NS의 장기 궤도에 미치는 영향은 미미하다는 결론을 도출했다. 이는 NS가 출생 직후 수백 km s⁻¹ 이상의 고속을 갖고 있어, 은하 구조의 작은 퍼텐셜 변동에 비해 운동 에너지가 크게 우세하기 때문이다.

시뮬레이션 결과는 현재 은하 내 NS의 공간 밀도와 규모 높이를 구체적으로 제시한다. 전체 NS 수를 10⁹개로 가정했을 때, 태양 근처(≈8 kpc 반경)에서의 NS 밀도는 2–4 × 10⁻⁴ pc⁻³이며, 이는 이전 추정치보다 약간 높은 편이다. 또한, NS의 수직 분포는 지수형으로 근사되며, 규모 높이는 초기 속도 분포에 따라 0.3 kpc(단일 피크)에서 0.6 kpc(이중 피크) 사이로 변한다. 이러한 결과는 은하계 내 고속 천체 탐색, 특히 라디오 트랜시언트(예: Fast Radio Bursts, Rotating Radio Transients)와의 연관성을 검증하는 데 유용한 기준을 제공한다.

제공된 카탈로그는 각 NS의 현재 위치, 속도 벡터, 그리고 태양 기준 적경·적위·거리·시분할 운동량 등을 포함하고 있어, 관측자들이 특정 방향이나 거리 구간에서 기대되는 NS 수를 직접 계산할 수 있게 한다. 이는 향후 대형 전파망원경(예: SKA)이나 X‑ray 관측기기의 목표 선정에 실질적인 도움을 줄 것으로 기대된다.