바륨‑라탄‑황 및 갈륨‑라탄‑황 시스템의 열역학적 평가와 상도 예측

초록

본 논문은 CALPHAD 방법과 AIMD 계산을 결합하여 Ba–La–S와 Ga–La–S 이진·준이진 시스템의 열역학적 모델을 구축하고, 실험 데이터와 비교해 신뢰성 있는 상도와 물성 데이터를 제공한다.

상세 분석

본 연구는 희토류 황화물계의 열역학적 데이터가 극히 부족한 상황에서, CALPHAD 프레임워크를 활용해 Ba–La–S와 Ga–La–S 시스템의 Gibbs 자유에너지 모델을 체계적으로 구축한 점이 가장 큰 특징이다. 먼저 기존 문헌에서 보고된 Ba–S, Ba–La, La–S, Ga–S 등 네 개의 이진계에 대한 상변화와 열역학적 특성을 정밀히 검토하고, 특히 Ba‑S와 Ga‑S 시스템에서 고온에서의 황 휘발성 및 용기 반응성으로 인해 실험 데이터가 제한적인 점을 강조한다. 이러한 데이터 부족을 보완하기 위해 저자들은 ab initio 분자동역학(AIMD) 시뮬레이션을 수행해 액체 상의 혼합 엔탈피를 직접 계산하였다. AIMD 결과는 특히 Ba‑S와 Ga‑S 액체의 짧은 거리 순서(short‑range ordering) 현상을 반영하는 Redlich‑Kister 형태의 교환 매개변수로 전환되었으며, 이는 기존 실험적 액체상 데이터와 일관된 경향을 보였다.

상도 모델링에서는 고체 용해도, 전이 상, 그리고 다형성(polymorph) 전이를 모두 포함하는 복합적인 상 모델을 도입하였다. 예를 들어 La₂S₃의 α, β, γ 3가지 다형성을 각각 독립적인 Gibbs 에너지 함수로 정의하고, β‑La₂S₃의 산소 함량에 따른 비정질 고체용액 모델을 적용하였다. 또한 BaS₂와 BaS₃ 같은 고온 상에 대해서는 전이 온도와 엔탈피 변화를 실험값에 맞추어 최적화하였다. Ga‑S 시스템에서는 논란이 많았던 Ga₂S와 Ga₄S₅ 같은 중간 화합물을 배제하고, GaS와 Ga₂S₃만을 주요 상으로 채택했으며, 두 개의 액체‑액체 상분리(Liquid‑Liquid miscibility gap)를 모델링하여 고온에서의 비균일 용해 현상을 재현하였다.

모델 파라미터 최적화는 Thermo-Calc 기반의 CALPHAD 소프트웨어를 이용해 수행되었으며, 각 상의 표준 상태 엔탈피, 엔트로피, 열용량을 문헌값과 AIMD 결과에 동시에 피팅하였다. 결과적으로 도출된 이진·준이진 상도는 기존 실험 데이터와 매우 높은 일치도를 보였으며, 특히 BaS‑La₂S₃, Ga₂S₃‑La₂S₃ 준이진 구간에서 예측된 공융점(eutectic)과 전이 반응이 실험적으로 보고된 온도와 거의 동일했다. 이러한 성공은 짧은 거리 순서 모델과 AIMD 기반 액체 혼합 엔탈피가 CALPHAD 모델에 효과적으로 통합될 수 있음을 증명한다.

본 연구는 희귀 금속 황화물 시스템에 대한 열역학 데이터베이스를 최초로 제공함으로써, 고온 적외선 광학재료, 반도체 소자, 보호 코팅 등 다양한 응용 분야에서 물성 예측과 합성 설계에 직접 활용될 수 있다. 또한 제한된 실험 데이터에 의존하는 전통적인 접근법을 보완하는 계산‑실험 융합 전략의 모범 사례로 평가될 수 있다.