쿼크-글루온 플라스마

물리 계산물리 유체역학 마지막 수정: 2026-09-05 04:31:47

1. 개요[편집]

쿼크-글루온 플라스마 Quark-Gluon Plasma (QGP)
정체	쿼크와 글루온이 하드론에 갇히지 않고 돌아다니는 QCD 물질 상태
생성 조건	$T \gtrsim 155\ \mathrm{MeV}$ · 에너지밀도 $\sim 1\ \mathrm{GeV/fm^3}$ 이상
만드는 법	상대론적 중이온 충돌 — RHIC(Au+Au, 200 GeV) · LHC(Pb+Pb, 5.02 TeV)
수명 · 크기	$\sim 10\ \mathrm{fm}/c$ ($10^{-23}$ 초) · 지름 $\sim 10\ \mathrm{fm}$
대표 성질	거의 완전 유체 — $\eta/s$ 가 KSS 하한 $1/4\pi$ 의 몇 배 이내
계산 파이프라인	초기조건 → 전평형 → 점성 유체 → 입자화 → 하드론 애프터버너

우주에서 가장 뜨겁고, 가장 작고, 가장 짧게 존재하며, 인류가 아는 가장 잘 흐르는 액체.

쿼크-글루온 플라스마(quark-gluon plasma, QGP)는 온도와 밀도가 충분히 높아져 쿼크와 글루온이 개별 하드론 안에 갇혀 있지 못하고 계 전체에 걸쳐 색 자유도로 돌아다니는 QCD 물질의 상태다. 빅뱅 이후 대략 마이크로초 이내의 우주가 이 상태였고, 지금은 중이온 가속기 안에서 매 충돌마다 $10^{-23}$ 초 동안만 다시 만들어진다.

이 문서는 시뮬레이션 관점에 무게를 둔다. QCD의 상전이 성질과 격자 계산 자체는 격자 QCD가, 강결합 수송계수의 홀로그래피적 해석은 AdS-CFT 대응이 이미 다루므로 여기서는 넘긴다. 여기서 볼 것은 격자에서 나온 상태방정식이 어떻게 유체 코드의 입력이 되고, 검출기의 각도 분포가 어떻게 점성 값으로 환산되는가 하는 파이프라인이다. 요약하면 QGP 연구는 물리학이면서 동시에 **하나의 거대한 역문제**다.

2. 왜 존재하는가 — 점근 자유성과 크로스오버[편집]

양자색역학의 결합상수는 에너지 눈금이 올라갈수록 약해진다. 그로스·윌첵·폴리처(1973)의 점근 자유성이며, 이 성질 때문에 온도를 올리면 유효 결합이 작아져 색 전하가 가려지지 않고 색 차폐(데바이 차폐) 길이가 짧아진다. 결과적으로 어떤 온도 이상에서는 가둠이 풀린다는 것이 1970년대 중반부터 예상되어 있었다.

문제는 그 전이가 어떤 종류인지를 섭동론으로 답할 수 없다는 것이었고, 답은 격자 QCD에서 나왔다. 현재 확립된 그림은 이렇다.

바리온 화학퍼텐셜이 0인 경우( $\mu_B = 0$ ), 실제 쿼크 질량에서의 전이는 1차 상전이도 2차도 아닌 매끄러운 크로스오버다. 아오키 등(2006)의 결과이며, 유한 부피 스케일링으로 판정했다.
유사임계 온도는 카이랄 감수율 등으로 정의하며 $T_{pc} \approx 155\ \mathrm{MeV}$ (약 $1.8\times10^{12}$ K). 관측량마다 값이 몇 MeV씩 다르다 — 크로스오버에는 유일한 전이 온도가 없다는 사실의 직접적 결과다.
유한 $\mu_B$ 에는 QCD 임계점이 존재할 수 있다는 예상이 있으나, 격자는 부호 문제 때문에 그 영역에 직접 들어가지 못한다. RHIC의 빔 에너지 스캔이 요동 관측량(순바리온수 누적률)으로 이 임계점을 찾는 프로그램이다.

크로스오버라는 사실이 시뮬레이션 쪽에 미치는 영향은 실용적이다. 잠열이 없고 상태방정식이 매끄러우므로, 유체 코드가 상변화 근처에서 특별한 처리를 할 필요가 없다. 1차 상전이였다면 혼합상 구간에서 음속이 0이 되는 구간을 다뤄야 했고, 충격파 처리와 수치 안정성이 훨씬 골치 아팠을 것이다.

3. 상태방정식 — 격자에서 유체로 넘어가는 지점[편집]

여기가 이 문서의 핵심이다. 상대론적 유체 방정식을 닫으려면 압력과 에너지밀도의 관계 $p(\varepsilon)$ 가 필요한데, 이 관계를 제1원리에서 주는 것이 격자 QCD다.

격자 쪽 계산의 표준 경로는 소위 적분법이다. 트레이스 어노멀리(상호작용 척도)

\Theta(T) \equiv \frac{\varepsilon - 3p}{T^4}

를 여러 온도에서 재고, 그것을 $T$ 에 대해 적분해 압력을 얻는다.

\frac{p(T)}{T^4} - \frac{p(T_0)}{T^4_0} = \int_{T_0}^{T} \frac{\Theta(T')}{T'}\,dT'

이후 $\varepsilon = 3p + T^4\Theta$ , $s = (\varepsilon+p)/T$ 로 나머지를 채운다. 2+1 맛깔 물리 질량에서 연속 극한까지 밀어붙인 결과가 2014년경 HotQCD와 부다페스트-부퍼탈 두 협력단에서 독립적으로 나왔고, 서로 일치했다. 이 합의가 이후의 모든 유체 시뮬레이션이 딛고 선 바닥이다.

읽는 법에서 자주 틀리는 대목이 둘 있다.

$\Theta$ 는 $T \approx 200\ \mathrm{MeV}$ 부근에서 봉우리를 만들고 고온에서 서서히 0으로 간다. $\Theta = 0$ 은 등각(비질량) 극한이며, 자유 기체라면 $p = \varepsilon/3$ 이 정확하다. QGP는 관측 가능한 온도 범위 어디에서도 자유 기체가 아니다.
슈테판-볼츠만 극한 대비 압력은 도달 가능한 온도에서 대략 80% 수준에 머문다. “가둠이 풀렸다”가 “상호작용이 없다”를 뜻하지 않는다. 이 20%가 QGP를 기체가 아니라 액체로 만든다.

실무에서는 격자 결과( $T \gtrsim 150$ MeV)와 하드론 공명 기체 모형( $T \lesssim 150$ MeV)을 매끄럽게 이어 붙인 표 형태의 상태방정식을 쓴다. 코드에 들어가는 것은 방정식이 아니라 보간표이며, 그래서 이 표의 보간 오차와 미분 가능성(음속 $c_s^2 = dp/d\varepsilon$ 이 튀지 않아야 한다)이 실제 수치 안정성 문제로 직결된다.

4. 유체 코드 — 왜 1차 이론을 못 쓰는가[편집]

QGP는 크기 $\sim10\ \mathrm{fm}$ , 수명 $\sim10\ \mathrm{fm}/c$ 짜리 물체다. 이런 것을 연속체 유체로 다루는 게 정당한지부터가 이 분야의 오래된 논쟁이었고, 결론적으로는 놀랄 만큼 잘 작동한다는 것이 실험 데이터가 준 답이다.

푸는 방정식은 에너지-운동량 텐서의 보존이다.

\partial_\mu T^{\mu\nu} = 0, \qquad T^{\mu\nu} = \varepsilon\, u^\mu u^\nu - (p+\Pi)\Delta^{\mu\nu} + \pi^{\mu\nu}

여기서 $\pi^{\mu\nu}$ 가 전단 응력, $\Pi$ 가 벌크 압력이다. 순진하게 하면 이것을 상대론적 나비에-스토크스(1차 이론)로 닫고 싶어지는데, 그러면 안 된다.

상대론적 1차 점성 이론(에카르트·란다우 형식)은 인과율을 위반하고 평형 해가 불안정하다. 확산항이 정보를 무한 속도로 전파시키기 때문이며, 히스콕-린들럼(1985)이 이 불안정성을 명시적으로 보였다. 수치적으로는 격자 간격을 줄일수록 발산하는 형태로 나타난다.

그래서 표준은 이스라엘-스튜어트 계열의 2차 이론이다. 응력 자체에 완화 시간 $\tau_\pi$ 를 준 전개 방정식을 하나 더 푼다.

\tau_\pi\, \Delta^{\mu\nu}_{\alpha\beta} u^\lambda \partial_\lambda \pi^{\alpha\beta} + \pi^{\mu\nu} = 2\eta\, \sigma^{\mu\nu} + \cdots

이렇게 하면 전파가 유한 속도가 되어 쌍곡형이 되고 수치적으로 안정해진다. 대가는 방정식 수와 미정 계수(2차 수송계수)의 증가다. 최근에는 유체 변수의 정의를 바꿔 1차에서도 안정·인과적인 형식을 만드는 BDNK 계열이 연구되지만, 대규모 데이터 분석의 주력은 여전히 이스라엘-스튜어트다.

수치 도식 자체는 전산유체역학의 표준을 그대로 쓴다. 보존형 유한체적법, 충격 포착이 가능한 고해상도 도식, 원시변수 복원(보존변수에서 $\varepsilon, u^\mu$ 를 뉴턴법으로 되찾는 단계 — 상대론 코드 특유의 병목), 밀러-라그랑주 좌표(고유시간 $\tau$ 와 공간 래피디티 $\eta_s$ )를 쓴 종방향 처리. 수치 상대론 쪽 코드와 기술적 조상이 같다.

5. 타원 흐름 — 각도 분포에서 점성을 읽어내기[편집]

QGP의 점성이 어떻게 측정되는지가 이 분야에서 가장 우아한 대목이다.

두 원자핵이 정면이 아니라 비스듬히 충돌하면, 겹치는 영역이 원이 아니라 아몬드 모양이 된다. 유체라면 압력 기울기가 짧은 축에서 더 크므로 그 방향으로 더 세게 밀어내고, 결과적으로 초기 공간 이방성이 최종 운동량 이방성으로 번역된다. 이 번역을 검출기에서 읽는 방법이 방위각 분포의 푸리에 전개다.

\frac{dN}{d\phi} \propto 1 + 2\sum_{n\ge1} v_n \cos\big(n(\phi - \Psi_n)\big)

$v_2$ (타원 흐름)가 아몬드 모양에 대응하는 주 성분이다.
$v_3$ (삼각 흐름)은 매끄러운 아몬드에서는 나오지 않는다. 그럼에도 관측되며, 그 기원은 핵자 위치의 사건별 요동이다. 2010년경 이 사실이 확립되면서 “매끈한 평균 초기조건” 시대가 끝나고 사건별(event-by-event) 시뮬레이션이 표준이 되었다.

점성이 들어오는 자리가 여기다. 전단 점성은 이방성 변환의 효율을 떨어뜨린다. 유체가 끈적할수록 초기 모양 정보가 운동량 공간으로 덜 전달되므로 $v_n$ 이 작아지고, 그 억제는 고차 조화( $n$ 이 클수록)에서 훨씬 강하다. 그래서 $v_2, v_3, v_4$ 를 동시에 맞추면 $\eta/s$ 가 강하게 제한된다.

결론적으로 나온 값이 유명한 그것이다. RHIC과 LHC 데이터를 점성 유체로 맞추면 $\eta/s$ 가

\frac{\eta}{s} \sim (1\!\sim\!3)\times\frac{1}{4\pi}

수준으로, AdS-CFT 대응에서 나온 KSS 하한 $1/4\pi \approx 0.08$ 의 몇 배 이내다. 물이나 헬륨보다 자릿수로 작은 값이고, 그래서 **“거의 완전 유체”**라는 표현이 붙었다.¹

다만 정직하게 붙일 단서가 있다. 추출되는 것은 상수 하나가 아니라 온도에 의존하는 함수 $\eta/s(T)$ 이고, 벌크 점성 $\zeta/s(T)$ 도 같이 들어간다. 벌크 점성은 등각성이 깨진 정도를 재는 양이라 $T_{pc}$ 근처에서 봉우리를 가질 것으로 예상되며, 이것을 무시하면 $\eta/s$ 추정이 편향된다. 즉 이 분야에서 ” $\eta/s$ 의 값”은 언제나 나머지 모형 가정과 묶인 조건부 진술이다.

6. 젯 소광 — 매질을 관통시켜 보기[편집]

유체가 집단 거동을 보는 창이라면, 젯 소광(jet quenching)은 매질의 미시적 성질을 보는 창이다.

충돌 초기에 큰 운동량 전달로 만들어진 고에너지 파톤은 QGP를 뚫고 나오면서 글루온을 방출해 에너지를 잃는다. 관측 신호는 두 갈래다.

핵 변형 인자 $R_{AA}$ . 중이온 충돌의 고운동량 하드론 수율을 양성자-양성자 충돌 수율에 충돌 횟수를 곱한 값으로 나눈 비. 매질 효과가 없으면 1이어야 하는데, RHIC과 LHC에서 고운동량 영역에 대해 1보다 크게 작은 값이 나온다.
디젯 비대칭. 반대 방향으로 나가는 두 젯 중 매질을 더 길게 통과한 쪽이 에너지를 더 잃어 에너지 균형이 깨진다. LHC 초기 Pb+Pb 데이터에서 눈으로 보일 만큼 뚜렷하게 관측됐다.

여기서 수송계수 $\hat q$ (단위 길이당 횡운동량 제곱의 증가율)를 뽑는데, 유체 쪽 $\eta/s$ 와는 완전히 다른 이론 틀에서 나온 값이라 두 결과의 정합성이 중요한 교차검증이 된다. 실무 계산은 유체 배경장 위에 파톤을 몬테카를로로 흘리는 하이브리드 방식이며, 유체 코드가 준 온도장이 그대로 젯 코드의 입력이 된다.

7. 전체 파이프라인과 베이즈 보정[편집]

지금 이 분야의 표준 계산은 다섯 단계짜리 조립 라인이다. 어느 한 단계도 단독으로 검증되지 않는다는 것이 방법론적 핵심이다.

단계	하는 일	미정 매개변수의 예
초기조건	충돌 직후의 에너지밀도 분포를 사건별로 생성	정규화, 핵자 폭, 요동 세기
전평형	유체 서술이 성립하기 전 구간의 처리	유체화 시작 시각
점성 유체	상태방정식 + 이스라엘-스튜어트 전개	$\eta/s(T)$ , $\zeta/s(T)$ , 완화 시간
입자화	등온 초곡면에서 유체를 하드론으로 변환(쿠퍼-프라이)	전환 온도, 점성 보정 분포
애프터버너	하드론 기체의 산란·붕괴를 수송 코드로	단면적 처방

매개변수가 십수 개이고, 계산 한 사건이 비싸며, 관측량은 수백 개다. 손으로 맞추는 것은 불가능하다. 그래서 2010년대 중반 이후 이 분야가 채택한 것이 베이즈 추론 + 대리 모델 조합이다.

매개변수 공간을 라틴 하이퍼큐브 등으로 성기게 뽑아 수백~수천 점에서 전체 파이프라인을 돌린다.
출력 관측량을 주성분으로 압축하고, 각 성분에 대해 가우스 과정 대리 모델을 학습시킨다.
대리 모델 위에서 마르코프 연쇄 몬테카를로를 돌려 $\eta/s(T)$ , $\zeta/s(T)$ 의 사후분포를 얻는다.

결과로 나오는 것이 “점 추정값”이 아니라 신뢰구간을 가진 온도 의존 곡선이라는 점이 중요하다. 실제로 이 방식이 정착한 뒤 밝혀진 사실 중 하나가, 벌크 점성을 넣지 않은 예전 분석이 $\eta/s$ 를 체계적으로 잘못 잡고 있었다는 것이다. 불확실성 정량화가 중심값을 바꾼 교과서적 사례이며, 이 분야가 검증 및 확인 문화를 진지하게 받아들인 계기이기도 하다.²

8. 남은 논쟁[편집]

작은 계에서의 집단 거동. 양성자-납, 심지어 고다중도 양성자-양성자 충돌에서도 $v_2$ 가 관측된다. 이렇게 작은 계에 유체 서술을 적용하는 게 정당한가, 아니면 초기 상태의 색 상관에서 오는 다른 기원인가가 미결이다. 유체가 예상보다 훨씬 빨리, 훨씬 작은 계에서도 성립한다는 “유체화 어트랙터” 연구가 여기서 나왔다.
초기조건의 축퇴. $v_n$ 은 초기 이방성과 점성의 곱에 민감하다. 초기조건 모형을 바꾸면 $\eta/s$ 가 배수로 움직인다. 사후분포의 폭 상당 부분이 이 축퇴에서 온다.
유한 밀도 영역. 임계점 탐색은 격자가 부호 문제로 막혀 있어 유효 모형에 의존한다. 낮은 충돌 에너지에서는 유체 서술 자체도 흔들려서, 수송 모형과 유체를 섞은 하이브리드가 필요하다.
자기장. 충돌 초기에 지구 자기장의 $10^{18}$ 배급 자기장이 생기고, 이것이 카이랄 자기 효과 같은 이색적 신호를 만들 수 있다는 예측이 있다. 동일 질량수 동중원소 충돌 비교 실험이 수행됐고, 해석은 아직 진행 중이다.

9. 여담[편집]

QGP 발견 발표에서 물리학자들이 가장 놀란 것은 틀린 예측이 맞았기 때문이 아니라, 맞을 줄 알았던 예측이 틀려서였다.³ 1990년대의 표준 기대는 “가둠이 풀리면 약결합 기체가 나온다”였다. 실제로 나온 것은 기체가 아니라 인류가 아는 가장 잘 흐르는 액체였고, 이 반전이 강결합 계산 도구(격자, 홀로그래피)에 대한 수요를 폭발시켰다.
“인공 태양”이나 “블랙홀 만들기” 류의 기사 제목이 이 분야에 자주 붙는데, 실제로 만들어지는 물체는 원자핵 하나 크기에 $10^{-23}$ 초를 사는 것이다. 위험한 것은 물리가 아니라 그 기사를 읽고 계산 자원을 신청하는 사람의 기대치다.
온도를 켈빈으로 환산하면 $\sim2\times10^{12}$ K 이고 태양 중심의 10만 배쯤 된다. 그런데 열용량이 워낙 작아서, 손에 닿아도 화상을 못 입힌다. 온도와 열량을 헷갈리면 안 된다는 열역학 첫 강의의 가장 극단적인 예시.

10. 관련 문서[편집]

격자 QCD · 양자색역학 · 게이지 이론 · 재규격화군
AdS-CFT 대응 · 등각장론 · 부호 문제
전산유체역학 · 유한체적법 · 충격파 · 점성
수치 상대론 · 상전이 · 볼츠만 방정식
마르코프 연쇄 몬테카를로 · 대리 모델 · 불확실성 정량화 · 민감도 해석
검증 및 확인 · 몬테카를로 방법 · 역문제

여기서 “완전 유체”는 일상어가 아니라 $\eta/s$ 가 $\hbar/4\pi k_B$ 근처라는 뜻이다. 무차원화가 $\hbar$ 로 되어 있다는 게 요점 — 이 값이 작다는 건 점성이 양자역학적 하한에 근접했다는 의미다. 참고로 KSS 하한은 정리가 아니라 추측이고, 고계 미분 보정을 넣은 중력 이론에서는 이 값을 밑도는 예가 구성된다. 자세한 단서는 AdS-CFT 대응에 정리되어 있다. ↩
이 분야의 베이즈 분석이 CFD 하는 사람에게 특히 흥미로운 이유가 있다. 우리가 난류 모델 계수를 손으로 튜닝하고 “그래도 검증은 했다”고 말하는 관행과 구조가 똑같은데, 저쪽은 그 튜닝 전체를 사후분포로 만들어 논문에 싣는다. 매개변수 십수 개짜리 파이프라인의 불확실성을 정직하게 보고하는 것이 가능하다는 실물 증거이며, 동시에 그러려면 대리 모델 학습용으로 계산을 수천 번 돌려야 한다는 청구서도 같이 보여 준다. ↩
그래서 2005년 RHIC의 네 실험군이 낸 종합 보고서들에 “이것은 기체가 아니라 액체다”라는 결론이 나란히 실렸을 때, 언론 요약은 “가장 뜨거운 물질”로 갔고 이론가들은 “가장 끈적하지 않은 물질”에 주목했다. 발표에서 조심스럽게 쓰인 표현도 “발견했다”가 아니라 “우리가 만든 물질은 예상했던 약결합 기체가 아니다” 쪽이었다. 이 분야가 자기 발견을 부정문으로 발표한 것은 데이터가 이론의 기본 가정을 깼기 때문이지 겸손해서가 아니다. ↩