인접 원소 기반 구조 탐색으로 새로운 칼코겐화합물 예측

초록

본 연구는 기존 ICSD 데이터베이스에 존재하는 구조 유형을 인접 원소(근접 원소) 시스템에서 차용하여, 황화물·셀레늄화물 계열의 새로운 이진 화합물을 효율적으로 탐색한다. 제한된 검색 영역에서 10개의 신규 화합물을 예측했으며, 근접 원소 기반 검색이 광범위한 전주기 탐색보다 계산 비용을 크게 절감한다는 점을 입증하였다.

상세 분석

이 논문은 재료 탐색의 핵심 병목인 “가능한 구조 후보의 폭발적 증가” 문제를 해결하기 위해, 구조‑조성 연관성을 이용한 휴리스틱 검색 전략을 제시한다. 저자들은 먼저 ICSD에 등재된 황화물·셀레늄화물·산화물의 이진 구조 유형을 정량적으로 분석하여, 서로 다른 화학계열 간에 33 % 정도의 구조 중복이 존재함을 확인하였다. 특히 같은 주기표 그룹에 속하는 원소쌍(예: Ti–Nb–Bi) 사이에서는 구조 공유 비율이 높아, 인접 원소 시스템에서 차용한 구조가 새로운 화합물의 안정성을 예측하는 데 유용함을 제시한다.

검색 절차는 세 단계로 구성된다. ① 최소 검색 영역에서는 목표 화합물(A‑S 혹은 A‑Se)과 동일한 A 원소를 갖는 반대 칼코겐(Se 혹은 S) 및 산화물(O) 구조를 차용한다. ② 근접 원소 영역에서는 A 원소를 같은 혹은 인접한 주기표 열에 있는 원소들로 교체한 구조를 추가한다. ③ 전주기 확장 영역에서는 모든 원소를 대상으로 구조를 차용한다. 각 단계마다 VASP‑SCAN 메타‑GGA를 이용해 전자 구조 최적화를 수행하고, 형성 엔탈피를 기반으로 Convex Hull를 구축해 열역학적 안정성을 평가하였다.

계산 결과, 첫 번째 제한된 영역만으로도 Ti₃S₄, Cu₃S₂와 같은 두 개의 새로운 안정 화합물과 10개의 메타안정 화합물을 발견하였다. 두 번째 근접 원소 영역에서는 Ta₃Se₄, CuSe, ZrSe, Zr₃Se₄, Zr₂Se₃ 등 다섯 개의 새로운 안정 화합물이 추가로 도출되었으며, 메타안정 후보도 37개가 확인되었다. 세 번째 전주기 확장에서는 Ti₅S₈ 한 종만이 추가로 발견되었으며, 이는 근접 원소 기반 탐색이 이미 대부분의 유망 후보를 포착했음을 의미한다.

효율성 측면에서, 첫 번째 단계는 48개의 구조 계산당 1개의 신규 화합물(≈ 2 % 성공률)이라면, 두 번째 단계는 205개의 계산당 1개의 신규 화합물(≈ 0.5 % 성공률)로 감소하였다. 전주기 전체를 탐색했을 때는 436개의 계산당 1개의 신규 화합물(≈ 0.23 % 성공률)로, 계산 비용 대비 수확률이 급격히 낮아짐을 보여준다. 이는 인접 원소 기반 구조 차용이 검색 효율을 크게 향상시킬 수 있음을 실증한다.

또한, 산화물 구조를 차용했을 때는 신규 화합물 발견이 거의 없었으며, 이는 황·셀레늄 계열과 산소 계열 사이의 구조적 연관성이 상대적으로 낮다는 기존 분석과 일치한다. 마지막으로, 칼코겐 자체를 다른 원소(예: Ti‑As 구조)로 교체하는 실험적 확장은 제한적인 성공을 보였으며, 원소 교체 거리가 멀어질수록 안정성 감소 경향이 뚜렷함을 확인하였다.

이러한 결과는 “구조‑조성 근접성”이라는 물리적 직관을 정량화하여, 대규모 고처리량 계산 없이도 새로운 물질을 효율적으로 탐색할 수 있는 실용적인 프레임워크를 제공한다는 점에서 의미가 크다. 향후 다른 화합물군(예: 탄소화물, 질화물)에도 동일한 전략을 적용한다면, 데이터베이스에 아직 등재되지 않은 다수의 잠재적 안정 화합물을 빠르게 발굴할 수 있을 것으로 기대된다.