플래티넘을 위한 차세대 각도‑의존 및 변형 테르소프 포텐셜 개발

초록

본 연구는 플래티넘(Pt) 원자 간 상호작용을 묘사하기 위해 각도‑의존 포텐셜(ADP)과 변형 테르소프 포텐셜(MT) 두 가지 새로운 모델을 제시한다. 두 포텐셜은 실험 데이터를 전혀 사용하지 않고, 광범위한 DFT 기반 레퍼런스 데이터베이스를 이용해 파라미터를 최적화하였다. 기존 EAM 포텐셜(Zhou, O’Brien) 대비 결정 구조, 탄성 상수, 결함 에너지, 표면·쌍곡면 에너지, 용융점, 포논 스펙트럼 등 다양한 물리량에서 향상된 정확도를 보이며, 특히 MT 모델은 금속‑공유 결합 혼합 시스템에의 확장 가능성을 논의한다.

상세 분석

이 논문은 플래티넘의 원자 간 상호작용을 정밀하게 기술하기 위해 두 가지 새로운 포텐셜 형식을 도입한 점에서 의미가 크다. 첫 번째인 ADP는 전통적인 EAM에 각도‑의존 항을 추가함으로써 비구면 전자밀도와 결합각의 영향을 반영한다. 이는 특히 전이 금속에서 관찰되는 비대칭 결함 구조와 전단 변형에 대한 정확한 에너지 예측에 유리하다. 두 번째인 MT는 원래 강한 공유 결합을 갖는 실리콘·탄소 계열을 위해 고안된 테르소프 포텐셜을 금속에 적용하기 위해 절단 반경을 다중 배위층까지 확대한 변형 형태이다. 이렇게 하면 MT가 MEAM과 유사한 다체 상호작용을 포괄하면서도, 파라미터화 과정에서 각도‑의존성을 자연스럽게 포함한다는 장점이 있다.

훈련 데이터베이스는 FCC, BCC, HCP, SC, A15, 다이아몬드 구조 등 6가지 결정구조에 대한 EOS, 비등방성 변형(축압·인장, 전단), 결함 경로(NEB 기반 빈자리 이동), 표면·쌍곡면 에너지, 그리고 고체·액체 상 coexistence를 위한 고온 구조까지 포괄한다. DFT 계산은 VASP와 PBE‑PAW를 사용했으며, 에너지·힘 수렴 기준을 엄격히 설정해 데이터의 신뢰성을 확보하였다.

포텐셜 파라미터 최적화는 목표 함수에 DFT 에너지·힘을 동시에 포함시켜, 구조적 안정성뿐 아니라 동역학적 특성(포논, 열팽창)까지 재현하도록 설계되었다. 결과적으로 ADP와 MT는 기존 EAM1·EAM2에 비해 (i) 격자 상수와 결합 에너지에서 <1 % 오차, (ii) C₁₁, C₁₂, C₄₄ 등 탄성 상수에서 실험값과 <5 % 차이, (iii) 빈자리·전위 결함 형성 에너지와 이동 장벽에서 DFT와 거의 일치, (iv) (111)·(100) 표면 에너지와 γ‑surface 에서도 실험·DFT와 일관된 값을 제공한다. 특히 MT는 용융점 예측이 1760 K(실험 1775 K)로, 기존 EAM이 1500 K 이하로 과소평가하던 문제를 크게 개선하였다.

전이 금속 특유의 전자밀도 변동을 각도‑의존 항으로 보정한 ADP는 고체 내 결함·전단 거동을 정밀히 모사하는 데 강점이 있다. 반면 MT는 절단 반경을 확대함으로써 다배위층 상호작용을 포괄하고, 금속‑공유 결합 혼합 시스템(예: Pt‑Si, Pt‑C 합금)에도 적용 가능성을 시사한다. 다만 두 포텐셜 모두 훈련 데이터에 포함되지 않은 고압·고온 비정질 구조에 대한 검증이 부족하며, 파라미터가 복잡해 계산 비용이 기존 단순 EAM보다 다소 증가한다는 점은 향후 연구에서 보완이 필요하다.