Platform for zero-field isolated skyrmions: 4$d$/Co atomic bilayers on Re(0001)

기술 요약: Re(0001) 상의 4d/Co 원자 이중층에서 제로 필드 고립 스카이르미온 플랫폼

문제 정의
자기적 스카이르미온은 위상적으로 보호되는 스핀 텍스처로서, 그 나노 규모의 크기와 안정성 덕분에 차세대 스핀트로닉스 소자, 뉴로모픽 컴퓨팅 및 양자 기술의 유망한 후보로 꼽힌다. 다양한 전이 금속 초박막 및 다층 구조에서 고립된 스카이르미온이 실험적으로 구현되어 왔으나, 이들은 일반적으로 안정화를 위해 외부 자기장을 필요로 한다. 장치를 구성하는 자기장 생성 부품을 제거하고 에너지 소비를 줄이기 위해서는 제로 필드(zero-field) 안정화가 매우 중요하다. 현재까지 전이 금속 초박막 계열 내에서는 교환 좌절(exchange frustration)에 의해 안정화되는 Rh/Co 원자 이중층이 Ir(111) 상에서 10 nm 미만의 제로 필드 고립 스카이르미온을 수용하는 것으로 관찰된 유일한 사례이다. 이론적 예측에 따르면 4d/Co 이중층이 Re(0001) 표면에서도 이러한 상태를 가질 수 있음을 시사하지만, 기존 연구들은 고차 다체 스핀 상호작용(HOI)을 간과한 유효 스핀 모델에 의존하였으며 안정성을 확인하기 위한 원자론적 스핀 시뮬레이션이 부족했다.

방법론
본 연구는 Re(0001) 표면 위에 성장한 세 가지 특정 시스템인 fcc-Rh/Co, hcp-Pd/Co, fcc-Ru/Co 원자 이중층을 조사하기 위해 제일원리 밀도 범함수 이론(DFT)과 원자론적 스핀 시뮬레이션을 결합한 접근 방식을 채택하였다.

1. DFT 계산: 전자 및 자기 구조는 전밀도 선형화 평면파 증강 평면파(FP-LAPW) 방법을 기반으로 하는 fleur 코드를 사용하여 계산되었다. 본 연구는 9개 층의 Re(0001) 상에 놓인 4d/Co 이중층을 모델링하였다.

   • 에너지 분산: 균일 스핀 나선(spin spiral)의 에너지 분산을 스핀-궤도 결합(SOC) 유무에 따라 고대칭 방향($\Gamma K M$ 및 $\Gamma M$)을 따라 전파 벡터 $q$의 함수로 계산하였다.
   • 상호작용 추출: DFT 에너지를 하이젠베르크(Heisenberg) 및 DMI 해밀토니안에 매핑하여 쌍별 교환($J_{ij}$) 및 자일로진스키-모리야 상호작용(DMI, $D_{ij}$) 상수를 추출하였다.
   • 고차 상호작용 (HOI): 표준 모델을 넘어서기 위해, 본 연구는 세 가지 HOI 상수인 이중 사각형(biquadratic, $B_1$), 3-사이트 4-스핀($Y_1$), 4-사이트 4-스핀($K_1$) 상수를 계산하였다. 이들은 특정 다중 Q 상태(uudd 및 3Q 상태)와 그에 대응하는 단일 Q 스핀 나선 상태 사이의 에너지 차이로부터 도출되었다.
   • 결정 자기 이방성 (MAE): 쉬운 자화축을 결정하기 위해 자가 일관적(self-consistently)으로 계산되었다.

2. 원자론적 스핀 시뮬레이션: 모든 DFT 유도 파라미터(하이젠베르크 교환, DMI, MAE 및 HOI)를 포함하는 확장된 원자론적 스핀 모델을 구축하였다.

   • 이완(Relaxation): 제로 온도 자기 상도표를 결정하기 위해 란다우-리프비츠(Landau-Lifshitz) 방정식을 사용하여 $150 \times 150$ 육방 격자에서 스핀 구조를 이완시켰다.
   • 안정성 분석: 지오데식 누지드 엘라스틱 밴드(GNEB) 방법을 사용하여 고립된 스카이르미온이 강자성(FM) 기저 상태로 붕괴되는 최소 에너지 경로(MEP)를 계산하여, 소멸 및 생성 에너지 장벽을 결정하였다.

주요 기여 및 결과

• 제로 필드 스카이르미온 플랫폼 식별: 확장된 스핀 모델 시뮬레이션 결과, fcc-Rh/Co/Re(0001) 및 **hcp-Pd/Co/Re(0001)**에서 고립된 스카이르미온이 제로 자기장에서 강자성 배경으로부터 자발적으로 출현함을 밝혔다. 반면, **fcc-Ru/Co/Re(0001)**에서는 불충분한 DMI 강도로 인해 제로 또는 유한한 자기장에서 고립된 스카이르미온을 안정화할 수 없었다.
• 스카이르미온 크기: 제로 필드 고립 스카이르미온의 반경은 fcc-Rh/Co/Re(0001)의 경우 약 6 nm, hcp-Pd/Co/Re(0001)의 경우 약 12 nm이다. 이 반경은 외부 자기장이 증가함에 따라 감소한다.
• 고차 상호작용의 역할: 본 연구는 근접 이웃 교환 및 DMI 상호작용을 넘어선 상호작용이 이 박막들에서 좌절을 유도한다는 점을 강조한다. 결정적으로, HOI(특히 3-사이트 4-스핀 상호작용)를 포함하는 것이 시스템의 정확한 묘사에 필수적이지만, HOI를 생략한 유효 스핀 모델 또한 기저 상태를 올바르게 예측하였다.
• 에너지 장벽 및 안정성: 전이 상태 이론 계산에 따르면, 안정적인 시스템의 고립된 스카이르미온은 제로 자기장에서 약 150 meV(소멸 장벽)의 상당한 에너지 장벽에 의해 보호된다.
  • 장벽 구성: DMI는 에너지 장벽에 대해 지배적인 양(+)의 기여를 하며, 3-사이트 4-스핀 상호작용은 작은 양의 기여를 한다. 쌍별 교환 및 MAE 항은 장벽을 감소시키는 방향(강자성을 선호하는 방향)으로 작용한다. 4-사이트 4-스핀 상호작용은 무시할 만한 기여를 한다.
  • 관측 가능성: 이러한 장벽의 크기는 이 스카이르미온들이 저온 실험에서 관측될 수 있음을 시사한다.
• 재료 특성:
  • 세 박막 모두 면외(out-of-plane) 쉬운 자화축을 가진다.
  • fcc-Rh/Co 및 hcp-Pd/Co에서의 DMI는 주로 Re 기판에 의해 구동되는 반면, fcc-Ru/Co의 경우 Ru와 Co 층이 서로 다르게 기여하여 총 DMI가 현저히 낮다.
  • Co와 4d 층의 자기 모멘트는 스핀 나선 파수(wave vector)에 따라 크게 변하는데, 이는 고정된 모멘트를 가진 단순 하이젠베르크 모델로는 포착되지 않는 특징이지만, 강자성 상태 근처에서는 모멘트가 상당히 일정하게 유지된다.

의의 및 주장
저자들은 4d/Co 이중층을 Re(0001) 표면 위의 새로운 유망한 플랫폼으로서 나노 규모의 제로 필드 고립 스카이르미온을 실현할 수 있다고 제안한다. 본 연구는 교환 좌절, DMI, 그리고 고차 스핀 상호작용의 복잡한 상호작용에 의해 외부 자기장 없이도 이러한 위상적 텍스처를 안정화할 수 있는 특정 4d 전이 금속(Rh, Pd)과 Co의 조합을 입증하였다.

저자들은 계산된 에너지 장벽이 저온에서 스핀 편극 주사 터널링 현미경(SP-STM)을 통해 실험적으로 관측하기에 충분히 높다고 주장한다. 또한, Re 기판이 1.7 K 미만에서 초전도성을 띠기 때문에, 이 시스템이 창발적 위상 초전도성과 관련될 수 있는 모델 자성-초전도 하이브리드 시스템 역할을 할 수 있음을 언급하였다. 본 연구는 유효 스핀 모델이 이러한 상태의 존재를 예측했을지라도, 스카이르미온의 안정성과 특성을 정확하게 평가하기 위해서는 HOI를 포함한 확장 모델이 필수적임을 강조한다.