[논문 리뷰] Core-collapse, evaporation and tidal effects: the life story of a self-interacting dark matter subhalo

Zhichao Carton Zeng, Annika H. G. Peter외 6명|arXiv (Cornell University)|2021. 10. 01.

Cosmology and Gravitation Theories참고 문헌 180인용 수 70

한 줄 요약

이 논문은 렌즈 질량을 갖는 주위 허브 홀로에서 자기상호작용 다크 물질(SIDM)의 하위구조가 코어 형성에서 코어 붕괴에 이르는 진화를 모의하기 위해 하이브리드 반분석적 및 N-body 방법을 개발한다. 연구 결과, 일정한 교차단면을 갖는 하위구조에서는 주위-하위구조 상호작용으로 인한 증발이 코어 붕괴를 크게 지연시키거나 방지하는 것으로 나타났으며, 이는 관측된 초박공한 다크 하위구조가 일정한 교차단면 모델이 아니라 속도에 의존하는 교차단면 또는 추가 물리적 메커니즘을 필요로 함을 시사한다.

ABSTRACT

Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass halos, accelerated if halos are subhalos of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of $\lesssim 10^8 M_\odot$ subhalos in lens-mass systems ($\sim 10^{13}M_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhalos with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhalos ($\lesssim 1/1000$ host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhalos to core-collapse, finding that it is nearly impossible to drive core-collapse in subhalos in SIDM models with constant cross sections. Any discovery of ultra-compact dark substructures with future substructure lensing observations favors additional degrees of freedom, such as velocity-dependence, in the cross section.

연구 동기 및 목표

• 렌즈 질량을 갖는 주위 허브 홀로 내에서 SIDM 하위구조의 전체 진화 과정—코어 형성에서 코어 붕괴에 이르기까지—를 연구하기 위해.
• 기존의 N-body 방법을 사용할 때 저질량 하위구조(≲10⁸ M⊙)를 거대한 주위 허브 시스템(∼10¹³ M⊙)에서 시뮬레이션하는 데 발생하는 계산적 과제를 해결하기 위해.
• 내부 열의 유출, 주위-하위구조 증발, tidal 효과가 하위구조 진화에 미치는 역할을 정량화하기 위해.
• 일정한 교차단면을 갖는 SIDM 하위구조에서 코어 붕괴가 발생할 수 있는 매개변수 공간을 규명하기 위해.
• 향후 하위구조 렌즈 관측에서 SIDM 모델을 구분하는 데 미치는 영향을 평가하기 위해.

제안 방법

• 하위구조 진화의 단계별 모의를 위해 반분석적 모델링과 고정밀 N-body 시뮬레이션을 융합한 하이브리드 시뮬레이션 프레임워크를 개발한다.
• 평균 자유로운 경로가 중력 스케일 높이와 유사해지는 Knudsen 수(𝐾𝑛 ≈ 1) 기반의 코어 붕괴 기준을 적용한다.
• 고정된 연화 길이를 갖는 연화된 중력 힘 법칙을 사용하며, 코어 붕괴 근처에서 수치 오차가 최소화됨을 검증하였다.
• 다수의 페리세트릭스 통과 동안 하위구조의 진화를 추적하여 입자 증발 및 tidal 스트리핑 효과를 해석한다.
• 계산 자원을 하위구조 진화의 핵심 단계에 집중시키기 위해 단계별 시뮬레이션 접근법을 채택한다.
• 다이나믹 프랙션 오류가 최소화되는 조건(하위구조 질량 ≤ 주위 허브 홀로 질량의 1/1000 이내)에서 방법을 검증하였다.

실험 결과

연구 질문

• RQ1 렌즈 질량 주위 허브 홀로 내에서 SIDM 하위구조의 진화는 어떤 요인에 의해 코어 형성에서 코어 붕괴에 이르게 하는가?
• RQ2 증발과 tidal 효과는 SIDM 하위구조의 코어 붕괴 시기와 가능성에 어떤 영향을 미치는가?
• RQ3 실제 천문학적 농도 조건에서 일정한 교차단면을 갖는 SIDM 하위구조에서 코어 붕괴가 가능할 수 있는가?
• RQ4 관측된 초박공한 다크 하위구조는 SIDM 교차단면의 속도 의존성에 어떤 제약을 가하는가?
• RQ5 열의 유출, 증발, tidal 가열이라는 경쟁적 과정들이 하위구조의 중심 밀도 변화에 어떻게 영향을 미치는가?

주요 결과

• 주위 허브 홀로 입자와의 상호작용로 인한 하위구조 입자의 증발은 SIDM 하위구조에서 코어 붕괴를 지연시키거나 방지할 수 있는 주요 과정이다.
• 실제 천문학적 농도 조건에서 일정한 교차단면을 갖는 하위구조에서는 코어 붕괴가 거의 불가능하며, 증발과 tidal 가열의 영향이 내부 열의 유출보다 일반적으로 더 크기 때문이다.
• 하위구조의 중심 밀도는 내부 열의 유출가 증발과 tidal 가열로 인한 에너지 증가를 초월할 때에만 증가한다.
• 하위구조 [𝑀sub = 10¹⁰.⁵ M⊙, 𝜎T/𝑚 = 6 cm²/g, 𝑟peri:𝑟apo = 1/10, 𝑐 = 75]는 냉각과 가열 사이에 거의 임계적인 균형을 이루며, 이로 인해 코어 붕괴는 시기와 초기 조건에 매우 민감하게 영향을 받는다.
• 향후 하위구조 렌즈를 통해 초박공한 다크 하위구조를 관측할 경우, 일정한 교차단면을 갖는 SIDM 모델은 강하게 배제되고, 속도에 의존하는 교차단면 또는 추가적인 자유도가 유리할 것이다.
• 이 방법은 코어 붕괴 근처에서 짧은 평균 자유로운 경로 영역으로의 전이를 성공적으로 캡처하였으며, 이는 N-body 시뮬레이션을 그 시점에서 종료하는 데 타당성을 부여하고, 유체역학적 또는 분석적 후속 연구의 필요성을 시사한다.