얇은 왜소 은하의 풍부함이 시뮬레이션을 위협한다

초록

**
관측된 은하들의 축비율 분포를 이용해 $10^{7}!-!10^{11},M_{\odot}$ 범위에서 얇은( $c/a<0.2$) 은하가 차지하는 비중을 추정하였다. $10^{9}!-!10^{10},M_{\odot}$ 은하의 40 %가 $c/a<0.2$이며, $10^{8},M_{\odot}$ 수준에서도 30 %가 같은 정도로 얇다. 반면 최신 TNG50, FIREbox, Romulus25 시뮬레이션은 $M_{\star}<10^{9},M_{\odot}$ 이하에서는 $c/a<0.2$ 은하를 전혀 만들지 못한다. 이는 해상도와 피드백 구현의 한계가 원인일 가능성을 시사한다.

**

상세 분석

**
본 연구는 GAMA, DESI, ALFALFA, 그리고 근거리 은하 카탈로그에서 추출한 10⁷–10¹¹ M⊙ 범위의 은하 1천만 개 이상을 대상으로 투영 축비율 $q=b_{\rm p}/a_{\rm p}$ 를 측정하고, 무작위 시점 가정을 통해 3차원 축비율 $c/a$ 의 내재 분포를 역추정하였다. 핵심 방법론은 (1) 고정된 관측 파장(주로 $r$‑밴드)에서의 Sersic 모델 적합을 통한 $q$ 확보, (2) $q$ 분포의 저값(예: $q<0.2$)이 실제 얇은 원반(oblates) 혹은 길쭉한 프로레이트(prolates) 구조의 존재를 강하게 암시한다는 점을 이용한 최소 얇은 은하 비율 추정, (3) 베이지안 혹은 최대우도 방식으로 $c/a$ 확률밀도 함수를 파라미터화하고 관측 $q$와 비교하여 최적 파라미터를 도출하였다.

관측 결과는 질량 의존적인 얇은 은하 비율을 명확히 보여준다. $10^{9}!-!10^{10},M_{\odot}$ 구간에서는 $c/a<0.2$ 은하가 전체의 약 40 %를 차지하고, $c/a<0.3$ 은 80 %에 달한다. 질량이 $10^{8},M_{\odot}$ 수준으로 내려가면 비율이 약간 감소하지만 여전히 30 %($c/a<0.2$)와 65 %($c/a<0.3$) 수준으로 높다. 반대로 $M_{\star}>10^{10.5},M_{\odot}$에서는 얇은 원반 비율이 급격히 떨어져, 전형적인 타원형 혹은 판상 구조가 지배한다.

시뮬레이션과의 비교에서는 TNG50‑1, FIREbox, Romulus25 세 모델 모두 $M_{\star}<10^{9},M_{\odot}$ 이하에서는 $c/a<0.2$ 은하를 전혀 생성하지 못한다는 점이 두드러진다. 이는 (i) 현재 해상도(별 입자 질량 $∼10^{5},M_{\odot}$ 이하)와 중력 연화 길이(∼100 pc)가 얇은 원반의 수직 구조를 충분히 포착하지 못함, (ii) 피드백 구현(에너지 주입, 바람, 방사선 압력 등)이 과도하게 난류와 수직 팽창을 유발해 원반을 두껍게 만든다는 물리적 가정, (iii) 초기 각운동량 분포와 은하 외부 환경(예: 서브헬로와의 충돌, 가스 공급)의 모델링 차이가 실제 관측과 일치하지 않음 등을 시사한다.

또한 논문은 관측적인 편향을 최소화하기 위해 HI‑선 선택(ALFALFA)과 광학 선택(GAMA, DESI) 두 가지 서로 다른 샘플을 결합했으며, 각각이 얇은 은하를 과소·과대 추정할 가능성을 논의한다. 특히 HI‑선 샘플은 면역적인 면에서 가장 얇은 원반을 놓치기 쉬우나, 결과적으로도 얇은 은하 비율이 높은 것을 확인했다.

결론적으로, 얇은 왜소 은하가 실제 우주에서 흔히 존재한다는 증거는 현재 ΛCDM 기반 수치 시뮬레이션이 아직 작은 질량 스케일에서 디스크 형성 메커니즘을 정확히 재현하지 못하고 있음을 강력히 시사한다. 향후 고해상도(∼10 pc)와 보다 정교한 피드백(다중 위상, 방사선, 화학적 냉각) 모델을 포함한 시뮬레이션이 필요하다.

**