로렌츠 부스트 확산의 초기값 문제와 정확 해법

본 논문은 “로렌츠 부스트된 확산 방정식”이 전통적인 형태(∂ₜn=∂ₓ²n)에서는 부스트 프레임에서 지수적 불안정성을 보여 초기값 문제가 ill‑posed가 된다는 오래된 문제를 재조명한다. 기존의 해결책으로는 (1) 원래 방정식의 해를 구한 뒤 로렌츠 변환을 적용하는 방법, (2) 방정식 자체를 부스트하고 Sobolev 공간에서 초기값을 설정하는 방법, (3) Cattaneo와 같이 2차 시간 미분을 도입해 인과성을 강제하는 방법 등이 있다. 그러나 각각은 물리적 제한(예: t>0만 정의, 고주파 불안정, 비공변성)이나 수학적 비현실성을 동반한다. 최근 연구에서 포커‑플랑크 동역학이 상대론적 입자들의 확산을 정확히 기술한다는 사실을 이용해, 저자들은 확산 방정식이 포커‑플랑크 방정식의 수소동역학 부문에 정확히 포함된다는 점을 강조한다. 포커‑플랑크 방정식은 (∂ₜ+V∂ₓ)f = β∂ₚ²f + …, 이며 여기서 V는 입자 속도, β는 외부 매질의 역온도이다. 해를 복소 파수 k와 주파수 ω에 대해 f∝e^{i(kx−ωt)} 형태로 가정하면, 동역학적 일치 조건은 ω=−ik²가 된다. 이는 바로 확산 방정식의 분산관계와 일치한다. 하지만 물리적으로 허용 가능한 해는 적분 ∫dp e^{−βε(p)} e^{−ikβp}이 수렴해야 하는 제약을 받는다. 비상대론적 경우는 Gaussian이므로 모든 k에 대해 수렴하지만, 상대론적 경우 ε(p)∼|p|이므로 수렴을 위해 |Im k|<1이 필요하다(식 (7)). 이 제한은 단순히 공변적 안정성(Im ω≤|Im k|)보다 강력하며, 실제로 포커‑플랑크 이론이 허용할 수 있는 파동수 대역을 정의한다. 이제 로렌츠 부스트( t̃=γ(t+vx), x̃=γ(x+vt) )를 적용하면 파수와 주파수는 k = γ( k̃ − v ω̃ ), ω = γ( ω̃ − v k̃ ). 부스트된 파동수 k̃에 대해 원래의 수렴 조건을 대입하면 |Im ω̃| < 1/(γv) (식 (9))가 얻어진다. 이는 부스트된 확산 방정식(∂_{t̃}+v∂_{x̃})n = γ^{-3}∂_{x̃}²n)의 불안정한 해(예: e^{t̃/(γv²)})를 자동으로 배제한다. 또한, ω̃ = −i k̃²를 부스트 관계에 대입해 |k̃|에 대한 상한 Λ를 구하면 Λ = 1 + 2v/(1−v) (식 (10)) 가 된다. v→0 혹은 v→1일 때 Λ→∞이며, v=¼에서 최소값 2√3≈3.46을 갖는다. 이러한 파동수 제한을 만족하는 초기 데이터는 Paley‑Wiener 공간 PW_Λ에 속한다. 구체적으로, 초기 프로파일은 δn(0, x̃) = ∫_{−Λ}^{Λ} (dk̃/2π) φ(k̃) e^{ik̃x̃}, φ∈L²(