Reactive Force Field for P/Sn/I System: Atomistic Insight into the Early Stage of Black Phosphorus and Phosphorene Synthesis Process

기술 요약: P/Sn/I 시스템을 위한 반응성 포스 필드 및 흑린 합성의 원자론적 통찰

문제 제기
흑린(Black phosphorus, BP)과 그 2차원 형태인 포스포린(phosphorene)은 높은 전하 이동도와 조절 가능한 밴드 갭으로 인해 차세대 전자 소자를 위한 유망한 재료이다. 그러나 화학 기상 수송법(Chemical Vapor Transport, CVT)을 통한 고품질의 대면적 BP 박막 합성은 여전히 도전적인 과제이다. 일반적으로 주석(Sn)과 요오드(I₂) 첨가제를 포함한 적린(red phosphorus)을 사용하는 이 공정은 원자론적 수준에서 완전히 이해되지 않았다. 첨가제의 구체적인 역할, 안정적인 기상 종(species)의 정체, 인(phosphorus) 수송 메커즘, 그리고 전핵 형성 중간체의 성질에 관한 핵심적인 불확실성이 지속되고 있다. 고온 및 반응성 환경으로 인해 초기 단계의 역학을 실험적으로 관찰하는 것은 어려우며, 표준 밀도 범함수 이론(DFT)은 작은 시스템과 정적인 0K 구성에 국한되어 있어, 성장 메커니즘을 모델링하는 데 필요한 동적이고 유한한 온도의 반응 네트워크를 포착하는 데 한계가 있다.

방법론
이러한 한계를 해결하기 위해, 저자들은 P/Sn/I-2512로 명명된 P/Sn/I 시스템 전용의 새로운 반응성 포스 필드(ReaxFF)를 개발하였다. 이 포스 필드는 DFT(PBE 범함수, TZ2P 기저 집합 및 Grimme의 D3 분산 보정 사용)를 통해 생성된 포괄적인 양자 역학(QM) 데이터셋을 바탕으로 학습되었다. 학습 데이터 세트에는 다음이 포함된다:

• Sn–Sn, Sn–I, Sn–P, I–P, I–I 결합의 결합 해리 에너지.
• 다양한 분자 구성에 대한 원자가 각도 및 비틀림 프로파일.
• 주요 종(SnI₂, SnI₄, PI₃)의 생성열.
• 응축상 주석(α-Sn 및 β-Sn)의 상태 방정식 및 주석 클러스터(Sn$_x$, x=2-6)의 형성 에너지.

포스 필드 파라미터는 Amsterdam Modeling Suite 내에서 공분산 행렬 적응 진화 전략(CMA-ES)을 사용하여 최적화되었다. 검증은 결합 해리, 상태 방정식 및 반응 에너지에 대한 ReaxFF 결과와 DFT 데이터를 비교함으로써 수행되었으며, 강력한 일치성을 보였다.

이어서, 저자들은 923 K에서 초기 단계의 전핵 형성 과정을 조사하기 위해 LAMMPS를 사용하여 반응성 분자 역학(ReaxFF-MD) 시뮬레이션을 수행하였다. 연구에는 Nosé-Hoover 열조절기를 사용하는 NVT 앙상블이 사용되었다. 시뮬레이션은 특정 변수를 격리하기 위해 다음과 같이 구조화되었다:

1. 시스템 검증: 다양한 밀도(유효 압력 0.01, 0.1, 1 atm)에서 Sn/I₂ 및 P/I₂ 시스템의 MD 시뮬레이션을 수행하여 SnI₂, SnI₄, PI₂, PI₃와 같은 주요 종의 형성을 확인하였다.
2. 요오드의 역할: 요오드가 인의 재결합에 미치는 영향을 격리하기 위해 세 가지 시스템(P 단독; P + Sn; P + Sn + I₂)을 비교하였다.
3. 밀도 효과: BP 형성이 실험적으로 관착되는 조건을 모사하기 위해, 희박한 증기에서 응축상으로의 전이를 모델링하고자 P/Sn/I₂ 시스템을 세 가지 밀도 영역(유효 압력 5, 10, 100 atm)에서 시뮬레이션하였다.

주요 결과

• 포스 필드 검증: P/Sn/I-2512 포스 필드는 결합 해리, 원자가 각도, 비틀림 장벽 및 주석 결정과 클러스터의 상태 방정식에 대한 DFT 참조 데이터를 정확하게 재현하였다. 또한 고온에서 SnI₂, SnI₄ 및 PI₃의 형성을 성공적으로 예측하였다.
• Sn/I₂ 및 P/I₂에서의 밀도 의존성: Sn/I₂ 시스템에서 밀도가 증가함에 따라 조성은 원소 종에서 복잡한 배위 환경으로 이동하였으며, SnI₂와 함께 SnI₄ 및 Sn₂I₄(SnI₂의 이량체)의 형성을 촉진하였다. P/I₂ 시스템에서는 높은 밀도가 PI₃ 및 고분자 PI$_x$ 종의 형성을 용이하게 하였는데, 이는 순수한 PI₃를 분리하기 어려운 실험적 어려움과 일치한다.
• 재결합에서의 요오드의 역할: 시뮬레이션 결과, 요iod가 없는 경우 인 원자들은 주로 원자 상태로 남아 있거나 최소한의 클러스터만을 형성하였다. I₂의 첨가는 P–P 재결합 속도를 크게 증가시켜, 이진 P$_x$I$_y$ 모티프와 일시적인 삼원계 Sn$_x$P$_y$I$_z$ 화합물의 형성을 유도하였다. 이는 인을 용해시키고 수송을 용이하게 하는 요오드의 역할을 확인시켜 준다.
• 고밀도에서의 전핵 경로: 낮은 밀도(5 atm)에서 시스템은 주로 원자 상태의 인과 작은 클러스터(P₂, P₃)로 구성되었다. 중간 밀도(10 atm)에서는 더 큰 클러스터(P₄, P₅)가 형성되기 시작하였다. 높은 밀도(100 atm)에서 시스템은 고차 P$_x$ 클러스터(x = 8–12) 및 거대한 삼원계 Sn$_x$P$_y$I$_z$ 화합물로 빠르게 응축되었다.
• 중간체의 식별: 가장 높은 밀도에서 시뮬레이션은 Sn$_x$P$_y$I$_z$ 클러스터의 가장자리에서 응축된 Hittorf형 인 구조를 생성하였다. 또한, 요오드에 의해 안정화된 BP 유사 육각형 모티프가 관찰되었으며, 이는 잠재적인 핵 생성 씨앗(nucleation seeds)으로 식별되었다.

의의 및 주장
본 논문은 BP/포스포린 증기 합성을 위한 성장 조건 및 첨가제 효과를 연구하기 위한 실용적인 원자론적 모델을 제공한다고 주장한다. 주요 기여는 다음과 같다:

1. 검증된 포스 필드 개발: P/Sn/I-2512 포스 필드는 DFT의 크기 및 시간 척도 제한을 극복하여, 인, 주석, 요오드를 포함하는 복잡한 반응계 시뮬레이션을 위한 신뢰할 수 있는 도구를 제공한다.
2. 메커니즘적 통찰: 본 연구는 요오드와 시스템 밀도가 인의 재결합을 조절하는 핵심 요인임을 밝혀냈다. 요iod가 단순히 수송제 역할만 하는 것이 아니라, 중간체인 P$_x$I$_y$ 및 삼원계 Sn$_x$P$_y$I$_z$ 종의 형성을 통해 P–P 재결합을 능동적으로 촉진한다는 것을 입증하였다.
3. 전핵 모델: 저자들은 고밀도 조건과 요듐의 결 combination이 인의 응축을 촉진하여 Hittorf형 중간체와 BP 유사 육각형을 형성하는 구체적인 전핵 경로를 제안한다. 이러한 종들은 냉각 단계 동안 이후의 흑린 성장을 위한 전구체 역할을 하는 것으로 제안된다.

본 연구는 새로운 실험 프로토콜을 제안하는 것이 아니라, 흑린의 CVT 합성에서 고압(밀도)의 필요성과 요iod의 결정적인 역할을 포함하여 기존의 실험적 관찰을 해석하기 위한 이론적 틀을 제공한다.