비상호작용 스위프트 헨든버그 모델의 패턴 전이와 차동 구조 분석

초록

본 연구는 1차원 비상호작용(Nonreciprocal) 스위프트‑헨든버그(NRSH) 방정식을 대상으로, 파라미터 ε와 α에 따라 나타나는 다섯 가지 시공간 패턴(무질서, 정렬, 스와프, 키랄‑스와프, 키랄)을 정량적으로 구분한다. 스펙트럼 기반 분류와 주요 1차 푸리에 모드만을 이용한 저차원 축소 모델을 구축하고, 선형·비선형 안정성 분석을 통해 Turing, 파동, 피치포크, 호프 분기선을 도출한다. 이론적 경계선은 직접 수치 시뮬레이션 결과와 일치한다.

상세 분석

논문은 먼저 비상호작용 스위프트‑헨든버그(NRSH) 모델을 두 개의 실수 순서 매개변수 ϕ(x,t), ψ(x,t) 로 정의하고, 비상호작용 계수 α와 상호작용 계수 χ를 도입한다. α≠0인 경우 자유에너지 함수로부터 유도될 수 없으며, 비평형 화학 포텐셜에 기인한 비대칭 상호작용을 의미한다. 모델은 ϕ와 ψ가 동시에 반전 대칭을 갖고, α와 χ의 부호 변환에 따라 변수 교환 대칭이 존재하므로, 연구는 χ≥0, α≥0 영역만을 다룬다.

시공간 패턴은 수치 시뮬레이션(L=2π, 주기적 경계조건)으로 5가지 유형을 확인한다. (1) D(무질서) – 전혀 구조가 나타나지 않음, (2) A(정렬) – 고정된 파동 형태, (3) S(스와프) – 진폭이 시간에 따라 진동하는 정지파, (4) CS(키랄‑스와프) – 진폭 진동과 동시에 한쪽 방향으로 전파, (5) C(키랄) – 일정 진폭을 유지하며 일정 속도로 전파하는 파동.

패턴 구분은 ϕ와 ψ의 시공간 푸리에 스펙트럼을 이용한다. k₀=1인 경우 n=±1 모드가 가장 불안정하며, 스펙트럼 피크의 위치(ω)와 대칭성에 따라 각 패턴을 정의한다. α가 작을 때는 ω=0 피크 하나만 존재해 A상태, α≈1.12에서 피크가 양쪽으로 분리되면 S상태(ω₊=−ω₋, 진폭 동등) 혹은 CS상태(ω₊+ω₋≠0)로 전이한다. α가 더 커지면 ω>0 피크 하나만 남아 C상태가 된다. 이러한 기준을 정량화한 뒤, ε–α 평면에 5가지 영역을 도식화하였다.

이론적 분석은 먼저 선형 안정성 검토로 시작한다. ϕ,ψ를 푸리에 전개하고 n=±1 모드만 고려하면 고유값 λₙ^{±}=ε−(1−n²)²±√(χ²−α²)이다. α≤χ이면 실수 고유값만 존재해 Turing 불안정(λ⁺>0) 조건은 √(χ²−α²)=−ε, 이는 빨간선으로 표시된다. α>χ이면 실수부가 ε이므로 ε>0이면 파동 불안정이 발생하고, 고유값에 허수부 ±i√(α²−χ²) 가 나타나 파동 bifurcation(초록선)으로 해석된다.

비선형 효과를 포착하기 위해 n=±1 모드만 남겨 ϕ₁, ψ₁의 복소 진폭 방정식(11)-(12)를 얻고, 위상 자유도를 제거해 ρ₁, ρ₂, δ(=θ₂−θ₁) 로 구성된 3차원 ODE 시스템(13)-(15)을 만든다. 여기서 δ의 고정점 조건 sinδ=0 은 정렬(A) 혹은 스와프(S, CS) 상태와 연결되고, (χ−α)ρ₁/ρ₂+(χ+α)ρ₂/ρ₁=0 은 키랄(C) 상태와 연결된다.

A상태의 고정점은 ρ₁,ρ₂가 양수이고 δ=0 혹은 π인 경우이며, 이를 만족하는 ρ₁,ρ₂는 4차 방정식(19)으로 정의된다. 방정식의 판별식 Δ(20)의 부호 변화가 고정점 수의 변화를 초래하고, Δ=0이 피치포크(주황선)와 호프(파란선) bifurcation을 구분한다. 특히 Δ=0은 두 고정점이 합쳐지는 saddle‑node(검정선)과 연관된다.

키랄(C) 고정점은 ρ₁∝ε³α/(α+χ), ρ₂∝ε³α/(α−χ)와 cosδ=−χε/α√(α²−χ²) 로 명시된다. 존재 조건 ε>0, α>χ, P≡α⁴−α²χ²−χ²ε²≥0 가 만족될 때만 실존한다. P=0이 바로 피치포크 경계이며, 이 선을 넘어가면 정렬(A)에서 키랄(C)으로 연속적인 전이가 일어난다.

마지막으로, 축소된 3차원 시스템의 수치 해석을 통해 고정점의 안정성, 주기해(호프) 발생, 그리고 두 고정점 사이의 연속적인 전이(피치포크)를 확인하였다. 이론적 경계선은 전체 1D PDE 시뮬레이션에서 관찰된 패턴 전이와 매우 높은 일치를 보이며, 비상호작용이 패턴 형성에 미치는 역할을 명확히 드러낸다.