Galaxy assembly bias in cosmological hydrodynamical simulations: a comparison between SIMBA and IllustrisTNG

1. 문제 제기 (Problem)

• 배경: 현대 우주론에서 은하의 공간 분포는 암흑물질의 대규모 구조를 추적하는 핵심 관측량입니다. 은하 - 헤일로 연결을 모델링하는 가장 기본적인 접근법인 **헤일로 점유 분포 (HOD, Halo Occupation Distribution)**는 은하의 점유율이 오직 헤일로 질량에만 의존한다고 가정합니다.
• 문제점: 그러나 헤일로의 질량 외의 속성 (형성 시간, 농도, 스핀, 우주적 환경 등) 이 은하 형성에 영향을 미친다면 GAB 가 발생하여 은하 군집 예측에 체계적인 편향을 일으킬 수 있습니다.
• 연구 필요성: GAB 의 물리적 기원과 그 크기는 여전히 불분명하며, 관측 데이터에 대한 해석에 차이를 낳고 있습니다. 특히, 서로 다른 물리 모델 (항성 및 활동성 은하핵 (AGN) 피드백 등) 을 적용한 시뮬레이션들 간의 GAB 진화 양상이 어떻게 다른지, 그리고 어떤 물리 과정이 GAB 를 주도하는지 규명하는 것이 시급합니다.

2. 방법론 (Methodology)

저자들은 두 가지 최첨단 시뮬레이션 (SIMBA, IllustrisTNG) 의 데이터를 비교 분석하기 위해 다음과 같은 기법을 사용했습니다.

• 데이터 소스:
  • SIMBA: MFM(Meshless Finite Mass) 유체역학 솔버 사용, 항성 및 AGN 피드백 모델 적용.
  • IllustrisTNG (TNG300-1): AREPO 이동 격자 솔버 사용, 다양한 AGN 피드백 모드 적용.
  • 두 시뮬레이션 모두 관측된 은하 질량 함수를 재현하지만, 피드백 메커니즘과 가스 역학 처리 방식에서 큰 차이가 있습니다.
• GAB 측정 (Shuffling 기법):
  • 헤일로를 질량 구간 (0.1 dex 간격) 으로 나눈 후, 각 질량 구간 내에서 은하를 무작위로 다른 헤일로에 재배치 (Shuffle) 하는 과정을 수행했습니다.
  • 원본 군집 vs 셔플링된 군집: 셔플링은 헤일로의 질량 외 모든 속성 (내부적, 환경적) 과 은하의 연결 관계를 끊습니다.
  • GAB 신호 정의: $[w_p(r_p) / w_p(r_p)_{shuffled} - 1]$로 정의하며, 여기서 $w_p(r_p)$는 투영된 상관 함수입니다. 양의 값은 GAB 가 은하 군집을 강화함을 의미합니다.
• 환경 변수 분석:
  • 과밀도 ($\delta_5$): 5 $h^{-1}$Mpc 반경으로 평활화된 물질 과밀도.
  • 조석 이방성 ($\alpha_5$): 조석 텐서의 고유값을 기반으로 한 환경의 비등방성.
  • 더블 셔플링 (Double Shuffling): 헤일로 질량뿐만 아니라 $\delta_5$ 또는 $\alpha_5$ 구간 내에서도 은하를 재배치하여, 특정 환경 변수가 GAB 에 기여하는 정도를 분리 측정했습니다.
• HOD 분석:
  • 과밀도 구간 (상위 25%, 하위 25%) 에 따라 헤일로 질량별 중심은하와 위성은하의 점유 수 ($N_{cen}, N_{sat}$) 변화를 분석했습니다.
  • 저질량 및 고질량 헤일로로 나누어 GAB 기여도를 분리했습니다.

3. 주요 결과 (Key Results)

A. 적색편이에 따른 GAB 진화의 분기

• 공통점: 두 시뮬레이션 모두 $z=5$에서는 GAB 가 거의 0 이었으며, $z=2$까지 약 5% 수준으로 증가했습니다. 이는 항성 피드백이 지배적인 시기입니다.
• 분기점 ($z < 2$):
  • SIMBA: $z=2$ 이후 GAB 가 급격히 감소하여 $z=0$에서 거의 0 이 됩니다.
  • IllustrisTNG: $z=2$ 이후 GAB 가 계속 증가하여 $z=0$에서 약 **10%**에 달합니다.
• 원인: 이 차이는 두 모델의 AGN 피드백 구현 방식 차이에서 기인합니다. SIMBA 의 안정된 제트 (jet) 피드백은 그룹 규모 헤일로의 가스 제거에 더 효과적이어서 GAB 를 억제하는 반면, TNG 의 랜덤화된 피드백은 GAB 를 유지하거나 증폭시킵니다.

B. 환경 변수의 기여도

• 과밀도 ($\delta_5$): 두 시뮬레이션 모두에서 GAB 의 대부분을 설명합니다. 더블 셔플링 결과, $\delta_5$를 고정하면 GAB 신호가 크게 감소했습니다.
• 조석 이방성 ($\alpha_5$):
  • SIMBA: GAB 에 거의 기여하지 않습니다.
  • TNG: 소량이지만 무시할 수 없는 기여를 합니다.
  • 결론적으로, GAB 의 주된 원인은 국소적인 과밀도 환경이며, 조석 이방성의 영향은 미미합니다.

C. HOD 및 저질량 헤일로의 역할

• 저질량 헤일로의 지배적 기여: GAB 신호의 대부분은 저질량 헤일로 ($M_h < 10^{11.5} M_\odot$) 에서 기원합니다. 고질량 헤일로는 환경에 관계없이 거의 일정하게 은하를 보유합니다.
• 환경 의존성: 고밀도 ($\delta_5$가 높은) 환경에 위치한 저질량 헤일로일수록 중심은하와 위성은하를 보유할 확률이 높아집니다.
• 은하 질량 임계값 변화: TNG 에서 은하 질량 임계값을 높여 ($10^{10} M_\odot$) 분석을 반복한 결과, GAB 의 진화 경향은 유지되지만 전체 크기는 약간 감소했습니다. 여전히 GAB 는 저질량 헤일로의 점유율 조절을 통해 발생합니다.

D. HOD 파라미터화 제안

• 저자들은 HOD 가 평균값에서 벗어나는 정도를 과밀도 ($\delta_5$) 에 따라 모델링하는 간단한 파라미터화 식을 제안했습니다.
• 특히, 중심은하와 위성은하의 '턴온 (turn-on)' 질량에 해당하는 기울기 ($\mu, \nu$) 가 밀도 순위 ($r_\delta$) 에 따라 선형적으로 변한다고 가정하여, HOD 기반의 은하 모의 (mock) 생성 시 GAB 효과를 반영할 수 있는 수식을 제시했습니다.

4. 주요 기여 및 의의 (Significance)

1. 물리적 기원에 대한 통찰: 서로 다른 유체역학 시뮬레이션 간의 GAB 진화 차이를 비교함으로써, AGN 피드백이 GAB 의 세기와 진화를 결정하는 핵심 물리 과정임을 규명했습니다. 이는 GAB 가 단순한 통계적 현상이 아니라, 은하 형성 물리 (특히 피드백) 와 밀접하게 연관되어 있음을 보여줍니다.
2. 환경 요인의 정량화: GAB 가 주로 **과밀도 ($\delta_5$)**에 의해 주도되며, 조석 이방성의 영향은 제한적임을 명확히 했습니다. 이는 기존 관측 연구들 간의 논쟁 (어떤 환경 변수가 GAB 를 주도하는가) 에 중요한 제약을 제공합니다.
3. 저질량 헤일로의 중요성 강조: GAB 는 주로 저질량 헤일로의 환경 의존적 점유율 변화에서 비롯된다는 사실을 확인했습니다. 이는 고질량 헤일로는 환경에 무관하게 행동하지만, 저질량 헤일로는 환경에 민감하게 반응함을 의미합니다.
4. 차세대 관측을 위한 모델링 개선: 정밀 우주론 (DESI, Euclid 등) 을 위한 은하 군집 모델링에서 GAB 를 무시하면 편향이 발생할 수 있습니다. 저자들이 제안한 밀도 의존적 HOD 파라미터화는 보다 정확한 은하 모의 (mock) 카탈로그 생성을 가능하게 하여, 관측 데이터 해석의 신뢰성을 높이는 데 기여합니다.

요약

이 연구는 SIMBA 와 IllustrisTNG 시뮬레이션을 비교하여 GAB 가 과밀도 환경에 의해 주도되며, 그 세기와 진화는 AGN 피드백 모델에 크게 의존함을 밝혔습니다. 특히 GAB 는 저질량 헤일로에서 주로 발생하며, 이를 정량화한 HOD 수정 모델을 제안함으로써 미래 정밀 우주론 관측 데이터의 정확한 해석을 위한 중요한 토대를 마련했습니다.