TRINITY: A coupled model of winds, radiation, and photoionised gas in molecular clouds. I. Methods and validation

기술 요약: TRINITY: 분자 구름 내의 바람, 복사, 그리고 광이온화 가스의 결합 모델 – I. 방법론 및 검증

문제 및 맥락
거대 분자 구름 (GMC) 은 대질량 별이 탄생한 후 1~5 Myr 내에 분산하며, 보통 첫 번째 핵붕괴 초신성이 도달하기 전에 분산합니다. 관측 조사 (예: PHANGS, LEGUS) 는 초신성 이전 피드백 메커니즘인 항성풍, 복사압, 그리고 광이온화 가스 압력 ($P_{\text{H II}}$) 이 구름을 비우는 주요 동력임을 시사합니다. 그러나 이러한 채널들의 상대적 우세성, 그 시기, 그리고 구름의 특성 (질량, 밀도 프로파일, 별 형성 효율) 에 대한 의존성은 여전히 불확실합니다. 기존 반해석적 도구들은 종종 이러한 피드백 메커니즘들을 고립적으로 다루거나 단순화된 밀도 프로파일에 의존하는 반면, 3 차원 시뮬레이션은 통계적 조사를 수행하기에는 계산 비용이 너무 큽니다. 따라서 다양한 구름 구조에 대한 모든 초신성 이전 피드백 채널을 결합하여 매개변수 공간 전반에 걸친 피드백 우세성을 매핑할 수 있는 효율적이고 자기 일관적인 순방향 모델이 필요합니다.

방법론
저자들은 분자 구름 내에서 피드백에 의해 구동되는 버블과 쓸려 올라간 껍질을 진화시키는 1 차원, 구대칭 얇은 껍질 코드인 Trinity를 소개합니다. Trinity 는 warpfield 프레임워크를 계승하지만 두 가지 중요한 구조적 변화를 도입했습니다:

1. 명시적 광이온화 가스 압력: $P_{\text{H II}}$가 경계 조건으로만 처리되거나 무시되는 것이 아니라, 껍질 운동 방정식에 명시적으로 포함됩니다.
2. 위상 인지형 구동 처방: 코드는 이중 계산을 방지하기 위해 진화 위상에 따라 구동 압력 항을 동적으로 선택합니다:
   • 에너지 구동 위상: 구동 압력은 뜨거운 버블의 열압 ($P_b$) 과 광이온화 가스 압력 ($P_{\text{H II}}$) 중 최댓값입니다.
   • 전이 위상: 구동 압력은 남은 $P_b$와 $P_{\text{H II}}$에 새로 주입된 바람과 초신성 잔해의 램압 ($P_{\text{ram}}$) 을 더한 값 중 최댓값입니다.
   • 운동량 구동 위상: 뜨거운 버블이 냉각되어 ($E_b \to 0$) 운동량 구동 위상에 도달하면, 구동 압력은 $P_{\text{H II}} + P_{\text{ram}}$의 합입니다.

이 모델은 버블 반지름 ($R_b$), 속도 ($v_b$), 그리고 열에너지 ($E_b$) 에 대한 결합된 상미분 방정식 시스템을 풉니다. 이는 다음을 포함합니다:

• 초기 조건: 균일, 조각별 멱함수 ($\rho \propto r^\alpha$), 그리고 Bonnor–Ebert (BE) 구 밀도 프로파일을 지원합니다.
• 물리 모듈: 항성풍, 초신성, 직접 및 먼지 재처리 복사압, 껍질과 성단의 자체 중력, 그리고 껍질에서 뜨거운 버블로의 전도성 증발을 포함합니다.
• 복사 - 껍질 결합: 정역학적 껍질 구조 계산을 통해 이온화면 위치 ($R_{\text{if}}$), 흡수된 복사 분율, 그리고 적외선 광학 두께를 자기 일관적으로 결정합니다.
• 냉각: 충돌 이온화 평형 (CIE) 에서 비-CIE 영역으로의 전이를 처리하고 에너지 - 운동량 전이를 결정하기 위해 냉각 손실을 계산합니다.

주요 결과
본 논문은 코드를 분석적 한계 (Weaver 바람, Spitzer 팽창, 눈썰매) 와 warpfield 실행 결과에 대해 검증한 후, $M_{\text{cloud}} = 10^5–10^{6.5} M_\odot$, 중심 밀도 $n_{\text{core}} = 10^3–10^4 \text{ cm}^{-3}$, 그리고 별 형성 효율 $\epsilon = 0.01–0.30$인 구름에 대한 결과를 제시합니다.

1. 광이온화 가스의 영향: $P_{\text{H II}}$를 포함하면 역학에 측정 가능한 변화를 줍니다. 기준 검증 실행 ($M_{\text{cloud}}=10^6 M_\odot, \epsilon=0.01$) 에서 $P_{\text{H II}}$를 포함하면 $P_{\text{H II}}$가 비활성화된 warpfield와同等한 실행에 비해 10 Myr 시점에서 껍질 반지름이 약 17% 증가합니다. $P_{\text{H II}}$는 운동량 지배적 진화로 전환된 후에도 역학적으로 중요합니다.
2. 밀도 프로파일 의존성: 항성 집단이 고정되어 있더라도 구름 구조가 결과를 결정합니다. 고정된 $M_{\text{cloud}}$, $n_{\text{core}}$, 그리고 $\epsilon$에 대해:
   • 균일 ($\rho \propto r^0$) 과 완만한 ($\rho \propto r^{-1}$) 프로파일은 재붕괴합니다.
   • 가파른 프로파일 ($\rho \propto r^{-2}$) 은 껍질이 밀도 기울기를 따라 가속화되어 계속 팽창하도록 합니다.
   • Bonnor–Ebert 구는 균일한 구름이 분산할 수 있는 경우에도 재붕괴하거나, 훨씬 나중에 분산합니다.
3. 분산 시간尺度: 분산이 보장되는 모델 그리드에서 Bonnor–Ebert 구름은 균일한 구름보다 약 55% 늦게 분산합니다. 이 체계적인 편차는 BE 구조가 동일한 질량을 더 큰 물리적 반지름에 분산시켜 껍질이 더 먼 거리를 이동해야 하기 때문에 발생합니다.
4. 초신성 이전 우세성: 테스트된 매개변수 공간에서 구름의 분산 또는 재붕괴는 첫 번째 핵붕괴 초신성이 도달하기 전 (보통 >3 Myr) 에 초신성 이전 피드백 (바람, 복사, $P_{\text{H II}}$) 에 의해 결정됩니다.
5. LyC 탈출: 코드는 Lyman-continuum (LyC) 탈출 분율 ($f_{\text{LyC}}^{\text{esc}}$) 을 역학적 비움과 별개의 "시계"로 추적하여, 피드백이 껍질 역학을 변경하는 시간과 광자가 희박한 성간 매질 (ISM) 로 누출되기 시작하는 시간을 구분합니다.

의의 및 주장
저자들은 Trinity 가 구름 매개변수 공간 전반에 걸친 피드백 우세성을 매핑하기 위한 매우 효율적이고 해석 가능한 프레임워크를 제공한다고 주장합니다. 주요 결론은 피드백 우세성을 항성 집단만으로 추론할 수 없다는 것입니다. 결과는 다음에 결정적으로 의존합니다:

• 피드백 채널이 껍질과 어떻게 결합하는지 (특히 $P_{\text{H II}}$의 역할).
• 뜨거운 버블 냉각 전이의 시기.
• 태생 구름의 밀도 구조.

본 논문은 정형화된 균일 구름이나 단일 멱함수 기울기를 가정하는 것이 추정된 출현 시간과 분산 임계값에 체계적인 오차를 도입한다고 주장합니다. 유연한 밀도 프로파일과 함께 바람, 복사, 그리고 광이온화를 자기 일관적으로 결합하는 도구를 제공함으로써, Trinity 는 분해된 H II 영역의 연구를 가능하게 하고, 이온화 광자 누출 및 구름 분산에 관한 은하 규모 시뮬레이션을 위한 서브-그리드 처방을 제공합니다.

한계점
저자들은 명시적으로 1 차원 구대칭이 각도적 구조 (예: 샴페인 흐름, 물집 파열) 를 평균화하며, 냉각을 증대시킬 수 있는 켈빈 - 헬름홀츠 또는 레이리 - 테일러와 같은 불안정성을 해결할 수 없다고 지적합니다. $P_{\text{H II}}$의 처리는 버블 - 이온화면 시스템을 단일 반지름으로 사용하며, $R_b$와 $R_{\text{if}}$를 별도의 역학적 변수로 진화시키지 않습니다. 또한, 현재 모델은 자기 지지, 우주선 압력 기울기, 그리고 덮개 분율 (누출 껍질) 의 진화를 무시합니다. 이러한 요소들은 향후 개발의 대상으로 식별되었습니다.