약한 무질서 격자에서 q‑브레이어의 비선형 특성과 제어 가능성

본 연구는 비선형과 무질서가 동시에 존재하는 격자 진동계의 핵심 현상을 이해하고, 이를 제어할 수 있는 방법을 제시한다. 먼저, 저자들은 고전적인 FPU‑β 체인을 변형하여 스프링 상수에 무작위 변동 κₙ을 도입한 모델을 정의한다. 이 모델은 N개의 동일 질량 입자를 갖고, 경계 조건 x₀=x_{N+1}=0을 만족한다. 해밀턴식(1)은 운동 에너지와 포텐셜 에너지(선형·비선형 부분)로 구성되며, 무질서 파라미터 D와 비선형 파라미터 β가 핵심 변수이다. 정규모드 변환을 수행하면 q‑공간에서의 정상모드 좌표 Q_q(t)와 고유진동수 ω_q가 도출된다. 무질서가 없는 경우, q‑브레이어(QB)는 특정 중심 모드 q₀에 에너지가 집중된 시간 주기적 해이며, 비선형 상호작용에 의해 q₀, 3q₀, 5q₀…와 같은 고조파 모드에 지수적으로 감소하는 에너지 분포를 보인다. 저자들은 이러한 QB 해를 무질서가 있는 경우에도 연속시킬 수 있음을 보인다. 연속 과정은 두 단계로 이루어진다. 첫 번째 단계에서는 무비선형(β≠0, D=0) QB 해를 정확히 알고 있다고 가정하고, 특정 κₙ 실현을 고정한 뒤 D를 작은 값으로 증가시킨다. ν와 d를 동시에 작은 파라미터로 두고 섭동 전개를 수행한다. 0차 항은 기존 QB 해이며, 1차 무질서 교정은 강제 진동자 형태의 선형 방정식(4)으로 나타난다. 이 식으로부터 모든 모드가 중심 모드와의 주파수 차이에 비례해 작은 진폭 A^{(1)}_q를 갖게 된다. 평균적으로는 ⟨A^{(1)}_q⟩=0이지만, 분산은 κₙ의 통계적 변동 σ_κ²에 비례한다. 그 결과, 모드 에너지 E_q는 비선형 기여 E_NL_q와 무질서 기여 E_DO_q로 분리된다. 식(5)에서 E_DO_q는 q≫q₀일 때 상수값 d²E_{q₀}σ_κ²/2을 갖으며, 이는 Fig.1에 나타난 플래토 형태와 일치한다. 반면 q≈q₀±1에서는 ω_q와 ω_{q₀}의 차이가 작아져 E_DO_q가 크게 증가하고, 이는 QB의 핵심 국소화가 약화되는 원인이다. 저자들은 q_c≈ln(Dσ_κ²(N+1))/lnλ+1이라는 전이점(플래토와 지수적 감소 사이)을 도출하고, 무질서가 충분히 강하면 모든 QB가 탈국소화될 수 있음을 보인다. 안정성 분석에서는 QB 주변을 선형화하고 심플렉틱 플로우 행렬의 고유값 θ_i를 계산한다. |θ_i|=1이면 안정이며, 무질서가 없을 때의 임계 비선형성 β*_0=π²/(6E_{q₀}(N+1))이 기준이 된다. 무질서가 도입되면 평균 β*는 변하지 않지만, 표준편차 σ_β*≈2σ_κ D√(N+1)/E_{q₀}가 D에 비례해 증가한다(식(8)). Fig.2(a)는 두 개의 서로 다른 κₙ 실현에 대해 β*가 증가하거나 감소하는 모습을 보여준다. 특히 저자들은 κₙ을 주기적 함수 κₙ=0.5 cos