Relativistic magnetohydrodynamics in the early Universe

기술 요약: 초기 우주의 상대론적 자기유체역학(Relativistic Magnetohydrodynamics)

문제 제기
초기 우주의 원시 자기장, 유체 섭동, 그리고 중력파 생성에 관한 연구는 흔히 자기유체역학(MHD) 시뮬레이션에 의존한다. 역사적으로 이러한 시뮬레이션은 아상대론적 벌크 유체 운동($u^2 \ll 1$)을 가정해 왔으며, 로런츠 인자 $\gamma^2 \approx 1$인 극한을 빈번하게 가정해 왔다. 이러한 근사는 물질 지배 시대의 유체에는 유효하지만, 복사 지배 시대나 상태 방정식 파라미터 $c_s^2$가 1에 가까운 경우(예: $c_s^2 = 1/3$)에는 검증되지 않은 오차를 유발한다. 기존의 공식들은 로런츠 인자의 시간 미분($\partial_\tau \gamma^2$)을 무시하는 경우가 많았으며, 이는 $U^2/\ell$($U$와 $\ell$은 각각 특성 속도와 길이 척도) 차수의 비선형 항들을 누락시키는 결과를 초래했다. 이러한 누락은 에너지 및 운동량 보존, 와도(vorticity) 생성, 그리고 MHD 파동의 전파 속도에 영향을 미치며, 적절한 교정이 이루어지지 않을 경우 초광속 알펜 속도(Alfvén speeds)를 초래할 수 있다.

방법론
저자들은 균질하고 등방적인 팽창하는 프리드만-로메트르-로버트슨-워커(FLRW) 배경에서의 MHD 보존 법칙을 도출한다. 본 연구는 다음과 같은 몇 가지 이론적 단계를 거친다:

1. 기하학적 프레임워크: 공형 시간(conformal time) $\tau$와 공변 좌표를 사용하여 FLRW 메트릭을 설정하며, 유체 변수의 다양한 스케일링 특성을 탐구하기 위해 일반적인 $\alpha$-시간 정의를 사용한다.
2. 완전 유체 역학: 일정한 상태 방정식 $p = c_s^2 \rho$를 갖는 완전 유체의 운동 방정식을 도출한다. 저자들은 이를 두 가지 형태로 제시한다:
   • 보존 형태: 응력-에너지 텐서 성분 $\tilde{T}^{0\mu}$를 진화시킨다.
   • 비보존 형태: 원시 변수(comoving energy density $\tilde{\rho}$ 및 특이 속도 $u$)를 진화시킨다.
     결정적으로, 이 유도 과정에서는 아상대론적 극한에서 흔히 버려지는 $\partial_\tau \gamma^2$ 및 $\partial_\tau u^2$를 포함하는 항들을 유지한다.
3. 불완전 유체: 고전적 비가역 열역학(Classical Irreversible Thermodynamics, CIT) 접근법을 사용하여 1차 불완전 유체 효과(점성과 열전도)를 포함하도록 프레임워크를 확장하며, 나비에-스토크스 점성과 푸리에 열전도에 대한 공변적 표현을 제공한다.
4. 전자기학: 팽창하는 배경에서의 맥스웰 방정식을 검토하고, 공변적 일반화 옴의 법칙(generalized Ohm's law)과 함께 공변적 전자기장($\tilde{E}, \tilde{B}$)을 도입한다. 변위 전류(displacement current)를 분석하여 초기 우주의 고전도도 영역에서 이상적인 MHD 극한을 정당화한다.
5. MHD 결합: 로런트 힘을 통해 유체 방정식과 맥스웰 방정식을 결합한다. 저자들은 음파, 알펜 파, 그리고 자기음파(magnetosonic waves)를 연구하기 위해 선형 섭동을 분석한다.
6. 교정: 변위 전류를 무시할 때 알펜 속도가 빛의 속도보다 빨라지는 것을 방지하기 위해, 상대론적 MHD 맥락에 맞게 조정된 "보리스 교정(Boris correction)"을 적용한다.

주요 기여 및 결과

• 교정된 아상대론적 방정식: 본 논문은 기존의 아상대론적 MHD 방정식(예: Pencil Code에서 사용되는 방정식)이 대류 항 $(u \cdot \nabla) \ln \tilde{\rho}$ 및 속도 항의 계수에서 오류를 포함하고 있음을 보여준다. 구체적으로, $\gamma^2$의 시간 미분은 복사 지배 유체($c_s^2 = 1/3$)에 대해 특정 항에서 $1/2$의 교정 인자를 도입하며, 일반적인 $c_s^2$에 대해 계수를 수정한다. 교정된 비보존 형태 방정식은 식 (1.2)와 (6.31)에 제시되어 있다.
• 와도 생성: 저자들은 초기 우주에서 와도($\omega = \nabla \times u$)가 더 이상 위상적 불변량이 아님을 보여준다. 외부 강제력이나 바로클리닉(baroclinic) 항이 없는 경우에도, $c_s^2 \neq 0$ 및 우주의 팽창과 관련된 상대론적 교정으로 인해 초기 회전 없는(curl-free) 속도장에서 와도가 생성될 수 있다.
• 스케일링 및 공형 불변성: 유체 변수의 다양한 스케일링(comoving, super-comoving)을 검토하여 허블 마찰 항을 최소화한다. 저자들은 $c_s^2 \neq 1/3$인 아상대론적 흐름의 경우, 특정 스케일링($\beta = 3(1+c_s^2)$)을 통해 에너지 방정식은 공형 평탄(conformally flat)하게 만들 수 있는 반면, 운동량 방정식은 허블 마찰 항을 유지한다는 것을 확인하였다.
• 상대론적 알펜 속도 및 보리스 교정: 이상적인 MHD에서 변위 전류를 무시하면 $B_0^2/\rho_0$가 클 때 초광속 알펜 속도($v_A > 1$)가 발생할 수 있음을 분석 결과가 보여준다. 저자들은 알펜 속도가 $v_A \to v_A \gamma_A \leq 1$ (여기서 $\gamma_A = (1+v_A^2)^{-1}$)이 되도록 운동량 방정식을 수정하여 보리스 교정을 상대론적 MHD에 적응시켰다.
• 수송 계수: 원시 플라즈마에서의 전단 점성, 벌크 점성 및 열전도 추정치를 검토하였으며, 대규모 척도에서 플라즈마가 효과적으로 완전 유체임을 확인하였다. 다만, 현실적인 난류 모델링을 위해서는 소산 효과(dissipative effects)가 필요하다.

의의
본 연구는 이전의 아상대론적 극한에서 간과되었던 교정 사항들을 포함하여, 팽창하는 배경에서의 MHD 방정식을 기술하는 최초의 완전한 상대론적 설명을 제공한다. 저자들은 이러한 교정들이 단순히 고차 상대론적 효과가 아니라, 비선형 유체 역학의 대류 미분과 동일한 차원의 효과라고 주장한다. 결과적으로, 본 논문은 초기 우주의 자기유체역학, 특히 복사 지배 유체나 고에너지 단계와 관련된 시뮬레이션이 비선형 역학, 난류 및 와도 생성을 정확하게 모델링하기 위해 반드시 이러한 교정된 항들을 포함해야 한다고 강조한다. 도출된 방정식들은 Pencil Code나 CosmoLattice와 같은 차세대 수치 코드 개발을 위한 이론적 토대가 된다.