자기장 구동 전진 영역에서 도메인벽의 결정론적 거칠어짐

초록

본 연구는 영온도에서 균일한 외부 자기장에 의해 구동되는 평면형 키랄 도메인벽이 Walker 붕괴장 위에서 선형 불안정성을 보이며, 길이가 일정 임계값을 초과하면 파동 모드가 성장해 결정론적 거칠어짐과 스페이스‑타임 혼돈을 일으킨다는 것을 수치와 이론으로 입증한다. 불안정 영역은 길이와 감쇠계수 α에만 의존하는 준보편적 스펙트럼 안정도 다이어그램으로 요약된다.

상세 분석

논문은 두 가지 접근법을 결합한다. 첫 번째는 마이크로스코픽 LLG(랜드우‑리프시츠‑길버트) 시뮬레이션을 이용해 실제 재료 파라미터(Co/Pt 다층)와 동일한 조건에서 도메인벽(DW)의 전진을 직접 계산한다. 두 번째는 u‑ϕ 집합좌표 모델, 즉 DW의 위치 u(x,t)와 내부 각도 ϕ(x,t)를 변수로 하는 축소된 연속 방정식(식 5)을 도입해 분석적·수치적 검증을 수행한다. 비차원화 과정을 통해 길이 L₀, 시간 T₀, 자기장 H_W 등 자연 단위를 정의하고, Walker 붕괴장 H_W를 기준으로 비차원 자기장 h를 도입하였다.

선형 안정성 분석에서는 평면 DW 해에 대한 작은 변동을 푸리에 변환하여 각 파수 κ에 대해 2×2 선형 시스템을 얻는다(식 15‑16). h<1(정상 영역)에서는 행렬 고유값이 모두 음의 실수부를 가져 안정함을 확인한다. h≥1(전진 영역)에서는 ϕ₁(t)이 시간에 따라 회전하므로 비자율 시스템이 되며, ϕ₁을 새로운 시간 변수로 삼아 Floquet 이론을 적용한다. 한 주기(π) 동안의 전이 행렬 M_κ를 수치 적분해 구하고, 그 고유값 μ_κ(플루켓 승수)를 계산한다. |μ_κ|>1이면 해당 모드가 기하학적으로 불안정함을 의미한다. 결과적으로 불안정 영역은 α에만 의존하는 “불안정 깃털” 형태로 나타나며, 최대 불안정 파수 κ_m와 임계 자기장 h_c는 α에 대해 κ_m∝α^{-1/2}, h_c∝α^{-2}와 같은 스케일을 보인다.

유한 크기 효과도 분석하였다. 시스템 길이가 L/L₀ < 2π/κ_m이면 모든 비영 파수가 안정하므로 강체 DW 움직임이 유지된다. 반대로 L가 충분히 크면 가장 낮은 파수 모드가 불안정해져 파동이 성장하고, 이후 비선형 상호작용에 의해 모드 간 에너지 카스케이드가 발생한다.

비선형 단계에서는 u‑ϕ 방정식의 sin(2ϕ) 항이 핵심 역할을 하며, 불안정 모드가 서로 결합해 스페이스‑타임 혼돈, 복잡한 Bloch‑line 핵 생성·소멸, 그리고 결정론적 거칠어짐을 유발한다. 시뮬레이션 결과는 특정 자기장 h_c에서 평면에서 거친 이동상으로의 동적 상전이가 일어나며, 이 전이점에서 전이 전후의 통계적 지표(예: 전방 속도 평균, 표면 거칠기 지수)가 급격히 변한다. 이러한 현상은 열이나 잡음이 없는 순수 결정론적 시스템에서도 보편적인 비평형 패턴 형성의 한 예로 해석될 수 있다.

마지막으로, LLG 시뮬레이션과 u‑ϕ 모델 사이의 정량적 차이를 보정하기 위해 실제 DW 폭 Δ′와 모델 폭 Δ의 차이를 스케일링(Δ′/Δ)으로 보정하였다. 이 과정을 거친 후 두 방법은 불안정 임계선, 성장 속도, 그리고 비선형 거칠어짐 패턴에서 뛰어난 일치를 보였다. 따라서 제시된 quasi‑universal stability diagram은 Gilbert 감쇠 α만 알면 다양한 재료와 기하학적 조건에 적용 가능한 실용적인 도구가 된다.