카멜라 호르모니 계층 방정식과 피크온 솔루션의 새로운 전개

본 논문은 카멜라‑호르모니(CH) 방정식과 그 계층 구조를 스펙트럼 이론과 일반화 푸리에 변환(GFT)이라는 통합적인 프레임워크로 재구성한다. 서두에서는 CH 방정식 u_t−u_{xxt}+2ωu_x+3uu_x−2u_xu_{xx}−uu_{xxx}=0이 물리적으로 얕은 물 파동, 비탄성 원통 등 다양한 현상을 모델링함을 언급하고, 라크스 쌍 Ψ_{xx}= (¼+λ(m+ω))Ψ, Ψ_t= (½λ^{−1}−u)Ψ_x+u_xΨ+γΨ (γ는 임의 상수) 를 제시한다. 여기서 m=u−u_{xx}이며, ω는 실수 상수이다. CH 방정식은 두 개의 호밀토니안 구조 J_1=(2ω∂+m∂+∂m), J_2=∂−∂^3 를 갖는 bi‑Hamiltonian 시스템이며, 무한 개의 보존량 H_n이 존재한다. 두 번째 장에서는 재귀 연산자 L=J_2^{-1}J_1을 도입하고, 그 고유함수가 바로 라크스 스펙트럼 문제의 제곱 고유함수 F^{±}(x,k)=(f^{±}(x,k))^2 임을 증명한다. f^{±}(x,k)는 각각 x→±∞에서 평면파와 일치하도록 정규화된 해이며, k는 복소 파라미터로 λ(k)=−(k^2+¼)/ω 로 정의된다. 제곱 고유함수는 스위치츠 급수 공간에서 완전성을 만족하므로, 임의의 실함수 q(x)=m(x)+ω 를 q(x)−ω = ±(1/2π i)∫_{ℝ} (R^{±}(k)/λ(k)) F^{±}(x,k) dk + Σ_{n=1}^N (κ_n/λ_n) R_n^{±} F_n^{±}(x) 와 같이 전개할 수 있다. 여기서 R^{±}(k)와 R_n^{±}는 연속·이산 스캐터링 데이터이며, κ_n은 이산 스펙트럼의 고유값이다. 변분을 추가 파라미터(예: 시간 t, 보조 변수 y)와 함께 취하면, 스캐터링 데이터는 선형 진화 방정식 R_t ∓ i k Ω(λ^{-1})R=0, R_{n,t} ± κ_n Ω(λ_n^{-1})R_n=0, λ_n,t=0 로 기술된다. 세 번째 장에서는 다항식 P_1(z), P_2(z) 로 정의된 Ω(z)=P_1(z)P_2(z) 를 이용해 계층 전반의 비선형 방정식 P_2(L) ∂^{-1}_±(√q)_t √q + P_1(L) (ω/q−1)=0 을 L과 ∂^{-1} 연산자를 통해 q_t + 2 q \tilde u_x + q_x \tilde u = 0, \tilde u = ½ Ω(L)(ω/q−1) 의 형태로 변환한다. Ω(z)=a_{−1}z^{−1}+a_0+a_1z 를 선택하면 구체적인 비선형 방정식 (23)이 도출되며, 계수 선택에 따라 원 CH 방정식, Dym 방정식, Dullin‑Gottwald‑Holm 방정식 등 다양한 알려진 모델이 재현된다. 이는 재귀 연산자와 GFT가 계층 전체를 하나의 연산자 Ω(L) 로 포괄한다는 강력한 결과이다. 해밀토니안 구조는 포아송 괄호 {A,B}=−∫(ω+m)(δA/δm)∂(δB/δm)−(δB/δm)∂(δA/δm)dx 로 정의하고, 스캐터링 데이터에 대한 해밀토니안 H_Ω = ∫_ℝ k^2 Ω(λ^{-1}) π^{-1} ω λ^{-2} ln(1−R^+(k)R^−(−k)) dk − 2 ω Σ_{n=1}^N κ_n^2 λ_n^{-2} Ω(λ_n^{-1}) dκ_n 을 제시한다. Ω(z)=z^{-1} 일 때는 H_Ω가 CH의 두 번째 해밀토니안 H_{CH,2}와 첫 번째 해밀토니안 H_{CH,1}의 선형 결합으로 나타나며, Ω(z)=z^{-2} 일 때는 (−1) 차 해밀토니안 H_{CH,−1}이 등장한다. 이는 스캐터링 데이터와 물리적 보존량 사이의 직접적인 연결고리를 제공한다. 네 번째 장에서는 (1+2) 차원 확장식 q_t + 2(U_{xy}+η U_{xx}) q + (U_y+η U_x+γ) q_x = 0, q = U_x − U_{xxx} + ω 을 고려한다. 여기서 η,γ는 임의 상수이며, 라크스 쌍은 비스펙트럼 형태로 변형된다. λ는 λ_t − (1/2)λ λ_y = 0 를 만족하도록 가정하고, 스캐터링 데이터가 t와 y에 동시에 의존함을 보인다. 결과적으로 R_t − (1/2)λ R_y ± 2 i k (γ−η^2/(2λ)) R = 0 와 같은 선형 방정식이 얻어지며, γ=η=0 인 경우 R(k,t,y)=R(k,t+2λ y) 형태의 이동파 해가 된다. ω→0 한계에서 제곱 고유함수는 δ-함수 형태의 피크온을 생성하고, q(x)=∑_{k=1}^N p_k(t,y) δ(x−x_k(t,y)) 라는 ansatz 로 N-피크온 해를 구성한다. 연산자 L을 이용해 x_k와 p_k에 대한 비선형 시스템 (36)·(37) 을 도출하고, N=1,2 경우에 대해 명시적인 해를 제시한다. 특히 N=2 경우에 p_1+p_2가 보존되는 등 기존 CH 피크온의 보존량이 그대로 유지됨을 확인한다. 다만, y-파라미터가 등장함에 따라 해밀토니안 구조가 퇴화하고, Ω(z)≡0 인 경우는 해밀토니안이 정의되지 않음에 주의한다. 마지막으로 저자는 스위치츠 급수 공간에서 제곱 고유함수의 완전성을 이용한 GFT가 CH 계층뿐 아니라 그 변형, 고차원 확장식까지 일관된 해석적 틀을 제공한다는 점을 강조한다. 스캐터링 데이터의 선형 진화와 재귀 연산자에 의한 비선형 방정식의 재구성은 기존의 역산술 변환 방법을 크게 일반화하며, 새로운 통합적 접근법으로서 향후 다른 비선형 파동 방정식에도 적용 가능함을 시사한다.