Local magnon modes studied by dynamic magnetic pair-density function analysis

기술 요약: 동적 자기 쌍 밀도 함수 분석을 통한 국소 마그논 모드 연구

문제 제기
비탄성 중성자 산란(INS)은 역격자 공간에서 동적 자기 구조 인자 $S_M(Q, E)$를 통해 자기 들뜸(마그논)을 연구하는 오랜 역사를 가지고 있지만, 이 방법은 실공간에서의 국소 스핀 역학을 조사하는 데에는 널리 활용되지 못했다. 액체 및 비정질 물질의 국소 원자 구조를 연구하기 위해 PDF(쌍 분포 함수)의 푸리에 변환이 표준적으로 사용되고, 동적 PDF가 국소 포논 역학에 적용되어 온 것과 달리, 그 자기적 대응물인 동적 자기 쌍 밀도 함수(DymPDF), 즉 $D_M(r, E)$는 아직 충분히 탐구되지 않았다. $D_M(r, E)$의 물리적 의미, 특히 실공간에서의 아쿠스틱(acoustic) 모드에서 옵티컬(optical) 모드로의 전이와 관련된 국소 마그논 모드에 대한 이해는 아직 미흡하다. 또한, 유한한 에너지 전달이 있는 비편광 INS 측정에서 자기 성분을 원자 격자 성분으로부터 분리하는 것은 기술적인 과제로 남아 있다.

방법론
본 연구는 비편광 비탄성 중성자 산란을 사용하여 측정된 동적 자기 구조 인자 $S_M(Q, E)$의 푸리에 변환을 통해 얻은 동적 자기 쌍 밀도 함수(DymPDF) 분석법을 채택한다.

• 실험 설정: J-PARC의 4SEASONS 초퍼 분광기를 사용하여 FeTiO$_3$ (면내 강자성체) 및 YBa$_2$Cu$_3$O$_6$ (2D 사각 격자 반강자성체) 다결정 시료에 대한 측정을 수행하였다.
• 데이터 처리: 자기 신호를 분리하기 위해, 유한한 에너지 전달에서의 $Q$ 의존성 차이법을 사용하여 포논 성분을 제거하였다. $S_M(Q, E)$ 데이터는 보스 인자(Bose factor) $\langle n(E) + 1 \rangle$로 나누어졌으며, 최대 운동량 전달($Q_{max}$) 근처의 오차를 억제하기 위해 창 함수(window function)를 사용하여 실공간으로 변환되었다.
• 이론적 프레임워크: 저자들은 단순 하이젠베르크 모델 자성체(1D, 2D, 3D 강자성체 및 반강자성체)에 대한 저에너지 한계에서의 $D_M(r, E)$에 대한 해석적 방정식을 유도하였다. 이 유도는 근접 이웃 상호작용을 갖는 스핀 해밀토니안을 활용하였으며, 국소 모드에 집중하기 위해 주기적 경계 조건 없이 동적 구조 인자를 근사하였다.
• 시뮬레이션: 다양한 격자 차원 및 자기 질서에 따른 $S_M(Q, E)$의 분말 패턴과 그에 대응하는 $D_M(r, E)$를 검증하기 위해 SpinW 소프트웨어를 사용하여 시뮬레이션을 수행하였다.

주요 기여 및 결과

1. $D_M(r, E)$ 방정식 유도: 저자들은 1D, 2D, 3D에서의 단순 강자성 및 반강자성 모델에 대한 DymPDF의 구체적인 방정식을 유도하였다. 이 방정식들은 신호를 마그논 상태 밀도(MDOS), 마그논 위상 인자 $F(r, E)$, 그리고 배위수 항 $C(r)$로 분리한다.
2. 마그논 모드 전이 특성화: 분석 결과, $D_M(r, E)$는 특정 에너지에서 스핀 쌍 상관관계의 특징적인 부호 변화를 나타낸다. 이 부호 변화는 마그논 모드가 아쿠스틱에서 옵티컬로 전이되는 것을 의미한다.
   • 강자성체의 경우, 이 전이는 중간 에너지 영역(예: 시뮬레이션 상에서 1D $\sim 2.4$ meV, 2D $\sim 4.4$ meV, 3D $\sim 6.4$ meV)에서 관찰된다.
   • 반강자성체의 경우, 부호 변화는 주로 마그논 밴드의 상단부 근처에서 관찰된다.
3. 차원성 및 위상 변화: 본 연구는 차원에 따른 마그논 위상 인자의 체계적인 경향성을 식별하였다. 강자성체의 경우, 1D에서 2D로 이동할 때 $-\pi/4$의 위상 변화가 관찰되며, 1D 기준점 대비 3D에서는 $-\pi/2$의 변화가 나타난다. 또한, 배위수를 포함할 때 반경 방향 강도의 스케일링은 1D, 2D, 3D 시스템에 대해 각각 $1$, $\sqrt{r}$, $r$로 나타난다.
4. 실험적 검증:
   • FeTiO$_3$: 실험적인 $D_M(r, E)$는 10.4 meV의 모드 전이 에너지까지 강자성 스핀 쌍 상관관계를 성공적으로 재현하였다. 이 에너지에서의 부호 변화는 아쿠스틱-옵티컬 전이를 확인시켜 주었다. 유도된 교환 파라미터는 시뮬레이션에 기반한 보정 계수를 적용한 후, 기존의 단결정 연구 결과와 일치하였다.
   • YBa$_2$Cu$_3$O$_6$: 넓은 자기 들뜸 스펙트럼(최대 300 meV)과 운동학적 제한에도 불구하고, 분석을 통해 최근접 이웃 Cu 결합 길이(3.9 Å)에서의 반강자성 상관관계를 검출하였다. 결과는 이 양자 스핀 시스템($S=1/2$)에서도 반강자성 상관관계가 존재함을 확인하였다.

의의
본 논문은 DymPDF 분석을 비주기적 조건이나 클러스터화된 물질에서 국소 자기 역학을 연구하기 위한 독특한 도구로 확립한다.

• 국소 역학: 이는 스핀 쌍 상관관계와 그 에너지 의존성을 시각화하는 방법을 제공하며, 역격자 공간 분석($S_M(Q, E)$)이 쉽게 제공하지 못하는 통찰력을 제공한다.
• 모드 식별: 실공간 상관관계의 부호 변화를 통해 아쿠스틱에서 옵티컬 마그논 모드로의 전이를 식별할 수 있는 능력은 새로운 발견이다.
• 복잡한 시스템을 위한 토대: 저자들은 이 방법이 국소 상관관계가 중요한 초전도 스핀 공명 모드와 같은 이색적인 동적 자기 클러스터 상태를 연구하기 위한 기초적인 토대를 형성한다고 주장한다.
• 한계 및 범위: 저자들은 현재의 저에너지 근사와 고전적 계산(양자 스핀 연속체를 무시함)이 보정 계수 없이 복잡한 실제 자성체의 파라미터를 추정하는 데 있어 한계가 있음을 겸허히 언급하였다. 그러나 정성적인 특징의 재현과 FeTiO$_3$ 및 YBa$_2$Cu$_3$O$_6$에 대한 성공적인 적용은 국소 자기 역학을 조사하는 데 있어 이 방법의 유용성을 입증한다.