Spin and Orbital Rashba response in ferroelectric polarized PtSe$_2$/MoSe$_2$/LiNbO$_3$ heterostructures

기술 요약: 강유전성 편극을 가진 $\text{PtSe}_2/\text{MoSe}_2/\text{LiNbO}_3$ 이종 구조에서의 스핀 및 궤도 라쉬바 응답

문제 제기
최근 스핀트로닉스 테라헤르츠(THz) 방출 분야의 발전은 페로마그네틱/비자성(FM/NM) 이종 구조에서 초고속 광펄스가 스핀-전하 변환(SCC)을 통해 THz 복사를 생성할 수 있음을 입증했다. 역 스핀 홀 효과(ISHE) 및 역 라쉬바-에델슈타인 효과(IREE)와 같은 메커니즘은 잘 확립되어 있으나, 특히 반데르발스(vdW) 이종 구조 내에서 궤도 자유도의 역할은 문헌에서 여전히 미개척 영역으로 남아 있다. 또한, $\text{LiNbO}_3$와 같은 강유전성 기판이 편극 의존적 전기장을 통해 THz 방출을 조절할 수 있음이 밝혀졌음에도 불구하고, 이러한 전기장이 스핀 및 궤도 라쉬바 텍스처와 그에 따른 SCC 효율에 어떻게 영향을 미치는지에 대한 제일원리 기반의 정량화 연구는 부족한 실정이다. 본 연구는 $\text{PtSe}_2/\text{MoSe}_2$ 계면에서의 라쉬바 효과를 정량화하고, 실험적 THz 신호와 강유전성 편극에 의한 방출 조절을 설명하기 위해 궤도 역학을 포함하는 필요성을 다룬다.

방법론
저자들은 밀도 범함수 이론(DFT)과 선형 응답 이론을 결합한 다단계 제일원리 접근 방식을 채택하였다:

1. 구조 완화 및 전기장 추출: VASP 코드(평면파 기저 함수, 500 eV 컷오프, D2 vdW 보정) 내의 DFT-PBE를 사용하여 $\text{PtSe}_2/\text{MoSe}_2/\text{LiNbO}_3$ 이종 구조를 모델링하였다. 저자들은 반대 방향의 편극($\hat{P}_{\uparrow}$ 및 $\hat{P}_{\downarrow}$)을 나타내기 위해 Z-컷 $\text{LiNbO}_3$ 슬래브의 상단과 하단에 2D 층을 배치한 두 가지 구성을 시뮬레이션하였다. 전체 이종 구조의 정전기적 퍼텐셜을 분리된 $\text{PtSe}_2/\text{MoSe}_2$ 이중층과 비교함으로써, 2D 계면에 생성되는 유효 수직 전기장($E_{eff}$)을 추출하였다.
2. 수송 및 응답 계산: SIESTA 코드(국소화된 DZP 기저 함수 세트)와 쿠보(Kubo) 형식론을 활용하여 비평형 수송 계산을 수행하였다. 추출된 $E_{eff}$ 값을 모사하기 위해 순수 $\text{PtSe}_2/\text{MoSe}_2$ 삼중층 시스템에 인공적인 전기장을 가하였다.
3. 선형 응답 분석: 저자들은 선형 응답 이론을 사용하여 인가된 면내 전기장($E_x$)에 의해 유도된 스핀 및 궤도 축적(밀도)을 계산하였다. 이는 비평형 스핀($\langle \delta \hat{s}_y \rangle$) 및 궤도($\langle \delta \hat{l}_y \rangle$) 밀도를 결정하기 위해 $200 \times 200$ k-포인트 메쉬에 대해 쿠보 공식을 적분하는 과정을 포함한다. 연구는 실험적 광 여기와 관련된 에너지 창(페르미 준위 대비 약 -2 eV ~ +2 eV)에 집중하였다.

주요 기여 및 결과

• 유효 전기장의 정량화: 본 연구는 $\text{LiNbO}_3$ 기판에 의해 2D 층에 작용하는 유효 전기장을 성공적으로 정량화하였다. $\text{MoSe}_2$가 $\text{LiNbO}_3$ 표면에 더 가까운 구성의 경우 전기장은 약 0.115 V/Å로 계산되었으며, $\text{PtSe}_2$가 더 가까운 구성에서는 ~0.047 V/Å로 계산되었다. 이 값들은 2D 시스템의 편극 의존적 조절을 시뮬레이션하는 데 사용되었다.
• 궤도 라쉬바 텍스처: 저자들은 $\text{PtSe}_2/\text{MoSe}_2$ 시스템의 궤도 텍스처를 매핑하였다. $\Gamma$ 점 근처의 카이랄리티(chirality)는 상대적으로 안정적으로 유지되는 반면, 편극 방향에 따라 $K$ 및 $K'$ 점 근처에서 상당한 왜곡이 발생함을 관찰하였다. 특히, 궤도 텍스처는 스핀 텍스처에 비해 더 육각형에 가까운 형태를 보였으며, $K/K'$에서는 반시계 방향, $\Gamma$에서는 시계 방향의 카이랄리티를 나타냈다.
• 스핀 vs 궤도 축적:
  • 스핀: 시뮬레이션 결과, 편극 방향이 계면을 향하는($\hat{P}_{\uparrow}$) 구간(약 -1.75 ~ -1.0 eV)에서 스핀 축적이 강화됨을 확인하였다. 이러한 강화는 THz 전기장 진폭의 ~30% 증가라는 실험적 관측 결과와 일치한다.
  • 궤도: 궤도 축적은 [-1, 0] eV 창 내에서 스핀 축적과 반대 부호를 나타냈다. 이러한 대립 관계는 에너지 창 전반에 걸친 음의 스핀-궤도 상관관계($\langle \hat{l} \cdot \hat{s} \rangle$)에 기인한다.
• 조절 및 부호 반전: 연구는 축적 부호가 반전되는 특정 에너지 레벨을 식별하였다. 스핀의 경우 0.9 eV 부근에서 부호 변화가 일어나며, 궤도 밀도의 경우 -1.5 eV 부근에서 발생한다. 이는 THz 신호의 부호를 여기 펄스의 광자 에너지나 게이팅을 통한 화학적 퍼텐셜 조절을 통해 이론적으로 튜닝할 수 있음을 시사한다.
• 라쉬바 결합 강도: $\text{PtSe}_2/\text{MoSe}_2$ 계면의 추정된 라쉬바 결합 상수($\alpha_R$)는 약 -300 meV·Å이며, 이는 주로 계면의 Se-p 상태에 의해 구동된다. 이는 $\text{PtSe}_2$/그래핀 시스템에서 발견된 값(-195 meV·Å)보다 높다.
• THz 방출에 대한 궤도 기여: 저자들은 편극 상태 간의 역 라쉬바-에델슈타인 길이($\Lambda_{REE}$)의 백분율 이득을 계산하였다. 스핀 기여가 <30%의 이득을 보인 반면, 궤도 기여는 특정 에너지 창에서 ~100%에 달하는 현저히 높은 민감도를 보였다. 이는 궤도 자유도가 강유전성 편극에 매우 민감하다는 것을 나타낸다.

의의 및 주장
본 논문은 궤도 자유도를 포함하는 것이 vdW 이종 구조에서의 각운동량-전하 변환 및 THz 방출을 심도 있게 이해하는 데 필수적이라고 주장한다. 저자들은 자신들의 제일원리 시뮬레이션이 강유전성 편극에 의한 방출 조절에 대한 실험적 THz-TDS 결과를 정량적으로 설명한다고 단언한다.

주요 주장은 다음과 같다:

1. 상보적 메커니즘: 스핀과 궤도 응답은 음의 스핀-궤도 상관관계로 인해 특정 에너지 범위에서 서로 반대 부호를 가지지만, 이들이 결합하여 실험적 THz 신호와 일치하는 순(net) 양의 축적을 생성한다.
2. 궤도 민감도: 궤도 기여는 스핀 기여보다 편극 전기장에 더 민감하며, 이는 THz 방출 효율을 향상시키거나 튜닝할 수 있는 잠재적인 경로를 제공한다.
3. 방법론적 프레임워크: 본 연구는 강유전성 전기장이 존재하는 환경에서 비평형 스핀 및 궤도 밀도를 계산하는 방법론을 확립하였으며, 이는 vdW 이종 구조의 물성 설계에 적용될 수 있다.
4. 향후 방향: 저자들은 전기적 구동 응답은 정량화되었으나, 페로마그네틱 측으로부터 주입되는 정확한 궤도 편극량은 계면 상세 구조에 따라 달라지는 미지수임을 언급하였다. 이들은 니켈(Ni)과 같은 원소를 사용하는 등의 재료 공학을 통해 궤도 기여를 극대화하여, 궤도 응답이 지배적인 스핀트로닉스 샘플을 구현할 수 있을 것이라고 제안한다.

결론적으로, 본 연구는 $\text{PtSe}_2$와 같이 스핀 퇴화(degeneracy)를 해제하여 높은 효율을 얻을 수 있는 중심 대칭 물질에서, 스핀과 궤도 자유도 모두의 역학을 이해하는 것이 차세대 THz 기술 및 스핀-궤도 토크 소자 개발에 필수적임을 밝히고 있다.