Resolving Star Cluster Formation in Galaxy Simulations with Cosmic Ray Feedback

기술 요약: 우주선 피드백을 통한 은하 시뮬레이션 내 별 클러스터 형성 문제 해결

문제 정의
별 클러스터는 거대 질량 별 형성이 일어나는 주요 지점이며, 이는 성간 물질(ISM)을 형성하고 향후 별 형성을 조절하는 피드백 과정을 주도한다. 은하 피드백에 대한 포괄적인 이해를 위해서는 파섹(parsec) 규모의 클러스터 형성과 다상(multiphase) 성간 물질을 분해능 있게 구현해야 한다. 우주선(Cosmic Rays, CRs)이 성간 물질의 압력과 가열에 기여한다는 사실은 알려져 있으나, 이것이 별 클러스터 형성률(SFRs), 클러스터 질량 함수 및 클러스터의 결합도(boundedness)에 미치는 구체적인 영향은 여전히 불확실하다. 기존 연구들은 은하 유출(outflows)의 맥락에서 우주선 피드백을 탐구하거나 클러스터에 대한 하위 격자(subgrid) 모델을 사용해 왔으나, 개별 별 입자, 동역학적으로 결합된 우주선, 그리고 피드백에 의해 조절되는 다상 성간 물질을 동시에 분해능 있게 구현한 연구는 거의 없다. 본 연구는 특히 자기 수렴 이론(self-confinement theory)에 기반한 우주선 수송 모델이 별 클서터의 형성과 특성에 어떻게 영향을 미치는지에 대한 이해의 공백을 다룬다.

방법론
저자들은 Arepio 이동 격자 코드(moving-mesh code)를 사용하여 수행된 고해상도 "tallbox" 시뮬레이션을 분석한다. 시뮬레이션에는 수소와 금속의 냉각/가열 과정을 추적하는 Crisp 비평형 열화학 및 피드백 모델이 적용된다. 연구는 우주선 수송의 효과를 격리하기 위해 세 가지 뚜렷한 물리 케이스를 활용한다:

1. MHD: 우주선을 균일한 가열/이온화 항으로 취급하는 자기유체역학(Magnetohydrodynamics).
2. CR-NL: 비선형 랜다우 감쇠(Nonlinear Landau Damping, NLLD)를 유일한 감쇠 기제로 사용하는 2-모멘트 우주선 자기유체역학(CRMHD). 이 경우 우주선과 가스의 결합이 균일하게 나타난다.
3. CR-NL-IN (Fiducial): NLLD와 이온-중성자 감쇠(Ion-Neutral Damping, IND)를 모두 포함하는 CRMHD. 이 케이스에서 우주선은 차가운 중성 가스에서 탈착되어 국지적인 영향력이 감소한다.

시뮬레이션은 질량 분해능 $10 M_\odot$, 공간 스케일 $\sim 0.5$ pc까지 가스를 분해한다. 별 형성은 가스 셀이 밀도 임계값($n_H = 10^3 \text{ cm}^{-3}$)을 초과할 때 국지적 효율 100%로 충돌 없는 별 입자로 변환되는 슈미트(Schmidt) 유형의 처방을 통해 구현된다. 항성 피드백은 초신성(SNe)을 통해서만 모델링되며, 복사 피드백과 항성풍은 우주선과 초신성의 영향만을 격리하기 위해 제외되었다. 초신성은 3~38 Myr 된 별들에 대해 확률적으로 주입되며, CRMHD 케이스에서는 초신성 에너지의 5%가 추가로 우주선으로 주입된다.

별 클러스터는 시간 좌표를 통합한 수정된 Friends-of-Friends (FoF) 알고리 알고리즘을 사용하여 식별되며, 이를 통해 식별된 그룹이 2 pc의 공간 연결 길 내에 있는 동시 생성(coeval) 별들임을 보장한다. 분석은 초기 특성을 규명하기 위해 티달 파쇄(tidal disruption)가 일어나기 전인 "젊은" 클러스터(연령 $\sim 10$ Myr)에 집중한다.

주요 결과

1. 클러스터 질량 함수: 세 가지 시뮬레이션 케이스 모두 질량 $> 10^3 M_\odot$ 영역에서 $-2.0$에서$-2.4$사이의 기울기를 갖는 멱법칙($dN/dM \propto M^{\alpha_{PL}}$)을 따르는 별 클러스터 질량 함수를 생성한다. 이는 질량 함수 기울기가 별 형성율 표면 밀도($\Sigma_{SFR}$)와 양의 상관관계를 갖는 관측 트렌드와 일치한다. 우선선(CRs)의 포함은 관측된 관계를 따라 데이터 포인트를 이동시키지만, 비물리적인 결과를 생성하지는 않는다.
2. 클러스터 형성 효율 (CFE): 거대 클러스터($M > 5000 M_\odot$)에 대해 계산된 CFE($\Gamma$)는 이론적 예측(예: Kruijssen 2012) 및 관측과 일치하게 $\Sigma_{SFR}$와 양의 상관관계를 보인다. 그러나 더 낮은 질량 임계값($M > 1000 M_\odot$)의 경우, CRMHD 케이스는 높고 일정한 CFE를 보이는 반면, MHD 케이스는 더 낮은 값을 보인다. 저자들은 CRMHD 케이스의 높은 CFE가 밀도 임계값과 효율을 높게 설정한 처방의 인위적 결과(artifact)일 수 있다고 제안한다.
3. 초신성 클러스터링: 세 케이스 모두에서 약 50%의 초신성이 희박한 사전 조건화된 가스($n_H < 10^{-2} \text{ cm}^{-3}$)에서 발생한다. 이 비율은 우선선 모델과 무관하게 일정하며, 이는 우주선 피드백이 전역적인 희박 성간 물질의 부피 분율을 변화시킴에도 불구하고, 첫 몇 개의 초신성이 폭발하는 국지적 환경을 크게 변화시키지는 않음을 나타낸다.
4. 클러스터 특성 및 비리알 매개변수:
   • 반경 및 속도 분산: MHD 케이스의 별 클러스터는 CRMHD 케이스에 비해 더 큰 반질량 반경과 높은 1D 속도 분산을 보인다.
   • 결합도 (Boundedness): MHD 클러스터는 평균적으로 결합되지 않은 상태($\alpha_{vir} > 1$)인 반면, CRMHD 클러스터는 평균적으로 결합된 상태($\alpha_{vir} < 1$)이다.
   • 메커니즘: 저자들은 이러한 차이를 성간 물질의 난류 에너지 예산(turbulent energy budget)의 차이로 돌린다. MHD 케이스는 별 형성율(SFR)이 가장 높아 초신성 에너지 주입률이 높고 성간 물질이 더 난류적이다. 이 높은 주변 난류가 형성되는 분자 구름에 상속되어, 결과적으로 더 높은 속도 분산을 가진, 중력적으로 덜 결합된 별 클러스터를 만든다. 우선선 피드백은 SFR을 감소시켜 난류 에너지 예산을 낮춤으로써, 더 단단히 결합된 클러스터를 생성한다.

의의 및 주장
본 논문은 동역학적으로 결합된 우선선, 다상 성간 물질, 그리고 거대 질량 별을 대표하는 개별 별 입자를 동시에 포함하는 최초의 별 클러스터 형성 시뮬레이션을 제시한다고 주장한다. 주요 의의는 다음과 같다:

• Crisp 피드백 모델은 초신성 이전의 복사 피드백이 없는 상태에서도 미세 조정 없이 관측된 질량 함수 기울기와 환경 의존성을 재현할 수 있음을 보여준다.
• 우선선 피드백의 주요 영향은 전역적인 SFR의 감소이다. 이러한 감소는 성간 물질의 난류 에너지 예산을 낮추며, 이는 결과적으로 효과적인 우선선 피드백이 없는 경우보다 더 결합된(낮은 속도 분산) 별 클러스터를 생성한다.
• 우선선은 개별 밀집 분자 구름의 붕괴를 동역학적으로 방해하거나 초신성의 국지적 사전 조건화(preconditioning)를 크게 변화시키지 않으며, 그 효과는 성간 물질의 전역적 조절을 통해 매개된다.

저자들은 본 연구의 설정이 복사 피드백과 항성풍을 생략한 이상화된 형태임에도 불구하고, 우선선의 영향에 관한 각 케이스 간의 상대적 차이는 견고하다고 결론짓는다. 이 결과는 우선선 수송 모델에 대한 새로운 제약을 제공하며, 우선선 피드백이 새로운 미관측 영역을 만드는 것이 아니라 관측된 별 클러스터 관계를 따라 시뮬레이션 결과를 이동시키는 역할을 한다는 점을 시사한다.