무독성 Mg₃ZBr₃ 페로브스카이트의 안정성 및 광전특성

초록

본 연구는 납 기반 페로브스카이트의 독성·불안정을 대체할 무독성 Mg₃ZBr₃(Z=As, Sb, Bi) 화합물을 첫 원리 계산으로 조사한다. Pm‑3m 구조의 입방체를 확인하고, HSE06 하이브리드 함수로 Mg₃AsBr₃(2.06 eV), Mg₃SbBr₃(1.65 eV), Mg₃BiBr₃(1.52 eV)의 간접 밴드갭을 예측하였다. 광학 흡수는 밴드갭 근처에서 급격히 상승하고, 정적 유전율은 As→Sb→Bi 순으로 증가한다. 포논 분산에서 Mg₃AsBr₃·Mg₃SbBr₃는 허수 모드가 없어 동적 안정성을 보이며, 큰 그뤼네이센 파라미터는 연성 격자와 강한 비조화성을 시사한다. 탄성 상수와 볼크스 모듈러스는 44 GPa→35 GPa로 감소해 적당한 강성을 나타낸다. p‑i‑n 드리프트‑디퓨전 시뮬레이션은 계산된 밴드갭과 일치하는 전류‑전압 특성을 보여, 이들 물질이 안정적인 얇은 필름 포토다이오드 및 태양전지 후보임을 제시한다.

상세 분석

본 논문은 Mg₂⁺와 할로겐(Br⁻)이 결합한 A₃ZBr₃(A=Mg, Z=As, Sb, Bi) 구조를 전산적으로 탐색함으로써, 기존 Pb‑기반 페로브스카이트의 환경·안전 문제를 해결하고자 한다. 첫째, VASP와 PAW 방법을 이용해 PBE와 HSE06 두 단계의 구조 최적화를 수행했으며, 최종적으로 모든 화합물이 Pm‑3m(221) 공간군의 단순 입방체 구조를 갖는 것을 확인하였다. 격자 상수는 As→Sb→Bi 순으로 5.48 Å→5.66 Å→5.72 Å(HSE06)로 증가하고, 이는 pnictogen 이온 반경이 커짐에 따른 자연스러운 팽창이다. Goldschmidt tolerance factor(t)는 0.78→0.74→0.63으로 감소해, 특히 Bi 함유체는 기하학적 불안정성 한계에 가까워 옥타헤드라틀의 회전·왜곡 가능성을 시사한다.

전자 구조 분석에서는 HSE06 하이브리드 함수가 25 % 정확한 교환을 포함함으로써 PBE 대비 밴드갭을 0.5 ~ 0.8 eV 정도 상승시켰다. Mg₃AsBr₃는 2.0645 eV, Mg₃SbBr₃는 1.6533 eV, Mg₃BiBr₃는 1.5226 eV의 간접 밴드갭을 보이며, 밴드꼭짓점은 각각 Γ와 R점 사이에 위치한다. 전도대는 주로 Mg‑Br σ* 결합과 pnictogen‑p 궤도 혼성에 의해 형성되고, 가전자대는 pnictogen‑5p(또는 6p)와 Br‑4p의 혼성으로 구성된다. 특히 Sb와 Bi의 ns² 비활성 전자쌍이 스테레오활성 lone‑pair 효과를 발휘해 전자구름을 비대칭적으로 왜곡시키며, 이는 밴드갭 감소와 유전율 상승을 동시에 유도한다.

광학 특성에서는 PBE 기반 독립 입자 근사를 사용해 복소 유전함수를 계산했으며, 실효 유전율(ε₁)은 As→Sb→Bi 순으로 5.2→6.1→7.3 정도로 증가한다. 이는 pnictogen‑Br p‑p 혼성 강화와 전자구름의 비대칭성 확대가 광전자 결합을 강화하기 때문이다. 흡수 스펙트럼은 밴드갭 직후 급격히 상승하여 10⁴ cm⁻¹ 수준의 강한 광흡수를 보이며, 가시광선 영역(≈1.8–2.5 eV)에서도 충분히 흡수 가능한 점이 포토디텍터 및 태양전지에 유리하다.

포논 계산에서는 DFPT와 Phonopy를 이용해 3×3×3 초셀에서 전자기적 힘 상수를 구하고, Mg₃AsBr₃와 Mg₃SbBr₃는 전 영역에 허수 주파수가 없어 동적 안정성을 확인했다. 반면 Mg₃BiBr₃는 Γ점 근처에 작은 허수 모드가 나타나며, 이는 저온에서 구조 변형(예: 옥타헤드라틀 회전) 가능성을 암시한다. 전체 포논 밴드폭은 As(≈8 THz)→Sb(≈6 THz)→Bi(≈5 THz)로 감소하고, Grüneisen 파라미터 평균값(γ̅)은 2.7→3.2→2.8으로 높은 비조화성을 나타낸다. 이는 격자 팽창에 민감한 저주파 모드가 많아 열전도도가 낮고, 열확산이 억제되는 소프트-라티스 특성을 부여한다.

탄성 분석에서는 PBE 레벨에서 6×6×6 변형을 가해 스트레스‑텐서 행렬을 추출, Hill 평균을 통해 Bulk modulus(K)와 Shear modulus(G)를 계산했다. K는 Mg₃AsBr₃≈44 GPa, Mg₃SbBr₃≈39 GPa, Mg₃BiBr₃≈35 GPa이며, G 역시 16–12 GPa 수준으로 비교적 연성이다. Born 안정성 기준(C₁₁−C₁₂>0, C₁₁+2C₁₂>0, C₄₄>0)을 모두 만족해 기계적 안정성을 확인했다. 이러한 중간 정도의 강성은 박막 성장 시 균열 방지와 동시에 유연한 기판에 적용 가능함을 의미한다.

마지막으로, COMSOL 기반 1 µm 두께 p‑i‑n 디바이스 모델을 구축해 전자·정공 확산 계수를 HSE06 밴드갭에 맞춰 설정하고, 전압‑전류 특성을 시뮬레이션했다. 결과는 전류가 광학 밴드갭에 의해 제한되는 형태를 보이며, Mg₃SbBr₃가 가장 높은 광전류 밀도와 낮은 역전압을 나타냈다. 이는 밴드갭이 1.65 eV로 태양 스펙트럼과 최적 매칭되고, 전하 운반체 이동도가 상대적으로 높기 때문이다. 전반적으로, 계산된 전자·광학·열·기계적 특성이 모두 실용적인 포토디텍터 및 태양전지 소재 요건을 충족함을 보여준다.