은하 형성

계산물리 시뮬레이션 물리 마지막 수정: 2026-09-05 04:49:12

1. 개요[편집]

은하 형성 Galaxy Formation
표준 틀	ΛCDM 계층적 병합 — 화이트-리스(1978)의 2단계 그림
1단계	암흑물질 헤일로가 중력으로 붕괴·병합
2단계	헤일로 안의 기체가 복사 냉각으로 가라앉아 별이 됨
중심 난제	과냉각 — 되먹임이 없으면 별이 너무 많이 생긴다
두 갈래 방법	준해석적 모형(SAM) · 유체 시뮬레이션(SPH · AMR · 이동격자)
주요 제약	은하 항성질량함수 · 툴리-피셔 · 질량-금속량 관계 · 우주 SFR 역사

중력은 쉽다. 기체가 어려운 게 아니라, 기체가 별이 되는 과정을 아무도 분해하지 못하는 것이 어렵다.

은하 형성(galaxy formation)은 초기 우주의 미세한 밀도 요동이 중력으로 성장해 암흑물질 헤일로를 만들고, 그 헤일로 안에 갇힌 바리온 기체가 냉각·수축해 별을 만들며, 그 별과 중심 블랙홀이 되먹임으로 다시 기체를 데우고 밀어내는 순환 전체를 다루는 분야다. 계산과학의 관점에서 이 문제가 특별한 이유는 딱 하나다 — 동역학 범위가 미쳤다. 우주 규모(Mpc)에서 별 형성 영역(pc)까지 여섯 자릿수 이상, 밀도로는 열 자릿수 이상을 한 계산에 담아야 한다.

이 문서는 무엇을 풀어야 하는가, 즉 물리 모듈과 그 교정 문제에 초점을 둔다. 팽창 상자·좌표계·중력 솔버·수렴 판정 같은 계산 인프라는 우주론 시뮬레이션이, 암흑물질 자체의 동역학과 소규모 문제들은 암흑 물질이 다루므로 여기서는 반복하지 않는다.

2. 계층적 그림 — 아래에서 위로[편집]

우주 마이크로파 배경이 보여 주는 초기 밀도 대비는 $\delta \sim 10^{-5}$ 수준이다. 이것이 물질 우세기에 $\delta \propto a$ 로 성장하고, $\delta \sim 1$ 을 넘으면 선형 이론이 깨지며 비선형 붕괴로 들어간다. 차가운 암흑물질의 파워 스펙트럼은 작은 스케일에서 더 큰 진폭을 가지므로 작은 것이 먼저 만들어지고, 병합으로 커진다. 이것이 계층적(bottom-up) 구조 형성이다.

화이트와 리스(1978)가 여기에 바리온을 얹은 2단계 그림을 세웠다.

암흑물질 헤일로가 붕괴한다. 이건 무충돌 문제이고, 원리적으로 해가 유일하며 해상도를 올리면 수렴한다.
헤일로의 퍼텐셜 우물에 갇힌 기체는 바이리얼 충격을 지나 $T_{\rm vir}$ 로 데워진 뒤, 복사로 에너지를 잃으면 가라앉는다. 각운동량이 있으면 원반이 되고, 밀도가 충분해지면 별이 된다.

여기서 도구가 되는 것이 병합나무(merger tree)다. 암흑물질 시뮬레이션의 각 시점 스냅숏에서 헤일로 찾기(FOF + 부구조 식별)를 하고, 시점 사이를 입자 ID로 연결해 “누가 누구의 조상인가” 그래프를 만든다. 이 자료구조가 준해석적 모형의 뼈대가 되고, 유체 시뮬에서도 은하의 성장 이력을 추적하는 데 그대로 쓰인다. 나무를 만드는 일 자체가 만만치 않아서 — 병합 직전에 부헤일로가 일시적으로 사라졌다 나타나는 등 — 나무 구성 알고리즘이 결과 통계를 흔든다는 보고가 꾸준히 있다.

3. 냉각 함수 — 은하를 만들 수 있게 하는 것[편집]

기체가 별이 되려면 먼저 에너지를 버려야 한다. 단위 부피당 복사 손실률을 다음처럼 쓰는 것이 관례다.

\mathcal{L} = n_{\rm H}^2\, \Lambda(T, Z, \ldots) \ [\mathrm{erg\,cm^{-3}\,s^{-1}}]

$\Lambda$ 가 냉각 함수이며, 이 곡선의 모양이 은하 형성의 많은 것을 결정한다.

$T \lesssim 10^4$ K 에서 원시 조성 기체는 사실상 냉각하지 못한다. 수소를 들뜨게 할 만큼 뜨겁지 않기 때문이다. 이 구간을 뚫으려면 수소 분자( $\mathrm{H_2}$ )나 금속·먼지가 필요하고, 그래서 최초의 별 형성이 $\mathrm{H_2}$ 화학에 목을 매는 것이다.
$10^4 \lesssim T \lesssim 10^6$ K 은 수소·헬륨의 선 냉각이 지배하며 $\Lambda$ 가 급격히 솟는다.
$T \gtrsim 10^7$ K 에서는 완전 이온화되어 제동복사만 남고 $\Lambda \propto T^{1/2}$ 로 완만해진다.
금속이 섞이면 $10^5$ – $10^7$ K 구간의 냉각률이 크게 올라간다. 은하가 별을 만들수록 다음 세대 기체가 더 잘 식는다는 양의 되먹임이 여기 숨어 있다.

여기에 자외선 배경복사가 광이온화 가열로 반대 방향의 항을 넣는다. 실무 코드는 $\Lambda$ 를 온도·밀도·금속량·적색편이의 4차원 보간표로 넣는다(CLOUDY 등으로 미리 계산). 평형 가정이 정당한지, 자기 차폐를 어떻게 근사할지가 이 표의 주요 불확실성이다.

냉각 시간과 자유낙하 시간의 비가 판정 기준이 된다.

t_{\rm cool} \simeq \frac{3 n k_B T}{2 n_{\rm H}^2 \Lambda}, \qquad t_{\rm ff} \simeq \sqrt{\frac{3\pi}{32 G \rho}}

$t_{\rm cool} < t_{\rm ff}$ 이면 기체가 압력 지지를 잃고 붕괴한다. 리스-오스트라이커와 실크가 1977년에 이 조건으로 은하의 특징적 질량 스케일을 논증했다. 관련해 냉각의 두 양상이 구분된다 — 무거운 헤일로에서는 기체가 바이리얼 충격을 지나 뜨거운 헤일로를 이루고 천천히 식는 뜨거운 모드, 가벼운 헤일로와 고적색편이에서는 충격이 안정하지 못해 차가운 흐름이 필라멘트를 타고 중심까지 곧장 들어오는 차가운 모드다. 후자는 시뮬레이션이 먼저 발견하고 나중에 해석적으로 정리된 대표 사례다.

4. 과냉각 — 이 분야를 만든 실패[편집]

여기가 은하 형성이 “중력 문제”에서 “미해결 문제”로 바뀌는 지점이다.

냉각 함수만 넣고 시뮬레이션을 돌리면 거의 모든 바리온이 식어서 별이 된다. 관측된 우주는 그렇지 않다 — 우주 바리온의 대부분은 여전히 기체이고, 헤일로 질량 대비 항성질량 비율은 어디에서도 우주 바리온 분율에 도달하지 못한다. 이 불일치가 과냉각 문제(overcooling problem, cooling catastrophe)다. 1990년대 초 나바로·벤츠·화이트 등의 초기 SPH 계산이 이 실패를 정면으로 보고했고, 동반 증상으로 각운동량 대재앙 — 기체가 너무 일찍 뭉쳐 작은 덩어리로 병합하면서 각운동량을 암흑물질에 넘겨줘, 만들어진 원반이 관측보다 압도적으로 작고 중심집중적이 되는 문제 — 가 따라왔다.

정량적으로 요약하면 이렇다. 존재비 맞추기(abundance matching)로 관측된 은하 항성질량함수와 헤일로 질량함수를 이으면, 별 형성 효율이 헤일로 질량에 대해 뚜렷한 봉우리를 그린다.

헤일로 질량	별 형성 효율	억제 주범
$\lesssim 10^{11}\,M_\odot$	매우 낮음	초신성 바람 · 재이온화
$\sim 10^{12}\,M_\odot$	최대 (바리온의 수십 % 수준)	—
$\gtrsim 10^{13}\,M_\odot$	다시 낮음	AGN 되먹임

양쪽 끝이 다 눌려 있다는 것이 은하 형성 이론이 설명해야 할 1번 사실이며, 두 억제 기구가 서로 다른 물리라는 것이 이 분야가 두 개의 서브그리드 모형을 동시에 안고 가는 이유다.

5. 되먹임 — 서브그리드가 사는 곳[편집]

초신성 되먹임. 얕은 퍼텐셜 우물에서 기체를 밖으로 날리는 것이 목표다. 수치적 문제는 냉혹하다. 초신성 잔해가 운동량을 실어 나르는 것은 세도프-테일러 단계인데, 그 단계를 분해하려면 별 형성 영역에서 pc 이하 해상도가 필요하다. 대규모 시뮬레이션의 해상도는 그보다 한참 위이고, 분해하지 못한 채 열에너지를 그냥 넣으면 그 기체가 즉시 냉각해 버려서 아무 일도 일어나지 않는다. 에너지를 넣었는데 냉각 함수가 그대로 다시 복사로 버리는 것이다. 처방들이 여기서 갈린다.

지연 냉각 — 에너지를 받은 입자의 냉각을 일정 시간 강제로 끈다. 조잡하지만 효과적이고, 오래 쓰였다.
확률적 열 주입 — 에너지를 넓게 퍼뜨리는 대신 소수 입자에만 몰아 큰 온도 상승( $\Delta T \sim 10^{7.5}$ K 급)을 준다. 그 온도에서는 냉각 시간이 충분히 길어져 팽창이 일어난다.
분리 바람 — 운동량을 실은 바람 입자를 쏘고 잠시 유체와의 상호작용을 끈다.
종단 운동량 주입 — 세도프 단계를 분해하는 대신, 소규모 계산에서 알려진 잔해가 최종적으로 남기는 운동량(초신성 하나당 대략 $3\times10^5\,M_\odot\,\mathrm{km/s}$ 규모, 밀도·금속량에 약하게 의존)을 직접 넣는다. 물리적으로 가장 정당화하기 좋은 축이며, 근래의 표준에 가깝다.

AGN 되먹임. 무거운 헤일로에서 냉각 흐름을 막고 은하를 소광(quenching)시키는 역할이다. 헤일로에 블랙홀 씨앗을 심고, 본디-호일형 또는 중력토크형 강착률 추정으로 성장시킨 뒤, 강착률이 높으면 열에너지(퀘이사 모드), 낮으면 운동량·제트(전파 모드)를 방출한다. 이것 없이는 은하 항성질량함수의 무거운 쪽 지수적 절단이 재현되지 않는다. 관측된 블랙홀 질량-벌지 속도분산 관계( $M_\bullet$ – $\sigma$ )도 자기조절 되먹임의 자연스러운 귀결로 나온다.

6. 두 갈래 접근[편집]

준해석적 모형(SAM)은 병합나무 위에 냉각·별 형성·되먹임·화학진화를 상미분방정식 다발로 얹는다. 유체를 풀지 않아 비용이 사실상 0이고, 그래서 매개변수 공간을 MCMC로 훑을 수 있다. 유체 시뮬레이션은 기체를 실제로 풀어 형태·내부구조·기체 분포를 직접 얻는다. 둘의 비용·성격 비교와 대표 프로젝트 표는 우주론 시뮬레이션에 정리되어 있고, 여기서는 은하 물리 쪽에서 갈리는 지점만 짚는다.

유체 기법의 선택은 은하 규모에서 실질적인 차이를 만든다.

SPH — 라그랑주식이라 밀도 적응이 공짜다. 고전 형식은 접촉 불연속면에서 유체 불안정을 억제하는 결함이 있었고, 이것이 은하 헤일로의 뜨거운 기체와 차가운 흐름의 혼합을 과소평가하는 문제로 직결됐다. 압력-엔트로피 형식과 개선된 인공점성으로 상당 부분 고쳤다.
AMR 오일러 격자 — 적응 격자 세분화 + 유한체적법. 충격 포착이 깔끔한 대신, 빠르게 움직이는 은하에서 이류 오차가 쌓이고 갈릴레이 불변성이 완전하지 않다.
이동 격자 · 무격자 유한질량 — 보로노이 다이어그램 격자를 유체와 함께 흘리거나(AREPO), 라그랑주 유한질량으로 리만 문제를 푸는 방식(GIZMO). 두 진영의 장점을 합치려는 시도이고, 최근 대형 프로젝트가 많이 채택한다.

여기에 은하 규모에서만 추가로 필요한 모듈이 붙는다 — 항성 진화와 화학 원소 방출(초신성 II형·Ia형·AGB별의 수율표), 자기유체, 우주선, 복사 되먹임. 완전한 복사 전달 방정식을 푸는 것은 여전히 대부분의 대형 시뮬레이션에서 감당 불가라, 감쇠 근사나 사후처리로 대체한다.

7. 교정, 그리고 재현성[편집]

이 분야의 정직한 자기소개는 이렇다. 서브그리드 매개변수는 유도된 것이 아니라 관측에 맞춰 교정된 것이다. 대표 프로젝트들이 논문에서 이 사실을 명시한다 — EAGLE(2015, SPH 기반, 100 cMpc 상자)은 되먹임 교정 절차를 논문의 한 절로 따로 서술했고, IllustrisTNG(2018, 이동격자, TNG50/100/300 상자군)도 $z=0$ 의 항성질량함수와 크기 관계에 맞춰 조정했음을 밝힌다. 그 밖에 확대(zoom-in)로 개별 은하를 고해상도로 파고드는 FIRE 계열, RAMSES 기반 AMR인 Horizon-AGN, GIZMO 기반 SIMBA 등이 각자 다른 처방으로 같은 관측을 맞춘다.

여기서 논리적 귀결이 따라온다.

교정에 쓰인 관측량은 예측이 아니다. 은하 항성질량함수를 맞춰 놓고 “우리 시뮬이 질량함수를 재현한다”고 말하는 것은 동어반복이다. 예측력은 교정에 쓰이지 않은 양 — 고적색편이 은하 개수, 은하 내부 구조와 형태, 은하단 기체의 X선 프로파일, 흡수선 통계 — 에서 평가해야 한다.

부정직해서 이런 게 아니라, 문제 구조가 그렇다. 그리고 여기서 이 분야 특유의 골칫거리가 하나 더 나온다. 서브그리드 매개변수가 해상도에 의존한다. 그래서 해상도를 올리면 결과가 바뀌고, 매개변수를 다시 맞추면 같아진다(약한 수렴). 수치해석 표준으로는 불편한 상태이며, 우주론 시뮬레이션 문서의 수렴 논의를 참고할 것.

재현성 쪽 대응은 코드 비교 프로젝트다. 같은 초기조건을 여러 코드에 주고 결과를 비교하는 방식으로, 은하단 규모의 초기 비교(1999)에서 시작해 AGORA 같은 은하 규모 비교로 이어졌다. 교훈은 반복적으로 같았다 — 중력만 있으면 코드들이 잘 일치하고, 냉각과 되먹임을 켜는 순간 갈라진다. 갈라지는 폭이 코드 간 수치 기법 차이보다 서브그리드 처방 차이에서 훨씬 크게 온다는 것도 확인됐다.

최근에는 아예 방향을 틀어, 서브그리드 매개변수를 조직적으로 바꾼 시뮬레이션 수천 개를 생산해 기계학습 대리 모델을 학습시키고, 우주론 매개변수와 되먹임 매개변수를 동시에 추론하려는 프로그램이 등장했다. 매개변수를 없애는 대신 주변화해 버리자는 접근이며, 불확실성 정량화 관점에서는 훨씬 건강한 태도다.

8. 관측이라는 채점표[편집]

시뮬레이션이 맞춰야 할 대표적인 제약들.

은하 항성질량함수(GSMF). 양 끝이 눌린 모양. 앞서 본 효율 곡선의 관측적 원본이다. 대부분 교정에 쓰인다.
툴리-피셔 관계. 나선은하의 광도(또는 바리온 질량)와 회전속도의 좁은 멱함수 관계. 1977년에 발견됐고, 바리온 형태로 쓰면 $M_b \propto v^4$ 에 가깝고 산포가 놀랄 만큼 작다. 시뮬레이션에서 이 관계의 기울기·절편·산포를 동시에 맞추는 것은 지금도 까다롭다. 산포가 작다는 사실 자체가 MOND 진영이 꾸준히 지적하는 지점이기도 하다(암흑 물질 참조).
질량-금속량 관계. 무거운 은하일수록 금속이 풍부하다. 유출(바람)이 얼마나 금속을 실어 나가는지를 직접 제약하므로, 되먹임 처방의 강력한 시험대다.
별 형성 주계열과 우주 SFR 밀도 역사. 우주 전체의 별 형성률은 적색편이 $z\sim2$ 에서 봉우리를 이루고 이후 한 자릿수 이상 감소한다. 봉우리의 위치와 이후 감쇠 기울기를 동시에 맞추는 것이 되먹임의 시점과 세기를 시험하는 축이다.
케니컷-슈미트 관계. 기체면밀도와 별 형성률 면밀도의 관계. 대부분의 코드는 이것을 재현하는 것이 아니라 입력으로 강제한다. 재현이 목표가 되려면 성간물질을 분해해야 한다.
고적색편이 은하 개수. 교정에 쓰이지 않으므로 진짜 예측이다. 근래 적외선 관측이 $z>10$ 에서 예상보다 밝은 은하를 다수 보고하면서 별 형성 효율·초기질량함수·먼지 소광 가정이 동시에 재검토되고 있다.

9. 여담[편집]

이 분야에서 “물리를 더 넣었더니 결과가 나빠졌다”는 경험은 흔하다. 냉각을 정교하게 만들면 과냉각이 심해지고, 해상도를 올리면 밀도가 높아져 또 심해진다. 되먹임 없이 해상도만 올리는 것은 문제를 악화시킨다 — 수치해석 직관이 배신당하는 드문 사례다.
초신성 되먹임 처방의 이름들이 하나같이 물리학 논문 같지 않다. “냉각 끄기”, “바람 입자 분리”, “확률적 가열” — 이것들은 물리 모형이라기보다 수치적 응급처치의 이름이고, 관계자들도 그렇게 부른다. 그럼에도 이 처방들 없이는 우주가 별 천지가 되므로 계속 쓴다.¹
은하 시뮬레이션 영상이 대중적으로 인기 있는 것은 그림이 예뻐서인데, 정작 그 그림에서 가장 눈에 띄는 나선팔 구조는 되먹임 처방에 가장 민감한 부분이라 논문에서 정량적으로 인용되는 일이 거의 없다. 예쁜 것은 슬라이드로, 숫자는 통계량으로 간다.²
시뮬레이션의 “별 입자” 하나는 별 하나가 아니다. 대형 상자 시뮬레이션에서 별 입자 하나의 질량은 $10^4$ – $10^6\,M_\odot$ 급이라, 실제로는 성단 하나 이상을 통째로 대표하는 표본점이다. 그래서 초기질량함수를 가정해 그 입자가 몇 개의 초신성을 낼지, 어떤 금속을 얼마나 뱉을지를 계산해 넣는다.³
은하 하나를 만드는 데 필요한 것을 순서대로 적으면 중력·유체·화학·복사·항성진화·블랙홀 강착이다. 이 목록이 곧 계산천체물리학이 다물리 연성 문제의 전시장이 된 이유이며, 동시에 어느 코드도 전부를 제1원리로 풀지 못한다는 고백의 목록이기도 하다.

10. 관련 문서[편집]

“냉각을 끈다”가 얼마나 과감한 처방인지 생각해 보면 재미있다. 물리 방정식의 한 항을 특정 입자에 대해 일정 시간 동안 강제로 0으로 만드는 것인데, CFD 쪽에서 같은 짓을 했다가는 리뷰어에게 살해당한다. 그럼에도 이게 통용된 이유는 대안이 “해상도를 세 자릿수 올리기”였기 때문이다. 근래 종단 운동량 주입 방식이 표준으로 자리 잡은 것은 이 응급처치를 조금 더 물리적으로 정당화된 근사로 바꾼 결과다. ↩
같은 이유로 이 분야의 논문 그림은 대부분 예쁜 렌더링이 아니라 지루한 산점도와 오차막대다. 화려한 이미지는 보도자료용이고, 심사에서 싸우는 것은 항성질량함수의 저질량 끝 기울기 같은 것들이다. 우주 전체를 렌더링해 놓고 정작 논쟁은 로그-로그 그래프의 기울기 하나로 하는 셈인데, 계산과학에서는 꽤 흔한 광경이다. ↩
이 사실을 놓치면 결과 해석이 통째로 어긋난다. 왜소은하 하나가 별 입자 열 개로 이루어져 있다면 그 은하의 항성질량은 유효숫자 한 자리도 못 되고, 형태 이야기는 아예 성립하지 않는다. 암흑 물질 문서의 “시뮬레이션 입자는 암흑물질 입자가 아니다”와 정확히 같은 종류의 경고이며, 이 분야에서 분해능 논쟁이 끝나지 않는 근본 이유이기도 하다. 요컨대 예쁜 은하 그림 속의 점 하나는 별이 아니라 통계량이다. ↩