라니오니트 초전도 얇은막의 초교환과 전하 이동

초록

본 연구는 압력 없이 40 K 이상 초전도성을 보이는 라니오니트(La₃Ni₂O₇) 얇은막의 초교환 상호작용과 전하 전달 특성을 11밴드 d‑p Hubbard 모델로 분석한다. 대규모 결정론적 양자 몬테카를로와 셀룰러 동적 평균장 이론을 이용해 얇은막에서의 수직 초교환 J⊥가 약 27 % 감소하고, 전하 전이 갭이 얇아짐을 확인하였다. 또한 전자와 정공 도핑 시 궤도별 전하 분포가 비대칭임을 밝혀, 얇은막과 벌크 사이의 물리적 차이를 정량적으로 규명하였다.

상세 분석

본 논문은 최근 보고된 La₃Ni₂O₇(LNO) 얇은막의 상온 초전도 현상을 이론적으로 검증하기 위해, 4개의 Ni‑3d(e_g) 궤도와 7개의 O‑2p 궤도를 포함하는 11밴드 d‑p Hubbard 모델을 구축하였다. 모델 파라미터는 압축된 격자를 갖는 얇은막에 대한 DFT 계산 결과를 직접 가져와, 벌크와 구조적 차이를 반영하도록 설계되었다. 전자 상호작용은 Ni‑d 궤도에만 적용되는 Hubbard U(기본값 7 eV)와 Hund’s coupling J_H(=0.15U)으로 기술했으며, 이중점유 보정(E_DC)은 Held식 공식으로 처리하였다.

수치 해석은 두 가지 상보적 방법으로 수행되었다. 첫째, 6×6 단위셀(총 396 오비탈)까지 확장 가능한 결정론적 양자 몬테카를로(DQMC) 시뮬레이션을 고온(T≈0.25 eV)에서 수행했으며, 평균 부호 ⟨s⟩≈0.79를 유지하였다. 둘째, 2×2 클러스터(44 오비탈)를 이용한 셀룰러 동적 평균장 이론(CDMFT)으로 저온(T≈0.08 eV)까지 접근했다. 두 방법 모두 스핀-스핀 상관함수 ⟨S_i·S_j⟩를 직접 측정해 초교환 상수 J⊥(수직)와 J∥(수평)를 추정하였다.

결과는 다음과 같다. 원자극한에서의 2차 섭동 이론에 따르면, U=7 eV일 때 J⊥≈0.135 eV, J∥≈0.084 eV이며, J∥/J⊥≈0.62이다. DQMC와 CDMFT 모두에서 수직 d_{3z²‑r²}–d_{3z²‑r²} 초교환이 가장 강하지만, 얇은막에서는 벌크 대비 약 27 % 감소하였다(⟨S·S⟩≈‑0.057→‑0.041 수준). 반면 수평 d_{x²‑y²}–d_{x²‑y²} 초교환은 거의 변하지 않아, 얇은막에서는 두 초교환이 비슷한 크기로 남아 있다. 온도 의존성 분석에서 얇은막의 J⊥는 온도가 낮아질수록 급격히 약해지는 반면, J∥는 비교적 완만하게 감소한다. 이는 얇은막에서 수직 초교환이 억제되어 초전도 전이온도 T_c가 벌크보다 낮아지는 원인으로 해석될 수 있다.

전하 전달 특성에 대해서는 µ‑n_h 곡선을 조사해 반전점 µ_h≈‑1.55에서 압축된 얇은막이 벌크보다 완만한 경사를 보이며, 전하 전이 갭이 약 0.2 eV 정도 작아짐을 확인했다. 도핑 시 전하 분포를 궤도별로 분석하면, 정공 도핑에서는 O‑p_x/p_y와 Ni‑d_{x²‑y²}에 약 1:1 비율로 정공이 들어가고, d_{3z²‑r²}에는 상대적으로 적게 들어간다. 전자 도핑에서는 Ni‑d 궤도, 특히 d_{x²‑y²}와 d_{3z²‑r²}에 전자가 주로 채워지며, d_{3z²‑r²}가 약 3배 더 많이 차지한다. 이러한 입자‑정공 비대칭은 Zaanen‑Sawatzky‑Allen 전하 전이 모델에서 중요한 역할을 하며, 얇은막에서 d_{3z²‑r²} 궤도의 점유율이 초전도 메커니즘(Hund’s coupling 주도 vs. e_g 혼성) 선택에 결정적임을 시사한다.

마지막으로, 얇은막의 초교환 감소와 전하 전이 갭 축소는 저에너지 t‑J 모델을 구축할 때 J⊥와 J∥의 비율을 재조정해야 함을 의미한다. 이는 얇은막 특유의 양자 임계 현상과 초전도 페어링 대칭성을 이해하는 데 필수적인 입력값이 된다.