On the Origin of the Ly$\alpha$ Damping Wing in Galaxies at $8\le z \le 10$: Explorations using the NINJA Simulations

기술 요약: $8 \le z \le 10$ 영역 은하의 Ly$\alpha$ 댐핑 윙 기원에 대하여

문제 제기
재이온화 시대(EoR)는 성간 매질(IGM)이 중성 상태에서 이온화 상태로 전이되는 시기로 특징지어진다. 제임스 웹 우주 망원경(JWST)을 통해 고적색편이($7 \le z \le 14$) 은하 스펙트럼에서 관측되는 Ly$\alpha$ 댐핑 윙은 이 시기를 탐사하는 강력한 도구를 제공하지만, 이러한 특징을 해석하는 것은 복잡하다. 댐핑 윙의 형태는 소스 은하 자체의 성간 매질(ISM), 은하 주위의 은하 주변 매질(CGM), 그리고 확산된 IGM에 의한 흡성이 결합되어 나타나는 퇴행성(degeneracy)을 가진다. 현재의 관측 결과는 높은 중성 수소 칼럼 밀도($N_{\text{H I}}$)를 가진 댐프드 Ly$\alpha$(DLA) 시스템의 높은 발생 빈도를 보여주지만, 이러한 특징이 전역적인 IGM 중성 분율($x_{\text{HI}}$)에 의한 것인지 아니면 은하 주변의 국부적인 가스 저장고에 의한 것인지는 여전히 불분명하다. 기존 모델들은 여러 규모에서의 가스 이온화 상태에 대한 물리적으로 동기화된 제약 조건 없이, 관측된 강한 흡수체의 발생 빈도와 중간 칼럼 밀도를 동시에 재현하는 데 어려움을 겪고 있다.

방법론
저자들은 별 형성, 항성 피드백, 블랙홀 성장을 포함하지만 자기 일관적인 복사 전달(RT)은 포함하지 않는 대규모 우주론적 유체 역학 시뮬레이션인 Ninja (NISER) 수트(suite)를 활용한다. EoR 환경을 모델링하기 위해, 저자들은 이 시뮬레이션과 이온화된 거품(bubble)의 토이 모델(toy-model)을 결합한 이상화된 후처리 접근 방식을 채택하였다.

본 연구는 서로 다른 가스 성분의 기여를 분리하기 위해 세 가지 뚜렷한 시나리오를 사용하여 $z=8$ 및 $z=10$에서의 모의 Ly$\alpha$ 흡수 스펙트럼을 생성한다:

1. Case I: 흡성이 균일하게 이온화되었거나(또는 부분적으로 이온화된) 확산된 IGM의 전역 중성 분율 $x_{\text{HI}}$로부터만 발생한다고 가정한다.
2. Case II: 은하들이 반지름 $R_{\text{H II}}$(0에서 400 pkpc 범위)를 갖는 구형의 완전히 이온화된 H II 영역 내에 매몰되어 있다고 가정한다. 흡성은 IGM과 거품 외부의 할로 주변 과밀 가스(CGM)의 기여를 포함하지만, 거품 내부에는 중성 가스가 없다고 가정한다.
3. Case III: Case II를 수정하여, 할로의 비리얼 반경($R_{\text{vir}}$) 내에 특정 중성 분율($x^{\text{V}}_{\text{HI}}$)을 가진 잔류 중성 가스를 포함시킨다. 이때 $R_{\text{vir}}$와 $R_{\text{H II}}$ 사이의 가스는 완전히 이온화된 상태를 유지한다.

시뮬레이션은 수치적 수렴성과 할로 질량 의존성을 평가하기 위해 세 가지 다양한 해상도 및 부피의 박스($L50N640$, $L50N1008$, $L140N1008$)를 사용한다. 모의 스펙트럼은 JWST/NIRSpec PRISM 모드의 낮은 분해능($R \sim 100-300$)에 맞추기 위해 가우시안 프로파일로 컨볼루션되며, 포아송 노이즈(SNR=30)가 포함된다. 흡성의 강도는 $N_{\text{H I}}$에 대응하는 등가 중성 수소 칼럼 밀도로 매핑되는 파라미터 $x_p$(10% 투과율에서의 속도 분리)로 정량화된다.

핵심 기여 및 결과
본 논문은 IGM, CGM, ISM의 기여가 Ly$\alpha$ 댐핑 윙을 어떻게 형성하는지에 대한 체계적인 정량화를 제시한다:

• 가스 과밀도 및 할로 질량: 시뮬레이션 결과, 은하들은 $\sim 60-100$ pkpc까지 확장되는 과밀 가스에 둘러싸여 있음이 드러났다. 이 과밀 영역($R_{\text{eq}}$)의 범위는 할로 질량에 따라 스케일링되며(저질량 할로의 $\sim 66$ pkp에서 고질량 할로의 $\sim 105$ pkpc로 증가), $z=8$과 $z=10$ 사이에는 큰 변화를 보이지 않았다.
• IGM 기여 (Case I): 균일하게 이온화된 IGM 또는 낮은 중성 분율($x_{\text{HI}} = 0.1$)을 가진 모델은 관측된 강한 DLA 시스템의 발생 빈도를 재현하는 데 실패한다. $x_{\text{HI}} = 1.0$인 완전 중성 IGM 환경에서도, 시뮬레이션된 IGM 시선(sightline)은 드물게 $N_{\text{H I}} > 10^{22} \text{ cm}^{-2}$를 초과한다. 이는 IGM 단독으로는 JWST 데이터에서 관측되는 가장 강한 흡수체들을 설명할 수 없음을 시사한다.
• H II 영역의 역할 (Case II): 이온화된 거품의 존재는 흡성 통계를 크게 변화시킨다.
  • 만약 $R_{\text{H II}} \approx 0$ (거품 없음)이라면, 모델은 $90\%$ 이상의 은하가 $N_{\text{H I}} > 10^{21} \text{ cm}^{-2}$를 가질 것이라고 예측하는데, 이는 관측 제약 조건을 초과한다. 이는 상당수의 은하 주변에 유한한 크기의 이온화 영역이 반드시 존재해야 함을 의미한다.
  • 반대로, $R_{\text{H II}} > 50$ pkpc라면, 모델은 중성 IGM 환경에서도 가장 강한 흡수체($N_{\text{H I}} > 10^{22} \text{ cm}^{-2}$)를 재현하는 데 실패한다. 이는 큰 이온화 거품이 강한 DLA 관측치를 맞추기에는 댐핑 윙을 너무 많이 억제함을 나타낸다.
• 잔류 중성 가스 (Case III): 비리얼 반경 내에 잔류 중성 가스($x^{\text{V}}_{\text{HI}} > 0$)를 포함하는 것이 가장 강한 댐프드 Ly$\alpha$ 흡수체를 생성하는 데 필수적이다. 그러나 비리얼 반경 내의 가스가 완전히 중성이라 하더라도, $N_{\text{H I}} > 10^{22} \text{ cm}^{-2}$ 시스템의 발생 빈도는 시뮬레이션에서 분해된 할로 질량에 대한 관측 추정치보다 낮게 나타난다.
• 질량 의존성: 추론된 $N_{\text{H I}}$는 항성 질량 및 할로 질량 모두와 양의 상관관계를 갖는다. 더 강한 흡수는 더 큰 가스 저장고와 밀집된 환경을 가진 더 질량이 큰 할로에서 우선적으로 발견된다.

의의 및 결론
본 논문은 $z \gtrsim 8$ 은하에서 관측되는 Ly$\alpha$ 댐핑 윙이 IGM, CGM, ISM에서의 흡성이 결합된 복합적인 특징이라고 결론짓는다. 저자들은 다음과 같이 주장한다:

1. IGM의 한계: 확산된 IGM의 중성 분율 단독으로는 JWST에서 관측되는 강한 댐프드 Ly$\alpha$ 흡수체의 높은 발생 빈도를 설명할 수 없다.
2. 국부 가스의 필요성: 가장 강한 흡수체를 재현하기 위해서는 은하의 비리얼 반경 내부와 그 즉각적인 과밀 환경으로부터의 기여가 필요하다.
3. 누락된 물리 현상: (비리얼 반경 내 중성 가스를 고려했음에도 불구하고) 강한 흡수체의 발생 빈도가 관측치와 일치하지 않는 점은, 관측된 은하들이 현재 시뮬레이션에서 분해된 것보다 더 질량이 큰 할로에 거주하고 있거나, 혹은 해결되지 않은 추가적인 중성 가스(예: 항성 탄생 구름 내 가스)가 흡성에 크게 기여하고 있음을 시사한다.
4. 방법론적 요구 사항: IGM 중성 분율($x_{\text{HI}}$)에 대한 견고한 제약을 얻기 위해서는 ISM, CGM, IGM의 이온화 상태를 동시에 고려하는 통합 모델이 필요하다. 본 연구는 저해상도 JWST 스펙트럼을 완전히 해석하기 위해 소규모 복사 전달 계산과 대규모 재이온화 시뮬레이션을 연결해야 할 필요성을 강조한다.

저자들은 자신들의 결과가 다양한 시뮬레이션 해상도와 부피에 대해 수치적으로 견고함을 강조하며, 이는 다중 규모 복사 전달을 통합하여 EoR 관측의 퇴행성을 해결하기 위한 향후 연구의 기초를 제공한다.