spin density·초기 자기모멘트·자기배열을 탐색해 magnetic ground state를 찾기

DFT · PART 5 · 숫자의 정확도와 이론의 적용 범위를 시험하는 동안03 / 08

Fe, Co, Ni와 전이금속 산화물은 spin을 무시하면 구조·에너지·전자상태가 크게 달라집니다. 그러나 초기 자기모멘트에 따라 서로 다른 SCF 해가 나올 수 있습니다.

어떤 자기상태가 안정한지 어떻게 탐색하고, 계산된 자기모멘트를 어디까지 해석할까요?

이번 장의 답spin DFT는 $n_\uparrow,n_\downarrow$ 또는 총밀도와 자화밀도를 기본 변수로 풉니다. 가능한 FM·AFM·ferrimagnetic 배열과 초기 모멘트를 비교하고 구조까지 완화한 뒤, SOC·비공선성·온도 효과가 필요한지 판단합니다.

01spin-polarized Kohn-Sham DFT

collinear spin 계산에서는 up과 down 전자 밀도를 따로 두고 spin별 유효 퍼텐셜을 풉니다. 두 밀도의 차이를 적분하면 셀의 spin magnetic moment를 얻습니다.

$$ M_s=\mu_B\int[n_\uparrow(\mathbf r)-n_\downarrow(\mathbf r)]d\mathbf r $$
원자별 moment는 공간 분할법에 의존하지만 총 spin moment는 더 명확합니다.

고립 원자와 분자의 spin multiplicity, 전이금속의 d 점유, 결함의 홀전자처럼 초기 spin을 물리적으로 설정해야 합니다. 비자성 초기값은 대칭 때문에 자성 해로 자발적으로 깨지지 않을 수 있습니다.

02자기배열을 구조 후보처럼 탐색하기

FM은 모멘트가 같은 방향, AFM은 이웃 모멘트가 반대 방향, ferrimagnetic은 반대 방향이지만 크기가 달라 순모멘트가 남습니다. 특정 AFM 파장은 원시셀보다 큰 magnetic supercell을 요구합니다.

자기 ground-state 탐색
  • 결정 대칭이 허용하는 대표 magnetic order 생성
  • 원자별 초기 모멘트의 크기와 부호를 여러 값으로 설정
  • 각 상태의 SCF·구조를 같은 기준으로 완화
  • 최종 moment·DOS·대칭과 에너지를 함께 비교
  • 더 복잡한 noncollinear/spiral 상태 가능성 평가

03spin-orbit coupling과 자기 이방성

spin-orbit coupling은 전자의 spin과 궤도 운동을 결합해 자화 방향에 따른 에너지 차이, orbital moment, band splitting을 만듭니다. 무거운 원소와 자기 이방성·topological band에서는 필수입니다.

magnetic anisotropy energy는 meV 이하로 작아 매우 촘촘한 k-grid와 엄격한 SCF, 동일 구조·설정이 필요합니다. collinear scalar-relativistic 계산의 에너지 차이보다 훨씬 민감한 문제입니다.

04이 글의 용어 사전

본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.

spin polarization
up과 down 전자 밀도가 달라 순 자화밀도가 생기는 상태입니다.
자화밀도magnetization density
위치별 spin-up과 spin-down 밀도의 차이입니다.
FMferromagnetic
주요 국소 자기모멘트가 같은 방향으로 정렬된 상태입니다.
AFMantiferromagnetic
이웃 모멘트가 반대 방향으로 정렬되어 전체가 상쇄될 수 있는 상태입니다.
noncollinear magnetism
모든 모멘트가 하나의 축의 양·음 방향으로만 정렬되지 않는 자성입니다.
spin-orbit coupling
전자 spin과 궤도 운동 사이의 상대론적 결합입니다.

05참고자료와 더 읽을거리

  1. Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
  2. Dudarev et al., Electron-energy-loss spectra and structural stability of NiO (1998) — 회전 불변 LSDA+U와 유효 매개변수 Ueff의 대표적 정식화입니다.
  3. Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996) — 재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.
NEXT DFT NOTE밴드갭 문제와 Kohn-Sham 고유값의 한계 →