암흑 물질

물리 계산물리 시뮬레이션 마지막 수정: 2026-08-30 04:31:07

1. 개요[편집]

암흑 물질 Dark Matter
정의	중력으로만 존재가 드러나는 비발광 물질
우주 조성	암흑물질 약 26%, 바리온 약 5%, 나머지 암흑에너지
첫 지적	츠비키 (1933, 머리털자리 은하단)
지배 방정식	무충돌 볼츠만 + 푸아송 (블라소프-푸아송)
계산 방법	N체 — 위상공간을 거대입자로 표본추출
표준 프로파일	NFW (내부 r⁻¹, 외부 r⁻³)
미해결	정체. 60년째 후보만 늘어나는 중

은하는 안 보이는 무언가에 붙들려 돈다. 그 무언가를 60년째 못 잡았는데, 그동안 그것으로 우주 구조를 예측했더니 자꾸 맞는다.

암흑 물질(dark matter)은 전자기적으로 거의 상호작용하지 않아 빛으로는 보이지 않지만, 중력 효과로 그 존재가 여러 독립적 관측에서 드러나는 물질이다. 강조점은 “여러 독립적”에 있다. 회전곡선 하나만 있었다면 이 개념은 벌써 폐기됐을 것이다.

심위키에서 이 문서가 하는 일은 두 가지다. 첫째, 증거가 어떻게 쌓여 지금의 무게를 갖게 됐는지를 정리한다. 둘째, 암흑물질이 계산 대상으로서 무엇인지 — 즉 무충돌 볼츠만 방정식의 몬테카를로 표본추출이라는 정체를 분명히 한다. 팽창 좌표·초기조건 생성·바리온 서브그리드 같은 시뮬레이션 운영 문제는 우주론 시뮬레이션이, 중력 솔버 알고리즘(트리·FMM·소프트닝)은 N체 문제와 반스-헛 알고리즘이 이미 다뤘으므로 여기서는 링크로 넘긴다.

2. 증거의 축적[편집]

암흑물질을 믿게 만드는 것은 서로 다른 물리, 서로 다른 스케일, 서로 다른 시대의 관측이 같은 양을 가리킨다는 사실이다.

은하 회전곡선의 평탄화. 원운동 조건 $v^2(r) = GM(<r)/r$ 에서, 빛나는 물질이 거의 다 중심에 있다면 바깥에서는 케플러적으로 $v\propto r^{-1/2}$ 로 떨어져야 한다. 실제로는 광학 원반 바깥까지 $v(r)\approx$ 일정이다. 이는 $M(<r)\propto r$ , 즉 $\rho\propto r^{-2}$ 인 헤일로를 요구한다. 1970년대 베라 루빈과 켄트 포드의 광학 분광, 그리고 원반보다 훨씬 멀리 뻗는 중성수소 21 cm 전파 관측이 이 그림을 굳혔다.

은하단 속도분산. 1933년 츠비키가 머리털자리 은하단 은하들의 속도분산에 비리얼 정리를 적용해, 광도에서 추정한 질량보다 훨씬 큰 질량이 필요하다고 지적했다. 그가 쓴 허블 상수 값이 지금보다 여러 배 컸기 때문에 불일치의 배율 자체는 과대평가였지만, 결론의 방향은 옳았다. 오늘날 은하단 질량은 속도분산·X선 기체 정역학·중력 렌즈 세 가지 독립 경로로 측정되고 서로 일치한다.

총알 성단. 두 은하단이 충돌한 계(1E 0657-558)에서, X선으로 보이는 고온 기체(바리온 질량의 대부분)는 충돌로 감속돼 가운데 남았는데 약중력렌즈로 재구성한 질량 중심은 양쪽 은하 무리를 따라 그대로 지나갔다. 중력의 원천이 보이는 물질과 공간적으로 분리된 것이라, “중력 법칙을 고치면 된다”는 대안에게 가장 아픈 관측이다.

CMB 음향 봉우리. 재결합 이전 광자-바리온 유체의 음향 진동이 남긴 각 파워 스펙트럼의 봉우리 높낮이가 조성을 직접 잰다. 바리온은 유체를 무겁게 만들어 압축 위상(홀수 번째 봉우리)을 선택적으로 키우고, 광자와 결합하지 않는 암흑물질은 퍼텐셜 우물의 깊이를 통해 전체 진폭과 감쇠를 조절한다. 2번째·3번째 봉우리의 상대 높이에서 $\Omega_b h^2 \approx 0.0224$ , $\Omega_c h^2 \approx 0.120$ 이 나온다. 즉 암흑물질은 바리온의 약 5배이며, 이 값은 은하 하나 보지 않고 재결합 시기 물리만으로 얻은 것이다.

구조형성의 타이밍. 이게 의외로 강력한 논증이다. 바리온은 재결합 이전에 광자와 묶여 있어 요동이 자라지 못한다. CMB에서 관측되는 요동 진폭은 $\delta\sim10^{-5}$ 이고 물질 우세기 이후 요동은 대략 $\delta\propto a$ 로 자라므로, 재결합( $z\approx1100$ ) 이후에만 자란다면 오늘 $\delta\sim10^{-2}$ 에 그친다. 은하도 은하단도 없어야 한다. 광자와 결합하지 않아 물질-복사 동등 시점부터 미리 자라기 시작한 성분이 있어야, 바리온이 나중에 그 우물로 굴러떨어져 따라잡을 수 있다.

여기에 하나 더 — 원시 핵합성이 중수소 존재비로 정한 바리온 밀도가 CMB 값과 일치한다. 두 완전히 다른 물리가 같은 $\Omega_b$ 를 주고, 그 값이 총 물질 밀도의 6분의 1밖에 안 된다. 나머지는 바리온이 아니다.

3. ΛCDM 에서의 역할과 후보들[편집]

표준 우주론 모형 ΛCDM 에서 암흑물질에게 요구되는 성질은 놀랍도록 적다.

차갑다 — 요동이 자랄 무렵 속도분산이 작아 자유흐름 길이가 짧다.
무충돌 — 자기들끼리도, 바리온과도 거의 상호작용하지 않는다.
안정 — 우주 나이보다 훨씬 오래 산다.

이 세 가지만 넣으면 CMB부터 은하단 개수, 라이먼-α 숲, 약중력렌즈까지 폭넓게 맞는다. 후보는 이 요구를 어떻게 만족하느냐로 갈린다.

부류	대표 후보	관측적 특징	상태
차가운 (CDM)	WIMP, 액시온, 원시 블랙홀	소규모 구조를 그대로 유지	표준
따뜻한 (WDM)	스테라일 중성미자 (keV급)	자유흐름이 소규모 파워를 절단	라이먼-α 숲이 강하게 제약
자기상호작용 (SIDM)	암흑 부문 매개입자	헤일로 중심을 산란으로 평탄화	총알 성단이 단면적 상한
뜨거운	표준 중성미자	자유흐름 길이가 너무 큼	주성분으로는 배제

뜨거운 암흑물질이 배제된 논리가 계산적으로 재미있다. 상대론적 입자는 자유흐름으로 작은 스케일의 요동을 지워 버리므로, 큰 구조가 먼저 생기고 쪼개지는 하향식 구조형성을 준다. N체 시뮬레이션이 그 시나리오를 돌려 보니 은하 분포가 관측된 것보다 훨씬 덩어리져 나왔고, 1980년대에 이 비교가 뜨거운 암흑물질을 사실상 끝냈다. 시뮬레이션이 모형을 죽인 초기 사례다.

4. 계산 대상으로서의 암흑 물질[편집]

여기서부터가 이 문서의 본론이다. 암흑물질의 동역학은 다음 한 쌍으로 완전히 결정된다.

\frac{\partial f}{\partial t} + \mathbf v\cdot\nabla_x f - \nabla_x\Phi\cdot\nabla_v f = 0, \qquad \nabla^2\Phi = 4\pi G\!\int\! f\,d^3v

위가 무충돌 볼츠만 방정식(블라소프 방정식), 아래가 포아송 방정식이다. 미지수 $f(\mathbf x,\mathbf v,t)$ 는 6차원 위상공간의 분포함수이며, 이것이 문제의 전부이자 문제 그 자체다. 6차원 격자에 한 축당 1000점만 깔아도 $10^{18}$ 셀이다. 직접 이산화는 원리적으로 가능하고 실제로 저차원 블라소프 코드가 있지만, 우주론 규모에서는 불가능하다.

그래서 쓰는 방법이 특성선을 따라가는 몬테카를로 표본추출이다. $f$ 를 유한 개의 거대입자로 표본추출하고, 각 입자를 특성선 방정식 $\dot{\mathbf x}=\mathbf v$ , $\dot{\mathbf v}=-\nabla\Phi$ 로 밀어낸다. 여기서 반드시 못 박아야 할 것.

시뮬레이션 입자는 암흑물질 입자가 아니다. 우주론 시뮬레이션의 입자 하나는 보통 $10^5\!\sim\!10^9\,M_\odot$ 짜리이며, WIMP라면 그 안에 $10^{60}$ 개 이상이 들어 있다. 입자는 위상공간의 표본점이지 물체가 아니다.

이 사실이 코드 설계를 전부 지배한다.

소프트닝은 근사가 아니라 필수다. 표본점끼리의 이체 산란은 물리가 아니라 표본추출 잡음이다. 뉴턴 힘을 $r<\varepsilon$ 에서 완화해 인공적인 이완을 억제해야 무충돌 극한에 수렴한다. 이완 시간은 대략 $t_{\rm relax}\sim \frac{N}{8\ln N}t_{\rm cross}$ 라 $N$ 이 크면 저절로 억제되지만, 헤일로 중심처럼 국소 입자수가 적은 곳에서는 여전히 문제다.
분해능 판정은 이완 기준으로 한다. 밀도 프로파일을 어디까지 믿을지는 “그 반지름 안쪽에서 이체 이완 시간이 우주 나이보다 짧은가”로 정한다. 은하 중심부 논쟁의 상당 부분이 실은 이 판정 기준 싸움이었다.
산탄잡음이 통계량에 직접 들어간다. 입자 표본추출은 유한 개수의 푸아송 잡음을 낳고, 이것이 측정된 파워 스펙트럼에 $1/\bar n$ 만큼 더해진다. 처리법은 파워 스펙트럼 문서 참고.

힘 계산 자체는 $O(N^2)$ 직접합에서 출발해 PM(격자 + 고속 푸리에 변환), 트리, TreePM, 고속 다중극자법으로 진화했고, 시간 적분은 심플렉틱 적분기(킥-드리프트-킥)를 쓴다. 이 부분은 N체 문제가 정리해 두었다.

입자 표본추출의 대안으로 위상공간 시트 방법도 있다. 초기의 3차원 라그랑주 시트를 사면체로 테셀레이션해 그 변형을 추적하면, 밀도를 입자 개수가 아니라 시트의 접힘 횟수로 얻는다. 껍질 교차(caustic) 구조가 선명하게 드러나고 저밀도 영역의 잡음이 극적으로 줄지만, 강한 비선형 영역에서 시트가 미친 듯이 접혀 사면체 수가 폭발하는 것이 한계다.

5. 시뮬레이션이 내놓은 예측들[편집]

ΛCDM 이 검증 가능한 이론인 이유는 N체 계산이 구체적이고 틀릴 수 있는 예측을 주기 때문이다.

NFW 프로파일. 나바로-프렝크-화이트(1996~97)는 질량·형성 시기·우주론을 바꿔 가며 헤일로를 만들어도 구면 평균 밀도 프로파일이 하나의 두 매개변수 꼴로 수렴한다는 것을 발견했다.

\rho(r) = \frac{\rho_s}{(r/r_s)\left(1+r/r_s\right)^2}

안쪽은 $r^{-1}$ , 바깥은 $r^{-3}$ 이고, 자유도는 질량과 집중도 $c=r_{200}/r_s$ 뿐이다. 이 보편성은 해석적으로 유도된 것이 아니라 순수한 수치 실험 결과이며, 지금은 안쪽 기울기가 조금 완만한 아인나스토 꼴이 더 잘 맞는다는 것이 정설에 가깝다. 중력 렌즈 모형화부터 직접 검출 실험의 국소 밀도 추정까지 이 프로파일 위에 얹혀 있다.

질량함수. 단위 부피당 질량 $M$ 헤일로의 개수는 프레스-셰흐터(1974)의 해석적 논증에서 출발한다. 선형 밀도장을 스케일 $R(M)$ 로 평활화한 분산 $\sigma(M)$ 을 구하고, 구형 붕괴의 임계 대비 $\delta_c\approx1.686$ 을 넘는 영역의 비율을 세는 방식이다. 이 이론은 형태는 맞히지만 정량적으로 어긋나고, 셰스-토먼(1999) 이후로는 N체 시뮬레이션에 직접 피팅한 형태를 쓴다. 요점은 여기서도 같다 — 이론이 뼈대를 주고 계수는 시뮬레이션이 준다.

소규모 문제들. 예측이 구체적이면 어긋남도 구체적으로 드러난다.

사라진 위성은하(missing satellites, 1999): 시뮬레이션이 예측한 은하수 주위 부헤일로 수가 알려진 위성은하 수보다 훨씬 많았다. 이후 SDSS·DES 등이 초희미 왜소은하를 대량 발견하고, 작은 헤일로에서는 별 형성이 재이온화로 억제된다는 계산이 붙으면서 상당 부분 해소됐다.
코어-첨점 문제(core-cusp): 시뮬레이션은 $\rho\propto r^{-1}$ 첨점을 주는데, 왜소은하 회전곡선은 평탄한 코어를 선호한다. 후보 설명은 (a) 초신성 되먹임에 의한 반복적 기체 유출이 퍼텐셜을 흔들어 암흑물질을 밖으로 밀어낸다, (b) 암흑물질이 자기상호작용을 한다, (c) 회전곡선 관측 자체의 계통 오차. 바리온 물리를 넣은 최근 시뮬레이션은 (a)만으로도 상당한 코어를 만든다.
너무 커서 실패(too big to fail, 2011): 가장 무거운 부헤일로들이 관측된 어떤 위성은하보다도 조밀해 보이는 문제. 위와 같은 후보들이 경쟁 중.

여기서 정직하게 짚을 것. 이 문제들의 표준 해법이 “바리온 되먹임”이라는 것은, 반박 가능성이 그만큼 약해졌다는 뜻이기도 하다. 서브그리드 처방에 자유도가 많으므로 어지간한 불일치는 흡수된다. 그래서 이 분야는 되먹임이 손댈 수 없는 관측량(고적색편이 왜소은하 개수, 헤일로의 존재 자체, 강한 렌즈의 미세 구조)으로 시험 축을 옮기고 있다.

6. 검출 시도와 MOND[편집]

직접 검출은 지하 검출기에서 암흑물질 입자가 원자핵과 탄성 산란하며 남기는 keV급 반발 에너지를 노린다. 제논 이중상 TPC(XENONnT, LZ, PandaX)가 주도하며, 스핀 무관 산란단면적 상한을 수십 GeV 질량대에서 $10^{-47}\,\mathrm{cm^2}$ 아래까지 밀어 내렸다. 감도 향상에는 원리적 종착점이 있다 — 태양·대기 중성미자의 결맞음 산란이 구별 불가능한 배경이 되는 중성미자 안개다. 여기 도달하면 노출을 늘려도 상한이 거의 개선되지 않는다.

측정값을 상한으로 바꾸는 과정 자체가 통계·계산 문제다. 국소 암흑물질 밀도와 속도분포(표준 헤일로 모형: 맥스웰-볼츠만 + 탈출속도 절단), 핵 형상인자, 검출기 반응 함수를 전부 넣은 우도를 세우고 최대우도추정이나 베이즈 방법으로 신뢰구간을 낸다. 속도분포를 어떻게 가정하느냐가 상한을 배수로 흔들기 때문에, 요즘은 시뮬레이션에서 얻은 헤일로 속도분포를 쓰거나 아예 분포에 무관한 상한을 제시한다.

간접 검출은 쌍소멸·붕괴 산물(감마선, 양전자, 중성미자)을 찾는다. 왜소구상은하는 바리온이 거의 없어 배경이 깨끗해 가장 좋은 표적이고, 은하 중심 감마선 초과는 신호일 수도 있고 밀리초 펄서 무리일 수도 있다. 가속기 쪽에서는 LHC의 단일 제트 + 소실 횡운동량 탐색이 상보적 제약을 준다. 세 경로가 서로 다른 가정에 의존하므로, 한 곳에서 나온 신호는 나머지 둘에서 확인돼야 한다.

MOND(수정 뉴턴 역학, 밀그롬 1983)는 방향을 바꾼다. 어두운 물질을 더하는 대신 가속도가 $a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}$ 보다 작을 때 중력 법칙이 달라진다고 놓는다. 성과를 과소평가하면 안 된다 — 개별 은하 회전곡선을 바리온 분포 하나와 보편 상수 하나로 예측한다. 바리온 툴리-피셔 관계의 좁은 산포와 방사 가속도 관계의 규칙성은 순수 ΛCDM 쪽에서 오히려 설명하기 껄끄러운 정도로 깔끔하다.

문제는 스케일을 넘어갈 때다.

은하단에서 MOND는 여전히 2배가량의 보이지 않는 질량을 요구한다. 암흑물질을 없애려는 이론이 암흑물질을 남긴다.
총알 성단의 질량-기체 분리는 국소적 수정 중력으로 설명하기 어렵다.
CMB 3번째 봉우리와 물질 파워 스펙트럼을 동시에 맞추려면 결국 결합하지 않는 성분이 필요하다.
상대론적 확장(TeVeS 등)은 여러 개가 제안됐지만, 중성자별 병합 GW170817에서 중력파와 빛이 사실상 같은 속도로 도착한 사실이 중력파 속도를 바꾸는 부류를 대거 정리했다.

정직한 요약은 이렇다. MOND는 은하 스케일의 현상론으로서 여전히 설명해야 할 규칙성을 지적하고 있고, 우주론적 이론으로서는 ΛCDM 을 대체하지 못한다.¹

7. 여담[편집]

암흑물질을 지지한 이른 논거 중에 순수 계산에서 나온 것이 있다. 1970년대 초 오스트라이커와 피블스가 차가운 원반 은하를 N체로 돌렸더니 막대 불안정으로 금방 뭉개졌고, 뜨거운 구형 헤일로가 감싸면 안정해진다는 결과가 나왔다. 관측이 아니라 시뮬레이션이 “은하는 헤일로 안에 있어야 한다”고 먼저 말한 셈이다.
“암흑”이라는 이름 때문에 신비한 물질처럼 들리지만, 물리학자 입장에서 이건 가장 지루한 종류의 물질이다. 빛과 안 놀고, 자기들끼리도 안 놀고, 그냥 중력으로 뭉치기만 한다. 지루하기 때문에 시뮬레이션이 잘 되는 것이고, 잘 되기 때문에 정밀 우주론이 성립한다.
원시 블랙홀이 후보로 주기적으로 부활한다. 중력파 관측이 예상보다 무거운 항성질량 블랙홀을 여럿 잡을 때마다 논문이 쏟아지는데, 미시렌즈·CMB 강착 제약 때문에 전체 암흑물질을 설명할 수 있는 질량 창은 상당히 좁다. 자세한 배경은 블랙홀.
검출 실험 쪽 밈 중 최고는 “DAMA/LIBRA”다. 20년 넘게 연간 변조 신호를 보고하는데 다른 어떤 실험도 재현하지 못했고, 같은 요오드화나트륨 결정을 쓰는 후속 실험들이 반박 결과를 내놓는 중이다. 재현성이 무엇인지에 대한 살아 있는 교보재.²

8. 관련 문서[편집]

우주론 시뮬레이션 · N체 문제 · 반스-헛 알고리즘 · 고속 다중극자법
파워 스펙트럼 · 은하 형성 · 중력 렌즈
NFW 프로파일 · 수정 뉴턴 역학
포아송 방정식 · 고속 푸리에 변환 · 심플렉틱 적분기
몬테카를로 방법 · 최대우도추정 · 마르코프 연쇄 몬테카를로
블랙홀 · 계산물리 · 검증 및 확인

이 분야에서 가장 소모적인 논쟁 방식이 “MOND가 회전곡선을 이렇게 잘 맞히는데?” 대 “총알 성단은?”의 무한 반복이다. 양쪽 다 사실이라는 게 함정. 생산적인 질문은 “왜 ΛCDM 안에서 바리온 분포가 회전곡선을 그렇게 잘 결정하는가”이고, 실제로 되먹임이 만드는 바리온-암흑물질 되먹임 고리로 방사 가속도 관계를 재현하려는 시뮬레이션 연구가 진행 중이다. ↩
연간 변조는 지구 공전 때문에 암흑물질 바람에 대한 상대속도가 6월에 최대가 되는 효과라, 신호로서는 매우 그럴듯한 서명이다. 문제는 계절에 따라 변하는 것이 암흑물질만이 아니라는 점 — 온도, 지하수, 라돈, 뮤온 플럭스가 전부 1년 주기를 갖는다. 계통 오차 후보 목록이 신호보다 길면 대개 신호가 아니다. ↩