일루스트리스에서 TNG까지 — 컴퓨터가 재창조한 우주 거대 구조

1941년, 컴퓨터로 우주를 시뮬레이션한다는 것은 누구도 상상하지 못하던 시대였습니다. 그로부터 반세기 뒤인 2005년, 독일 가르힝(Garching)에 있는 막스 플랑크 천체물리학 연구소의 볼커 스프링겔(Volker Springel)이 이끄는 팀은 수퍼컴퓨터에 빅뱅 직후의 조건을 입력하고, 138억 년의 우주 역사를 재현하는 시뮬레이션을 완성했습니다.

프로젝트의 이름은 밀레니엄 시뮬레이션(Millennium Simulation). 이 시뮬레이션 하나가 천문학계에서 1,000편이 넘는 논문을 낳았습니다.

밀레니엄 시뮬레이션(2005): 첫 번째 우주

밀레니엄 시뮬레이션은 당시 기준으로 전례 없는 규모였습니다.

밀레니엄 시뮬레이션 제원 (2005)
────────────────────────────────────────────────────
항목                    수치
────────────────────────────────────────────────────
시뮬레이션 박스 크기    685 Mpc (한 변)
암흑물질 입자 수        21억 6천만 개
암흑물질 입자 질량      8.6 × 10^8 M☉ (태양 질량)
시간 해상도             최소 수백만 년 단위
저장 스냅샷 수          64개
소요 계산 시간          ~350,000 CPU 시간
────────────────────────────────────────────────────

밀레니엄은 암흑물질만 포함한 시뮬레이션이었습니다. 보통 물질(별, 가스, 블랙홀)은 반경험적(semi-empirical) 모형으로 추가되었습니다. 즉, 컴퓨터가 직접 계산한 것이 아니라, 관측 데이터를 기반으로 한 경험칙으로 보정한 것이었습니다.

그럼에도 밀레니엄이 만들어낸 우주 웹(cosmic web)의 모습은 충격적이었습니다. 필라멘트와 공동, 그리고 필라멘트 교차점의 초은하단이 실제 관측과 놀랍도록 닮아 있었습니다. "우주는 암흑물질의 뼈대 위에 세워진다"는 이론이 처음으로 시각적으로 확인된 순간이었습니다.

스프링겔은 당시 논문에서 이렇게 썼습니다: "우주 구조의 계층적 성장은 이제 더 이상 추론이 아닙니다. 우리는 그것을 컴퓨터 안에서 재현할 수 있습니다."

일루스트리스(Illustris, 2014): 가스와 별을 더하다

밀레니엄의 한계는 분명했습니다. 암흑물질은 중력만으로 움직이지만, 보통 물질은 훨씬 복잡합니다. 가스는 압력과 온도에 반응하고, 별이 폭발해 주변 가스를 날려버리고(초신성 피드백), 블랙홀은 에너지를 방출해 은하 형성을 억제하거나 촉진합니다.

이 복잡성을 처음으로 우주론적 규모에서 구현한 것이 2014년 하버드-스미소니언 천체물리학 센터의 마르크 포겔스베르거(Mark Vogelsberger)와 스프링겔이 이끄는 팀이 발표한 일루스트리스(Illustris) 시뮬레이션이었습니다.

일루스트리스 제원 (2014)
────────────────────────────────────────────────────
항목                    수치
────────────────────────────────────────────────────
시뮬레이션 박스 크기    106.5 Mpc (한 변)
암흑물질 입자 수        18억 2천만 개
가스 셀 (Voronoi 격자)  18억 2천만 개
별/항성풍 입자          ~10억 개
블랙홀 입자             수십만 개
물리 과정               중력, 유체역학, 복사 냉각, 별 형성,
                        초신성 피드백, AGN 피드백
소요 계산 시간          ~1,900만 CPU 시간
────────────────────────────────────────────────────

일루스트리스가 만들어낸 가상 우주에는 실제와 비슷한 타원 은하와 나선 은하가 자연스럽게 형성되었고, 은하들의 색깔 분포(청색-적색 이분 분포, blue-red bimodality)와 크기-광도 관계가 관측과 상당히 일치했습니다.

그러나 일루스트리스에도 중요한 문제점이 있었습니다.

첫째, 시뮬레이션이 예측한 은하군과 은하단 안의 가스 비율이 실제 X선 관측보다 너무 높았습니다.

둘째, 강한 AGN(활성 은하핵) 피드백 모형을 사용했음에도, 적색 타원 은하의 수가 관측보다 너무 많았습니다.

셋째, 은하 크기와 형태의 세부 분포가 관측과 미묘하게 어긋났습니다.

IllustrisTNG(2018): 다음 세대

이러한 한계를 극복하기 위해 포겔스베르거와 스프링겔 팀은 피드백 메커니즘을 전면 재설계했습니다. 특히 **자기장(magnetic field)**과 새로운 블랙홀 피드백 모형을 포함시켰습니다. 그 결과가 2018년 공개된 **IllustrisTNG(The Next Generation)**입니다.

IllustrisTNG는 세 가지 해상도 버전으로 구성되었습니다.

버전	박스 크기	주요 특징
TNG50	50 Mpc	초고해상도, 개별 은하 세부 구조 연구
TNG100	100 Mpc	일루스트리스와 비교 가능한 크기
TNG300	300 Mpc	대규모 통계, 은하단 연구

TNG에서 개선된 주요 결과는 다음과 같습니다.

은하 색깔 이분 분포: 새로운 AGN 피드백 모형 덕분에 청색 별 형성 은하에서 적색 수동 은하로의 전이가 훨씬 현실적으로 구현되었습니다.

성간 자기장: 은하 내부와 은하군 간 공간의 자기장 구조가 Faraday 회전 측정과 잘 일치했습니다.

은하 형태: 나선 은하의 원반과 팽대부, 타원 은하의 광도 분포가 관측과 크게 개선되었습니다.

TNG 데이터는 2019년 완전히 공개되어, 현재 전 세계 수백 개 연구팀이 활용하고 있습니다. 2023년 기준으로 TNG 관련 논문은 1,000편 이상이 출판되었습니다.

FLAMINGO(2023): 1조 셀의 우주

2023년, 영국 더럼 대학교와 레이던 대학교 팀이 새로운 기준을 세웠습니다. **FLAMINGO(Full-hydro Large-scale structure simulations with All the MapINGs)**는 IllustrisTNG보다 수십 배 큰 부피를 동일한 해상도로 시뮬레이션했습니다.

FLAMINGO 제원 (2023)
────────────────────────────────────────────────────
항목                    수치
────────────────────────────────────────────────────
최대 시뮬레이션 박스    2,800 Mpc (한 변) — L2p8_m9 모델
입자/셀 총 수           ~3조 개 (가장 큰 버전)
소요 계산 시간          수억 CPU 시간
주요 물리 과정          중력, SPH 유체역학, 복사 냉각,
                        별 형성, 초신성, AGN, 뉴트리노
보정 방법               기계학습(ML) 에뮬레이터 + 베이지안 교정
────────────────────────────────────────────────────

FLAMINGO의 특별한 점은 **통계적 보정(calibration)**에 있습니다. 이전 시뮬레이션들은 피드백 파라미터를 수동으로 튜닝했습니다. FLAMINGO는 대신 **기계학습 에뮬레이터(ML emulator)**를 활용해 관측 데이터(CMB 렌즈, X선 기체 분율, 은하 질량 함수)와 시뮬레이션을 베이지안 방식으로 자동 정합했습니다.

이 덕분에 FLAMINGO는 우주론적 파라미터의 불확실성을 시뮬레이션 내에서 직접 탐색할 수 있습니다. 예를 들어, 뉴트리노 질량이나 암흑에너지 상태 방정식이 달라지면 우주 웹이 어떻게 달라지는지를 정량적으로 예측합니다.

시뮬레이션 세대 비교

시뮬레이션	연도	물리 과정	박스 크기	핵심 기여
밀레니엄	2005	암흑물질만	685 Mpc	우주 웹 최초 재현
일루스트리스	2014	DM + 가스 + 별 + BH	106.5 Mpc	유체역학 포함 첫 대규모 시뮬
IllustrisTNG	2018	+ 자기장 + 개선 피드백	50~300 Mpc	은하 색깔 이분포 구현
FLAMINGO	2023	+ 뉴트리노 + ML 보정	최대 2,800 Mpc	우주론 파라미터 공간 탐색

시뮬레이션과 현실 사이의 남은 간극

2023년 현재, 우주 시뮬레이션은 놀라운 성공을 거두었지만 아직 해결되지 않은 문제들이 있습니다.

소규모 구조 문제(Small-Scale Structure Problem): ΛCDM 시뮬레이션은 큰 은하 주위에 수백 개의 위성 은하(satellite galaxy)가 있어야 한다고 예측합니다. 하지만 실제로 관측되는 위성 은하의 수는 훨씬 적습니다. 이를 **위성 문제(missing satellites problem)**라고 합니다. 피드백 물리학이 개선되면서 이 차이가 줄었지만, 완전히 해결되지는 않았습니다.

형태-밀도 관계: 실제 관측에서 은하단 중심부에는 타원 은하가 많고 바깥쪽에는 나선 은하가 많습니다. 시뮬레이션은 이 경향성을 대체로 재현하지만, 세부적인 분포는 여전히 차이가 있습니다.

허블 긴장(Hubble Tension): 관측으로 측정한 허블 상수(H₀ ~ 73 km/s/Mpc)와 CMB 기반 이론 예측(H₀ ~ 67 km/s/Mpc)이 맞지 않는 문제입니다. 이 긴장이 시뮬레이션의 구조 형성 예측에도 영향을 줄 수 있습니다.

밝고 큰 은하(BCG)의 크기-광도 관계: 은하단 중심의 가장 밝은 은하(Brightest Cluster Galaxies, BCG)의 크기와 광도가 시뮬레이션의 예측보다 일관되게 큰 경향이 있습니다.

미래의 시뮬레이션: 엑사스케일 컴퓨팅

2030년대에는 엑사스케일 컴퓨터(exascale computer, 초당 10^18 부동소수점 연산) 위에서 돌아가는 다음 세대 시뮬레이션이 등장할 것입니다. 이 시뮬레이션은 현재보다 수십~수백 배 큰 부피를 현재보다 훨씬 높은 해상도로 처리할 수 있게 됩니다.

무엇보다 중요한 것은, 실제 관측과의 대조입니다. 유클리드 망원경(Euclid), 루빈 천문대(Rubin/LSST), SKA(Square Kilometre Array)가 수십억 개 은하의 위치와 형태를 측정하면, 시뮬레이션은 "이렇게 생겼을 것"이 아니라 "정확히 이것이어야 한다"는 기준으로 검증받게 됩니다.

스프링겔은 2022년 강연에서 이렇게 말했습니다: "우리는 처음으로 이론과 관측이 동시에 성숙하는 시대에 살고 있습니다. 관측이 더 정밀해지고, 시뮬레이션이 더 현실적이 될수록, 남은 차이들은 새로운 물리학의 창문이 됩니다."

밀레니엄이 보여준 암흑물질의 뼈대, 일루스트리스가 불어넣은 가스와 별의 생기, TNG가 추가한 자기장과 피드백의 복잡성, FLAMINGO가 열어 보인 통계적 우주론. 컴퓨터 안에서 다시 만들어진 우주는 매 세대 조금씩 더 진짜를 닮아가고 있습니다. 그리고 그 차이가 줄어드는 속도만큼, 우리가 우주에 대해 이해하는 것들이 늘어가고 있습니다.

빅뱅으로부터 138억 년. 그 역사를 컴퓨터 안에서 몇 주 만에 재생하는 시대가 왔습니다. 그 재생이 완벽해지는 날, 우리는 비로소 우주가 왜 지금의 모습인지를 완전히 이해할 수 있을 것입니다.