상압 60K 초전도 리튬 도핑 LiMgZr2H12 설계

초록

본 연구는 MgZrH6 구조에 리튬을 도핑하여 LiMgZr₂H₁₂를 설계하고, 첫 원리 계산을 통해 상압에서 열역학·기계·동적 안정성을 확인하였다. 전자‑음자 결합 분석과 전자‑포논 결합(EPC) 계산 결과, λ = 2.22, ω_log ≈ 396 K, N(E_F) ≈ 2.34 states/eV·atom을 얻어 Allen‑Dynes 식으로 Tc ≈ 60.8 K를 예측하였다. 또한 초전도 성능 지표 S = 1.56으로 기존 MgZrH₆(36 GPa)보다 우수함을 입증하였다.

상세 분석

LiMgZr₂H₁₂는 Pmmm 대칭을 갖는 사중 원소 수소화물로, MgZrH₆의 1 × 1 × 2 초셀에 Mg 원자를 하나 Li로 치환해 얻었다. DFT‑PBE 계산에서 600 eV 평면파 컷오프와 2π × 0.02 Å⁻¹ k‑점 격자를 사용해 구조 최적화를 수행했으며, 형성 에너지가 음수이고 Born 안정성 기준을 만족함으로써 열역학·기계적 안정성을 확인하였다. PHONOPY를 이용한 포논 스펙트럼에서는 모든 파동수가 양의 값을 보여 동적 안정성도 확보되었다.

ELF와 Bader 전하 분석 결과, Li, Mg, Zr가 각각 약 0.86 e, 1.63 e, 1.66 e를 전자와 공유하고, H는 0.41–0.55 e를 받아 금속‑수소 이온 결합이 주된 결합 형태임을 알 수 있다. COHP 분석에서는 Zr‑H 결합이 가장 강하고, Li‑H·Mg‑H 결합도 음의 피크를 나타내어 전자‑음자 결합에 기여한다. H‑H COHP는 거의 0에 가까워 H 원자 간 직접 결합이 없으며, 이는 수소가 단원자 형태로 존재함을 의미한다.

전자 구조에서는 Fermi 레벨 근처에 Y, S, X 점에서 뱅 호프 특이점이 나타나며, 전체 DOS는 H와 Zr 기여가 크게 증가한다. 특히 H‑유래 상태가 N(E_F)의 약 70 %를 차지해 전자‑포논 결합을 강화한다. EPC 계산에서는 저주파(0–9 THz) 영역이 Li, Mg, Zr 진동에 의해 λ의 48 %를, 중·고주파(9–40 THz) 영역이 H 진동에 의해 52 %를 담당한다. λ = 2.22는 MgB₂(λ ≈ 0.61)와 MgZrH₆(λ ≈ 1.13)보다 현저히 크며, ω_log = 396 K는 가벼운 H 원자에 기인한다. Allen‑Dynes 수정 McMillan 식에 μ* = 0.10을 적용하면 Tc ≈ 60.8 K가 도출된다.

Fermi surface는 네 개의 밴드가 교차하는데, n = 3 밴드의 rhombic 형태와 n = 4 밴드의 평면이 거의 평행해 nesting이 형성된다. 이는 Kohn anomaly을 유발해 특정 포논 모드가 연화될 가능성을 높이며, EPC 강화와 Tc 상승에 기여한다.

마지막으로, 기존 고압 합성 사례(6 MgH₂ + ZrH₂ + n LiH)와 비교해 LiMgZr₂H₁₂는 약 170 meV/atom의 메타안정성을 보이지만, 촉매와 급냉, 불활성 분위기 등으로 실험적 합성이 가능할 것으로 제시한다.