Accretion of Primordial Black Holes in Stellar Interiors

기술 요약: 항성 내부에서의 원시 블랙홀의 강착 (Accretion of Primordial Black Holes in Stellar Interiors)

문제 제기
원시 블랙홀(PBH)은 특히 제약이 적은 아달(sub-lunar) 질량 범위($10^{-16} M_\odot \lesssim M_{BH} \lesssim 10^{-10} M_\odot$)에서 유효한 암흑 물질 후보로 남아 있다. 별은 PBH를 탐사하기 위한 독특한 실험실을 제공하는데, 이는 포획된 PBH가 항성 핵에 박혀 있을 경우 주변 가스를 강착하여 항성의 진화를 변화시킬 수 있기 때문이다. 이전 연구들(예: Bellinger et al. 2023)은 1차원 항성 진화 코드를 사용하여 이 시나리오를 모델링했으나, 얇은 강착 원반 이론(thin accretion disk theory)에서 유도된 가정에 의존했다. 구체적으로, 그들은 고정된 복사 효율 $\eta \approx 0.08$을 가정하고 모든 질량에 대해 에딩턴 광도 한계를 적용했다. 그러나 이 질량 범위의 PBH의 경우, 본디 반경($r_B$)이 항성의 원형화 반경(circularization radius)보다 수십 배 작기 때문에 원반 형성이 불가능하다. 대신 강착은 구형(spherical)으로 일어나며, 가스는 슈바르츠칠트 반경($r_S$) 근처에서 상대론적 온도($T \sim 10^{10}–10^{12}$ K)까지 압축되고 가열된다. 이 영역에서는 미시물리(상대론적 제동 복사, 쌍생성 과정, 클라인-니시나 산란)와 흐름의 광학적 두께가 얇은 원반 패러다임과 근본적으로 다르므로, 강착률($\dot{M}$), 복사 효율($\eta$), 그리고 흐름 구조에 대한 자기 일관적인(self-consistent) 계산이 필수적이다.

방법론
저자들은 태양형 항성의 핵에 박힌 PBH로의 강착 흐름을 모델링하기 위해, 암시적 냉각 소스 항(implicit cooling source term)을 포함한 시간 의존적 구형 오일러 방정식을 해결한다.

• 미시물리: 냉각 항 $\epsilon(\rho, T)$는 상대론적 제동 복사(전자-이온 및 전자-전자)와 전자-양전자 쌍 소멸을 포함한다. 불투명도는 톰슨 산란, 자유-자유 흡수, 그리고 높은 광자 에너지에서의 클라인-니시나 억제를 고려한다.
• 수치적 접근법: 저자들은 초기 단열 본디 프로파일로부터 자기 일관적인 정상 상태(steady state)로 흐름을 진화시키기 위해 오픈 소스 파이썬 패키지인 radbondi를 사용한다. 이러한 시간 의례적 접근 방식은 필수적인데, 냉각이 지배적인 영역에서는 음속점(sonic point)이 안장점(saddle)이 아닌 초점(focus, 복소 고윳값)이 되어 표준적인 정상 상태 ODE 슈팅 방법(shooting methods)을 사용할 수 없기 때문이다.
• 영역 처리:
  • 광학적으로 얇은 영역 ($M_{BH} \lesssim 5 \times 10^{-13} M_\odot$): 오일러 솔버가 물리적인 복사 효율과 강착률을 직접 계산한다.
  • 광학적으로 두꺼운 영역 ($M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$): 본디 구체가 광학적으로 두꺼워져 광자 트래핑(photon-trapping) 영역에 진입한다. 저자들은 Begelman (1979)의 해석적 처방을 채택하여, 트래핑 반경($r_{trap}$) 내부에서 국부적으로 생성된 광자는 블랙홀로 다시 내반(advected)되는 반면, 탈출하는 광도는 $L_{esc} \lesssim 0.6 L_{Edd}$로 스스로 조절된다고 가정한다.
• 피드백: 본 연구는 본디 구체 외부로 전달되는 복사 피드백을 자기 일관적으로 모델링한다. 여기서 복사로 인해 가열된 에너지는 대류 수송(혼합 길이 이론, Mixing Length Theory)에 의해 완화되어, 순수 복사 확산이 예측하는 극단적인 온도 상승을 방지한다.

주요 결과
본 연구는 세 가지 뚜렷한 구형 강착 영역을 식별하고 다음 결과를 도출한다:

1. 세 가지 강착 영역:

   • Hot Bondi ($M_{BH} \lesssim 10^{-14} M_\odot$): 냉각이 역학적으로 무시할 만하다. 흐름은 거의 단열적이며, $r_S$ 근처에서 $T \sim 10^{11}$ K에 도달한다. 내부 흐름은 쌍생성 온도에 도달하기 전에 충돌이 없는(collisionless) 상태가 되어 쌍 평형을 억제한다. 제동 복사가 냉각을 지배한다.
   • 제동 복사 냉각 ($10^{-14} \lesssim M_{BH} \lesssim 5 \times 10^{-13} M_\odot$): 냉각이 효율적이 되어 가스가 쌍생성 온도에 도달하기 전에 열에너지를 배출한다. 흐름은 등온 상태($T \approx T_\infty$)에 가까워지며, 강착률은 단열 본디 rate에 비해 약 2~7배 증가한다.
   • 광자 트래핑 ($M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$): 본디 구체가 광학적으로 두껍다. 국부적으로 방출되는 광도는 질량에 따라 급격히 증가하지만, 탈출하는 광도는 Begelman 한계($L_{esc} \lesssim 0.6 L_{Edd}$)에 의해 제한된다. 결정적으로, 구형 에딩턴 한계가 유입되는 질량 흐름에는 적용되지 않으므로 강착률은 본디 값에 가깝게 유지된다 ($\dot{M} \propto M_{BH}^2$).

2. 복사 효율: 자기 일관적으로 계산된 복사 효율은 광학적으로 얇은 범위에서 $\eta \approx 10^{-2}$이다. 이는 기존 문헌에서 가정된 표준 얇은 원반 값인 $\eta \approx 0.08$보다 한 자릿수 낮은 수치이다. 효율은 $10^{-15} M_\odot$ 근처에서 정점을 찍은 후 냉각 및 등온 한계에서 감소한다.

3. 성장 궤적: 구형 에딩턴 한계가 적용되지 않기 때문에, PBH의 성장은 광자 트래핑 영역 전체에서 초지수적(super-exponential, $\dot{M} \propto M_{BH}^2$) 성장을 유지한다. 성장 병목 현상은 가장 낮은 질량에서의 느린 Hot Bondi 단계에 의해 결정된다.

4. 임계 초기 질량: PBH가 허블 시간 내에 태양 질량의 별을 소멸시키는 데 필요한 임계 초기 질량은 $M_{0,crit} \sim 10^{-16} M_\odot$이다. 이는 더 높은 효율과 에딩턴 병목 현상을 가정했던 이전 추정치($\sim 10^{-11} M_\odot$)보다 다섯 자릿수 낮은 수치이다.

5. 자기장: 저자들은 자기 선속(magnetic flux) 축적이 미미함을 발견했다. 태양 핵의 주변 자기장은 흐름을 역학적으로 저지하기에는 너무 약하며, 자기 재결합(magnetic reconnection)은 선속 축적 수준을 역학적으로 중요하지 않은 수준($P_{mag}/P_{ram} \ll 1$)으로 제한한다.

의의 및 주장
본 논문은 아달 질량 범위 전반에 걸쳐 구형 PBH 강착에 대한 복사 효율과 강착률을 최초로 자기 일관적으로 계산했다고 주장한다. 얇은 원반 가정을 대신하여 구형 흐름의 미시물리적 처리를 도입함으로써, 저자들은 다음을 입증한다:

• 복사 효율이 기존에 가정되었던 것보다 현저히 낮다.
• 구형 강착에는 에딩턴 제한 단계가 존재하지 않으며, 흐름은 복사압에 의해 억제되지 않고 질량 흐름 측면에서 초에딩턴(super-Eddington) 상태를 유지한다.
• 결과적으로, 항성 관측으로부터 얻은 PBH 암흑 물질에 대한 제약은 이전보다 훨씬 강력하다. 왜냐하면 초기 질량이 $10^{-16} M_\odot$만큼 낮은 PBH라도 허블 시간 내에 모항성을 집어삼킬 만큼 성장할 수 있기 때문이다.

저자들은 본 연구가 강착 미시물리에 초점을 맞추고 있으며, 항성의 전역적 구조에 대한 역작용(back-reaction)은 포함하지 않았음을 언급한다(이는 MESA와 같은 항성 진화 코드와의 결합이 필요함). 그러나 저자들은 향후 연구를 가능하게 하는 데 필요한 $\eta(M_{BH})$ 및 $\dot{M}(M_{BH})$ 관계식을 제공한다. 또한, 원반 형성(더 높은 질량 또는 더 빠르게 회전하는 별에서 발생)에 대한 논의는 동료 논문으로 미룬다.