다크 물질이 만든 오르트 구름 폭격: 원시 블랙홀과 혜성 발생

초록

이 논문은 달 질량(≈10⁻⁷ M☉)의 원시 블랙홀(PBH)이나 유사한 거시적 다크 물질 입자가 오르트 구름을 교란시켜 내부 태양계로 혜성을 보낼 수 있음을 시뮬레이션으로 보여준다. PBH가 전체 다크 물질의 10 % 정도를 차지하면 관측되는 혜성 유입률과 일치한다는 결과를 제시한다.

상세 분석

본 연구는 다크 물질이 거시적 입자 형태(주로 달 질량의 PBH 혹은 자유 부유 소행성, FFM)로 존재한다는 가정 하에, 이러한 입자들이 태양계 외곽의 오르트 구름을 통과하면서 원시 혜성(프로코멧)을 탈출시켜 내부 태양계로 유입시키는 메커니즘을 정량화한다. 저자는 먼저 오르트 구름의 질량 분포를 Dehnen(1993)의 반경 의존 밀도 n∝a r²/(r+a)² (a=10⁵ AU)와 기존 문헌의 r⁻³·⁵ 법칙 등 여러 모델과 비교해 제시하고, 전체 혜성 수 N≈4×10¹³(직경 > 10 km)이라고 가정한다. 시뮬레이션에서는 2.5×10⁵개의 프로코멧을 무작위로 배치하고, 달 질량 DM 입자가 최대 충돌 반경 b_max≈M AU(여기서 M은 입자 질량을 태양 질량 단위로 표현) 내에서 균일하게 접근한다고 가정한다. 충돌 후 전달되는 운동량 Δp=2GMm/(b v) 를 이용해 각 혜성의 궤도를 적분했으며, 시간 간격은 0.1–0.3 yr, 총 시뮬레이션 시간은 수십만 년에 이른다.

핵심 계산식은 두 단계로 나뉜다. 첫 번째는 오르트 구름을 통과하는 DM 입자의 유입률 rate₁=4πρa²v/m이며, 여기서 ρ=0.01 M☉ pc⁻³, v=220 km s⁻¹, a는 구름 반경이다. 두 번째는 시뮬레이션에서 실제로 혜성을 내부(300 AU)로 끌어들인 횟수를 이용해 보정된 혜성 유입률 rate₂를 구한다. 보정 인자는 (1) b_max보다 큰 충돌 파라미터를 무시한 면적 손실, (2) 시뮬레이션에 포함된 혜성 수와 전체 혜성 수의 비율이다. 최종적으로 rate₂는 DM 질량 10⁻⁹–10⁻⁵ M☉ 구간에서 b_max를 적절히 조정하면 연간 0.1–10개의 혜성이 내부 태양계에 도달한다는 결과를 얻는다. 특히 f_PBH=0.1(전체 DM 중 10 %가 PBH)일 경우 관측된 혜성 유입률(≈1–2 yr⁻¹)과 일치한다.

논문은 또한 기존의 별·은하 조석 이론과 비교해 DM 입자에 의한 충돌률이 별 충돌률보다 10⁷배 높다고 주장한다. 이는 DM 입자의 질량이 작아도 수밀도가 높아 단위 부피당 입자 수가 많기 때문이다. 그러나 저자는 다음과 같은 제한점을 인정한다. (1) 혜성 간 상호작용과 행성의 중력적 스캐터링을 무시했으며, (2) b_max를 임의로 설정했는데 실제 우주에서는 더 큰 충돌 파라미터도 비무시할 수 없는 영향을 미칠 수 있다. (3) 시뮬레이션에 사용된 혜성 수가 전체 대비 10⁻⁹ 수준이라 통계적 불확실성이 크다. (4) DM 입자의 속도 분포를 단순히 220 km s⁻¹와 110 km s⁻¹의 이방성 분산으로만 모델링했으며, 실제 은하역학적 복잡성을 반영하지 못한다.

마지막으로 저자는 루빈 LSST와 로만 우주망원경을 이용한 미세중력 렌즈링 탐색이 달 질량 PBH를 검출하거나 상한을 제시하는 데 유용할 것이라고 제안한다. 현재 감도는 ≈10⁻¹¹ M☉ 수준이며, 초자외선 파장에서 관측하면 한 단계 더 낮은 질량까지 탐색 가능하다고 주장한다. 향후 연구에서는 행성의 중력 효과, 더 큰 b_max, 그리고 혜성 질량 분포를 포함한 정교한 N‑body 시뮬레이션이 필요하다고 강조한다.