강한 외부 자기장에서도 유지되는 방사형 스키르미온 안정화

초록

본 논문은 2차원 자성 시스템에 스키르미온 밀도 $q$의 제곱($q^{2}$) 항을 도입해, 역전 대칭을 보존하면서도 강한 외부 자기장 하에서 스키르미온을 안정화시키는 모델을 제시한다. Landau‑Lifshitz‑Gilbert 방정식을 이용한 수치 시뮬레이션과 스케일링 분석을 통해 $q^{2}$ 항이 Zeeman 에너지와 경쟁하여 스키르미온의 크기를 결정함을 확인하고, 작은 방사형 퍼터베이션에 대해 에너지 최소구조가 위상적으로 보호됨을 증명한다.

상세 분석

이 연구는 기존 스키르미온 모델이 비대칭 Dzyaloshinskii‑Moriya(DM) 상호작용에 의존하는 한계를 극복하고자, 스키르미온 밀도 $q=\mathbf{n}\cdot(\partial_x\mathbf{n}\times\partial_y\mathbf{n})/4\pi$ 의 제곱에 비례하는 새로운 항 $H_q\propto q^{2}$ 를 Hamiltonian에 추가한다. $q^{2}$ 항은 3차원 Hopfion 모델에서 등장하는 Skyrme 항의 2차원 버전으로, Manton이 제시한 변형 텐서 $D=JJ^{T}$ 의 두 번째 불변량 $\det D$ 와 동일하다. 중요한 점은 $q^{2}$ 항이 삼각형 형태의 세 스핀을 필요로 하여 2차원 격자에서만 의미를 갖고, 전통적인 Heisenberg 교환($\partial_i\mathbf{n}\cdot\partial_i\mathbf{n}$)과 달리 스케일링 변환 $a\rightarrow Ka$ 에 대해 $K^{-2}$ 로 스케일링된다. 따라서 강한 외부 자기장($H_Z$)이 지배적인 영역에서는 $H_q$와 $H_Z$ 사이의 경쟁이 스키르미온의 안정성을 결정한다.

저자들은 먼저 연속적인 스케일링 분석을 수행해, 전체 에너지 $H_S=H_Z+H_H+H_{DM}+H_q$ 의 임계 스케일 $K_S$ 를 정의하고, $H_Z$ 가 충분히 클 때 $K_S$ 가 $H_q$와 $H_Z$ 의 비율에만 의존함을 보였다(그림 2a). 이때 $H_H$와 $H_{DM}$ 은 무시해도 된다. 이러한 스케일링 결과는 Derrick‑Hobart 정리와 일치하며, 2차원에서 $q^{2}$ 항만으로도 위상적 안정성을 확보할 수 있음을 의미한다.

이후 격자 모델을 도입해 $q$ 를 삼각형 셀에 대한 면적 보존 형태인 arctangent 식(식 2)으로 이산화하고, LLG 방정식(식 4)을 이용해 동역학적 이완 과정을 시뮬레이션한다. Gilbert 감쇠항을 포함함으로써 장시간 진화 후 정적 해에 도달하도록 설계했으며, 정적 조건식(식 6)에서 $q^{2}$ 항이 Zeeman 항에 의해 발생하는 유효 장과 균형을 이루는 것을 확인한다. 중심부에서는 $\mathbf{n}\parallel -\hat z$, 외곽부에서는 $\mathbf{n}\parallel +\hat z$ 로 전이하면서 $q=0$ 선을 경계로 형성한다.

연속극한에서는 Belavin‑Polyakov(BP) Ansatz를 변형해 $\lambda(r)$ 를 반경에 따라 변하도록 두고, 방사형 대칭을 유지한 채 Euler‑Lagrange 방정식을 풀었다. 결과적으로 라디얼 프로파일은 $r_s$ (즉, $\lambda(r_s)=r_s$) 로 정의되는 유한한 크기를 갖고, $q^{2}$ 항의 강도 $\Lambda$ 가 클수록 $r_s$ 가 감소한다(그림 4). 또한 헬리시티 $\gamma$ 에는 자유도가 남아 있어 Neel형($\gamma=0$)과 Bloch형($\gamma=\pm\pi/2$) 스키르미온이 동등하게 존재한다는 점을 확인했다.

마지막으로 에너지 하한을 분석해 전체 Hamiltonian 이 진공 에너지보다 항상 높으며, 위상 전하 $Q$ 가 정수값을 유지함을 보였다. 이는 작은 방사형 퍼터베이션에 대해서도 에너지 변동이 2차 이상이므로, 위상적으로 보호된 안정 상태임을 수학적으로 증명한 것이다.

핵심 기여는 (1) 역전 대칭을 보존하면서도 $q^{2}$ 항만으로 스키르미온을 안정화시키는 새로운 메커니즘을 제시, (2) 강한 외부 자기장 하에서 전통적인 DM‑exchange 의 필요성을 없애고, (3) 스케일링 및 수치 시뮬레이션을 통해 위상적 보호와 에너지 하한을 명확히 입증한 점이다. 이 모델은 DM‑interaction 이 억제된 재료나 고자기장 환경에서 스키르미온 기반 스핀트로닉스 소자를 설계하는 데 직접적인 활용 가능성을 제공한다.