쌍 유도 이방성으로 보는 무질서 자성 시스템

초록

이 논문은 희석된 자기 반도체 Ga₁₋ₓMnₓN에서 최근접 이웃 Mn–Mn 쌍이 개별 이온의 대칭을 깨뜨려 새로운 단축축 이방성을 만든다는 ‘쌍 유도 이방성(pair anisotropy)’ 개념을 제시한다. DFT 계산으로 쌍에 의한 구조 변형과 스핀‑궤도 결합을 분석하고, 이를 원자 스핀 시뮬레이션에 포함시켰을 때 실험적 자화 곡선과의 일치도가 크게 향상됨을 보여준다.

상세 분석

본 연구는 무질서한 희석 자성 시스템에서 전통적으로 사용되어 온 단일 이온 이방성 모델이 근본적인 한계를 가지고 있음을 지적한다. 특히 Ga₁₋ₓMnₓN과 같은 wurtzite 구조의 희석 자기 반도체에서는 Mn 이온이 12개의 인접 Ga 자리 중 하나를 차지할 확률이 높아, 평균 Mn 농도 x ≈ 6 %일 때 약 52 %의 Mn이 최소 하나의 최근접 이웃 Mn과 쌍을 형성한다. 이러한 쌍은 두 가지 물리적 메커니즘을 통해 새로운 이방성을 유도한다. 첫 번째는 두 번째 Mn이 존재함에 따라 국부 결정장 환경이 변형되어 기존의 삼각형(trigonal) 변형과 Jahn‑Teller(JT) 변형이 재조정되는 것이다. DFT 결과는 단일 Mn이 존재할 때는 삼각형 변형(t ≈ −0.045 Å)과 JT 변형(j₁ ≈ 0.015 Å, j₂ ≈ 0.055 Å)이 동시에 나타나는 반면, Mn–Mn 쌍에서는 JT 변형이 사라지고 Mn–N 결합 길이가 쌍 축 방향으로 비대칭적으로 짧아지는 구조적 비대칭이 주된 원인임을 보여준다. 두 번째 메커니즘은 비등방성 교환( anisotropic exchange)의 대칭적인 부분이 직접적인 스핀 방향 의존성을 부여한다는 점이다. 논문은 비등방성 교환을 별도의 항으로 분리하기보다, 쌍 축을 따라 발생하는 추가적인 단축축(uniaxial) 이방성 K_pair S₁·e_pair S₂·e_pair 형태로 효과적으로 모델링한다.

스핀‑궤도 결합을 포함한 비콜리니어 DFT 계산을 통해 에너지 지형 E(θ, φ)을 얻고, 이를 H = −½K_TR S_z² − K_JT(S·e_JT)² − K_pair(S₁·e_pair)(S₂·e_pair) 형태의 유효 스핀 해밀토니안에 피팅한다. 결과적으로 K_pair은 쌍의 방향에 따라 양(또는 음)값을 가지며, 특히 in‑plane 쌍과 out‑of‑plane 쌍에서 서로 다른 크기의 K_pair이 도출된다. 이러한 파라미터를 원자 스핀 시뮬레이션에 적용했을 때, 실험에서 관측된 저온 자화 곡선(M(H) vs H)과의 차이가 기존 단일 이온 모델 대비 2배 이상 감소한다. 이는 쌍 유도 이방성이 실험적 데이터 재현에 필수적임을 강력히 시사한다.

또한, 논문은 DFT + U나 하이브리드 함수와 같은 전자 상관 보정이 절대 에너지 정확도에는 도움이 되지만, 비콜리니어 스핀‑궤도 계산의 일관성을 해칠 수 있기에 PBE‑GGA를 선택한 이유를 명확히 설명한다. 이는 구조적 변형과 이방성 에너지 차이를 비교적 작은 에너지 스케일(수 meV)에서 정확히 포착하기 위한 전략이다. 마지막으로, 저자는 쌍 유도 이방성이 희석 자성 반도체뿐 아니라 2D van‑der‑Waals 자성체, 스핀 글래스, 클러스터 기반 자성체 등 다양한 무질서 시스템에 보편적으로 적용될 수 있음을 제안한다.