온도에 따라 변하는 결함 전하 상태와 전하 전이 레벨

초록

본 연구는 머신러닝 기반 인터액션 포텐셜과 열역학적 적분을 이용해 MgO, LiF, CsSnBr₃의 빈자리 결함에 대한 온도 의존적 전하 전이 레벨(CTL)을 계산한다. 진동 엔트로피와 밴드 엣지 이동이 CTL을 크게 이동시키며, 특히 CsSnBr₃에서는 60 K 이상에서 중성 전하 상태가 열역학적으로 안정화되어 새로운 Fermi‑레벨 구간이 생성된다.

상세 분석

이 논문은 전하 전이 레벨(CTL)이 0 K에서만 정의되는 전통적 결함 모델의 한계를 지적하고, 온도에 따른 자유에너지 변화를 정량화하는 새로운 워크플로우를 제시한다. 핵심은 DFT‑학습 신경진화 포텐셜(NEP)을 이용해 대규모 분자동역학(MD) 시뮬레이션을 수행하고, 두 가지 열역학적 적분 경로—Frenkel‑Ladd와 직접 전하‑전하 경로—를 통해 Gibbs 자유에너지 차이를 얻는 것이다. NEP 모델은 결함 주변 원자를 전하 상태에 따라 라벨링함으로써 동일한 원자 배열에서도 전하 상태를 구분하도록 설계되었다. MgO와 LiF는 PBEsol, CsSnBr₃은 r2SCAN 기반 데이터로 학습했으며, 각각 +2/0, –1/0/+1, –1/0/+1 전하 상태를 포함한다.

열역학적 적분 결과, MgO의 ε⁺²/0 레벨은 0 K에서 5.15 eV에서 시작해 1000 K까지 약 –0.785 eV 이동한다. 이 이동은 VBM 상승(≈0.414 eV)과 진동 엔트로피 기여(≈–0.371 eV/K)로 나뉘며, 진동 효과가 전체 변위의 절반을 차지한다. LiF에서도 ε⁺¹/0와 ε⁰/–1 레벨이 각각 –0.994 eV, –1.035 eV 정도 이동해, 높은 온도에서 전하 상태 전이가 크게 변한다. 가장 흥미로운 결과는 CsSnBr₃에서이다. ε⁺¹/0, ε⁰/–1, ε⁺¹/–1 전이 레벨이 각각 약 0.111 eV, 0.324 eV, 0.218 eV 이동했으며, 60 K 이상에서 중성 전하 상태가 열역학적으로 가장 안정화된다. 이는 0 K에서는 존재하지 않던 새로운 Fermi‑레벨 윈도우를 형성한다는 의미다.

계산 효율성 측면에서 직접 전하‑전하 적분 경로는 Frenkel‑Ladd보다 적은 샘플링으로 수렴하며, 두 경로가 제공하는 CTL 값은 거의 일치한다. 또한, 동일한 MD 궤적에서 전하 상태 q와 q′의 에너지 차이를 추출해 온도 의존적 광학 전이 레벨(E_opt)도 계산했으며, 이들 역시 VBM 이동과 진동 엔트로피에 의해 유사하게 변한다. 제로‑K 광학 전이는 GW₀와 Bethe‑Salpeter 계산과 좋은 일치를 보이며, 온도 상승에 따라 전이 에너지가 넓은 분포를 갖는 것이 확인되었다.

이 연구는 (i) 머신러닝 포텐셜을 이용한 대규모 MD를 통해 정확한 자유에너지 계산이 가능함을, (ii) 진동 엔트로피와 밴드 엣지 변동이 CTL에 미치는 영향을 정량화함을, (iii) 온도에 따라 전하 상태 자체가 변할 수 있음을 입증한다. 따라서 전자소자, 광전지, 고체 이온 전도체 등에서 결함 특성을 예측할 때, 0 K 정적 모델만으로는 충분하지 않으며, 온도 효과를 포함한 전반적인 자유에너지 해석이 필수적임을 강조한다.