NbSe₂ 이중층에서 Mn 삽입 및 확산 메커니즘 DFTB 시뮬레이션

초록

본 연구는 밀도 함수 기반 타이트 바인딩(DFTB) 방법을 이용해 NbSe₂ 2‑층 구조에 망간(Mn) 원자의 삽입(인터칼레이션)과 확산 경로를 원자 수준에서 조사하였다. 결합 에너지 계산 결과 Mn은 표면 흡착보다 층 사이(intercalated)와 층 내부(embedded) 위치에서 더 안정적이며, NEB 계산을 통해 층 사이로 이동하는 에너지 장벽이 0.68 eV임을 확인했다. 525 K에서 20 ps MD 시뮬레이션은 Mn 농도가 증가할수록 표면에서 내부로 점진적으로 이동해 최종적으로 완전한 인터칼레이션이 이루어지지만, 과도한 농도에서는 클러스터링이 발생한다는 점을 보여준다.

상세 분석

본 논문은 전이금속 디칼코게나이드(TMD)인 NbSe₂ 이중층에 Mn 원자를 삽입했을 때의 구조적·열역학적 안정성 및 동역학적 거동을 체계적으로 규명한다. 먼저 DFTB+와 최신 PTBP 파라미터셋을 활용해 NbSe₂ 3‑차원 및 2‑차원 구조를 최적화했으며, 실험값(a ≈ 3.608 Å, c ≈ 13.44 Å)과 거의 일치하는 격자 상수를 얻어 계산 방법의 신뢰성을 확보하였다. 이후 4 × 4 초셀에 Mn을 다섯 가지 가능한 위치(표면 hollow, Nb‑top, Se‑top, 층 사이 인터칼레이션, 층 내부 embedding)에서 배치하고, 결합 에너지(E_coh)를 식(1)로 정의해 비교하였다. 결과적으로 인터칼레이션(−5.573 eV)과 embedding(−5.639 eV)이 표면 흡착(Se‑top −3.129 eV, hollow −4.483 eV, Nb‑top −5.060 eV)보다 현저히 낮은 에너지를 보여, Mn이 층 사이 혹은 내부에 자리 잡는 것이 열역학적으로 우선함을 확인했다.

동역학적 장벽은 NEB 계산으로 평가했으며, hollow → 인터칼레이션 전이 경로에서 최대 0.68 eV의 장벽을 관찰했다. 이는 Mn/graphene·Ge(001)·GaAs(001) 등 기존 보고된 시스템(0.1–0.8 eV)과 비슷한 수준으로, 결함이나 에지 사이트가 없어도 자연스럽게 van der Waals 간극으로 확산될 수 있음을 시사한다.

MD 시뮬레이션에서는 525 K(≈250 °C)에서 20 ps 동안 다양한 Mn 농도(1, 4, 8, 10, 12 원자)를 조사하였다. 저농도(1 원자)에서는 Mn이 표면에 머무르면서 서서히 내부로 이동하려는 경향을 보였고, 4–8 원자에서는 부분적인 내부 침투가 관찰되었지만 완전한 인터칼레이션은 이루어지지 않았다. 10 원자에서는 두 개의 Mn이 완전히 층 내부로 이동해 실험적 인터칼레이션 결과와 일치했으며, 12 원자에서는 과포화에 따른 클러스터링 현상이 나타나 추가 확산을 억제했다. 이러한 농도 의존적 거동은 실제 성장 과정에서 결함이나 초과 공급량이 인터칼레이션 효율에 미치는 영향을 설명한다.

전반적으로 연구는 (1) Mn이 NbSe₂ 이중층에 삽입될 때 가장 안정한 위치는 층 사이와 내부이며, (2) 전이 장벽이 0.68 eV로 비교적 낮아 열적 활성화가 용이하고, (3) 충분한 Mn 농도가 확보될 경우 자연스럽게 인터칼레이션이 진행되지만 과도한 농도는 클러스터링을 초래한다는 세 가지 핵심 인사이트를 제공한다. 이러한 결과는 Mn 도핑을 통한 NbSe₂의 전자·자성 특성 조절, 스핀트로닉스 및 vdW 이종구조 설계에 실질적인 가이드라인을 제공한다.