재이온화가 만든 킬로파섹 규모 난류와 우주 자기장 성장

초록

본 연구는 100 pc 수준의 방사선‑수력학 시뮬레이션을 통해 재이온화 전후의 IGM에서 압력 불균형이 작은 규모의 난류를 유발한다는 사실을 최초로 입증한다. 해상도가 0.5 kpc 이하일 때 Kolmogorov k⁻⁵ᐟ³ 스펙트럼이 나타나며, k≈1 kpc⁻¹에서 난류의 회전시간은 1 Gyr 미만이다. 이러한 난류는 터뷸런트 다이너모를 통해 자기장을 증폭시켜, TeV 블레이저 관측이 요구하는 B > 3×10⁻¹⁶ G(코히런스 길이 ≈ 10 kpc) 수준을 충분히 설명한다. 사전 X‑레이 가열이 클수록 난류 발생이 억제되며, 전반적인 난류 특성은 사전 가열 온도에만 의존한다.

상세 분석

이 논문은 기존에 재이온화 과정이 IGM의 열·압력 균형을 회복하는 동안 발생할 수 있는 미세 난류를 전혀 고려하지 않았던 점을 지적한다. 저자들은 RadHydro 기반의 방사선‑수력학 코드를 수정해, 250 h⁻¹ kpc 크기의 박스에 1024³ 셀(Δx≈244 h⁻¹ pc)까지 해상도를 높였으며, 2048³ 셀(Δx≈122 h⁻¹ pc)까지 확장해 수렴성을 검증했다. 이때 I‑front가 통과한 후 약 100 Myr 이내에 과밀 필라멘트와 미니‑halo가 팽창·압축되면서 온도 차이가 10⁴–5×10⁴ K 수준으로 커지고, 압력 구배가 형성된다. 압력 구배는 라일리‑테일러 불안정과 함께 전단을 일으켜, k≈1 kpc⁻¹ 규모에서 eddy turnover time이 <1 Gyr인 Kolmogorov 난류를 유도한다. 해상도가 0.5 kpc 이상이면 수치 점성에 의해 난류가 억제됨을 확인했으며, 이는 전리된 가스의 Jeans 길이와는 별개의 제한이다.

또한 X‑ray 사전 가열을 모사해 Tmin=100 K와 1000 K 두 경우를 비교했는데, Tmin=100 K에서는 여전히 대규모 난류가 관찰되지만, Tmin=1000 K에서는 거의 난류가 사라진다. 이는 사전 가열이 IGM의 미세 구조를 매끄럽게 만들어, 재이온화 후 압력 평탄화 과정에서 서로 겹치는 구조가 감소함을 의미한다.

난류가 자기장을 증폭시키는 메커니즘은 비헬리컬 터뷸런트 다이너모로, 에너지 스펙트럼이 k⁻⁵ᐟ³까지 유지되는 경우, 난류 속도 Vdr≈20 km s⁻¹와 재이온화 시점 z≈7을 가정하면 현재 B≈1.5×10⁻⁹ G·f_B^{1/2}가 생성된다. 여기서 f_B는 난류 동에너지의 자기장 전이 효율이며, 시뮬레이션 결과는 f_B≈1에 가까운 값을 시사한다.

이러한 자기장은 우주 팽창에 따라 (1+z)⁻²로 감쇠하지만, 코히런스 길이는 (1+z) 배로 늘어나며, 최종적으로 λ≈10 kpc 규모의 볼륨 충전 자기장을 형성한다. 이는 TeV 블레이저의 γ‑ray halo 부재가 요구하는 B > 3×10⁻¹⁶ G(λ≈10 kpc) 제한을 충분히 만족한다. 따라서 재이온화‑구동 난류는 은하·은하핵 주변이 아닌 전체 IGM에 걸쳐 자기장을 증폭시킬 수 있는 보편적인 메커니즘으로 제시된다.