바륨지르코네이트에서 프로톤 확산을 위한 반응성 힘장 모델

초록

**
본 연구는 BaZrO₃와 Y‑도핑된 BaZrO₃ 내 프로톤 확산을 기술하기 위해, 양자역학 계산으로부터 파라미터를 추출한 비반응성 힘장을 구축하고, 이를 경험적 원가 결합(EVB) 모델과 결합한 새로운 반응성 힘장을 제안한다. EVB‑MD 시뮬레이션을 통해 0.42 eV의 활성화 에너지와 실험값에 근접한 확산 계수를 얻었으며, 저농도에서 Y 도핑이 확산에 미치는 영향은 미미함을 확인하였다.

**

상세 분석

**
이 논문은 고체산화물 연료전지(SOFC) 핵심 물질인 BaZrO₃의 프로톤 전도 메커니즘을 원자 수준에서 정량적으로 이해하고자 하는 시도이다. 기존의 전통적인 고전력장(예: ReaxFF) 모델은 파라미터 수가 방대하고, 특히 BaZrO₃가 DFT 계산에서 비정방형(삼방형) 구조로 최적화되는 문제 때문에 실제 실험에서 관찰되는 입방체 구조를 재현하지 못한다는 한계가 있었다. 저자들은 이러한 문제를 회피하기 위해, 먼저 고정 전하와 Buckingham 포텐셜을 이용한 단순 강체 이온(force‑field) 모델을 구축하였다. 여기서 Ba²⁺, Zr⁴⁺, Y³⁺는 공식 전하를 부여하고, O²⁻와 H⁺는 부분 전하를 할당해 전기적 상호작용을 기술한다. 단거리 상호작용은 A·exp(−r/ρ) 형태의 Buckingham 포텐셜으로, 분산항(C₆)은 모두 0으로 설정해 전자구름의 장거리 인력을 배제하였다. 이는 DFT‑GGA(AM05) 계산에서 얻은 격자 상수와 결합 길이, 그리고 실험적인 탄성 상수를 동시에 만족시키도록 least‑squares 최적화가 수행되었다.

수소와 산소 사이의 수산기(OH⁻)를 묘사하기 위해, 저자들은 O‑H 결합을 조화 진동자와 H‑O 간의 순수 반발 Lennard‑Jones 포텐셜으로 분리하였다. 특히 Y‑도핑이 존재할 경우, Y에 인접한 산소(O₃)와 수산기 사이의 반발 파라미터를 강화해, Y‑인접 산소가 프로톤을 “덫어두는” 현상을 재현하였다. 이러한 비반응성 파라미터는 이후 EVB(Empirical Valence Bond) 모델에 의해 동적으로 전이한다. EVB는 두 개의 전자 상태(프로톤이 각각 인접한 두 산소에 결합된 상태)를 2×2 해밀토니안으로 연결하고, 오프다이아고날 원소 H₁₂=H₂₁를 거리 의존적인 Gaussian 형태로 정의해 프로톤 전이 장벽을 조절한다. 이는 물속 Grotthuss 메커니즘에서 다중 상태 EVB가 필요하다는 기존 연구와 달리, 고체 격자 내에서는 단일 전이 경로만을 고려해도 충분함을 전제로 한다.

MD 시뮬레이션은 확장된 등온등압(isostress‑isothermal) 앙상블을 사용해, stochastic thermostat과 barostat을 동시에 적용하였다. 이는 기존 NPT 구현이 고체의 비등방성 응력 텐서를 제대로 반영하지 못하는 문제를 해결한다. 시뮬레이션은 300 K–800 K 범위에서 수행됐으며, 프로톤의 평균 제곱 변위(MSD)를 통해 확산 계수를 추출하였다. 결과적으로, 저농도(프로톤 1 % 이하)에서 Y‑도핑이 존재해도 확산 계수에 유의미한 차이가 없었으며, 이는 프로톤이 Y‑인접 산소에 잠시 머무르지만, 탈출 빈도가 높아 전체 확산에 큰 영향을 주지 않기 때문이다. 활성화 에너지는 0.42 eV로, 실험값(≈0.45 eV)과 매우 근접했지만, 전형적인 Arrhenius 프리팩터는 약 30 % 정도 낮게 예측되었다. 이는 EVB 파라미터가 전이 상태의 진동 자유도를 완전히 포착하지 못했을 가능성을 시사한다.

마지막으로, 저자들은 ReaxFF와 비교해 모델의 장점과 한계를 논의한다. ReaxFF는 복잡한 파라미터화와 비정방형 구조 예측 문제를 안고 있지만, 다중 결합 파괴·재결합을 자연스럽게 다룰 수 있다. 반면 현재 제안된 EVB‑FF는 파라미터 수가 적고, 계산 비용이 낮으며, 실험적인 입방체 구조를 그대로 유지한다는 점에서 고체 전해질 시뮬레이션에 적합하다. 다만, 고농도 도핑이나 grain boundary와 같은 구조적 결함을 포함한 경우, 다중 상태 EVB 혹은 추가적인 반응 좌표가 필요할 수 있다.

**