SiAs 및 SiAs₂ 단일층의 변형 조절 밴드 갭과 탁월한 캐리어 이동도: 제1원리 통찰

초록

상당한 밴드 갭과 높은 캐리어 이동도를 겸비한 원자 수준 얇은 2차원(2D) 반도체 발굴은 재료 과학의 핵심 과제입니다. 최근 층상 SiAs 단결정 합성 성공으로 소층 SiAs 기계적 박리가 가능함이 밝혀져 단일층 제작 가능성을 시사했습니다. 밀도범함수이론(DFT)을 사용해 SiAs와 SiAs₂ 단일층이 동역학적으로나 열역학적으로 안정한 반도체이며, 각각 2.39 eV와 2.07 eV의 간접 밴드 갭을 가짐을 입증했습니다. 쌍축 변형 적용 시 이 간접 갭이 직접 갭으로 전환됩니다(SiAs 1.75 eV, SiAs₂ 1.60 eV). 더 큰 변형에서는 금속 전이가 유도됩니다. 캐리어 이동도 계산 결과 SiAs 정공 이동도가 최대 3.9 × 10³ cm² V⁻¹ s⁻¹로 MoS₂를 능가하며, 흑린과 유사한 이방성을 보입니다. 이 발견은 SiAs와 SiAs₂가 차세대 광전자 소재로서 유망함을 보여줍니다.

배경

2D 소재는 현대 응집물질 연구의 초석이 되어 전자적, 광학적, 기계적 특성을 조절할 수 있게 해줍니다. 그래핀 발견 후 실리신, 질화붕소, 전이금속 칼코게나이드(TMDs), 흑린, 보로펜, 아르세넨, 텔루렌 등 다양한 2D 결정군이 실험적으로 실현되었으며 각각 고유한 장점과 과제를 안고 있습니다[1-23].

그래핀은 탁월한 캐리어 이동도를 보이지만 영 밴드 갭으로 디지털 응용이 제한됩니다. TMDs는 고유 밴드 갭을 가지나 낮은 이동도를 보입니다[26-28]. 흑린은 밴드 갭과 높은 이동도를 겸비했으나 공기 중 화학적으로 불안정합니다[13,29]. 최근 SiAs와 SiAs₂ 층상 결정 성장 성공으로 안정적이고 박리 가능한 2D 반도체의 새로운 길이 열렸으며 잠재적으로 우수한 특성이 기대됩니다[30-32].

본 연구에서는 제1원리 DFT로 SiAs와 SiAs₂ 단일층의 구조적, 전자적, 수송 특성을 예측했습니다. 기계적 안정성을 확인하고 밴드 구조를 계산했으며, 변형 유도 밴드 갭 엔지니어링과 캐리어 이동도를 평가해 이 소재들을 고성능 광전자 소자 유력 후보로 자리매김했습니다.

계산 방법

모든 계산은 Vienna Ab-initio Simulation Package(VASP)로 수행했으며, 일반화 기울기 근사(GGA) 내 Perdew-Burke-Ernzerhof(PBE) 교환-상관 범함수를 사용했습니다. 투영 증강파(PAW) 방법으로 코어-가전자 상호작용을 기술했습니다. 20 Å 진공 간격으로 층간 결합을 방지했고, 500 eV 평면파 컷오프로 수렴을 보장했습니다. 브릴루앙 존 샘플링은 2D 시트에 15 × 5 × 1 Monkhorst-Pack 그리드를 사용했습니다. 구조 완화는 총 에너지 10⁻⁴ eV, 힘 0.01 eV/Å까지 수렴시켰습니다.

포논 분산은 PHONOPY의 초셀 접근법으로, 힘 상수는 VASP의 밀도범함수 섭동 이론(DFPT)에서 도출했습니다. 열적 안정성 확인을 위해 3×3×1 초셀(시간 단계 2 fs)로 300 K에서 6 ps NVT 분자동역학 시뮬레이션을 수행했습니다. 라만 스펙트럼은 CASTEP[39-41]으로 PBE 수준에서 계산했습니다.

결과 및 토의

SiAs와 SiAs₂ 단일층의 완화된 기하 구조는 그림 1에 나와 있습니다. 벌크 SiAs(C2/m 대칭)와 SiAs₂(Pbam 대칭)는 약한 반데르발스 힘으로 결합하며 층간 간격은 각각 3.06 Å, 1.66 Å입니다. SiAs 단일층은 마름모꼴 격자(a₁ = 3.69 Å, b₁ = 10.83 Å, φ = 99.81°)를 가지며 단위 세포당 Si 6개, As 6개 원자를 포함합니다. Si-Si 결합 2.35 Å, Si-As 결합 2.39-2.43 Å, 버클링 높이 4.86 Å입니다. SiAs₂는 직사각형 격자(a₂ = 3.68 Å, b₂ = 10.57 Å)를 가지며 단위 세포당 Si 4개, As 8개 원자를 포함합니다. As-As 결합 2.50 Å, Si-As 결합 2.41-2.45 Å, 버클링 높이 5.09 Å입니다. 전자 밀도 맵에서 As 원자가 더 높은 전기음성도로 Si를 희생시키며 전자 밀도를 축적해 강한 공유 결합을 확인했습니다.

라만 스펙트럼(그림 S2)은 단일층과 벌크 상 간 명확한 이동을 보이며, 이는 층간 반데르발스 상호작용 소실에 기인합니다.

안정성은 결착 에너지로 정량화했습니다: SiAs 5.13 eV/atom, SiAs₂ 4.98 eV/atom으로 아르세넨(2.99 eV/atom)과 실리신(3.71 eV/atom)보다 높습니다[18,43]. 포논 계산(그림 2a)은 Γ 근처 연화된 횡음향 분기를 제외하고 허수 모드가 없어 동역학적 안정성을 확인했습니다. 300 K 분자동역학(그림 2b)은 최소 에너지 변동과 구조 보존으로 열적 견고성을 입증했습니다.

밴드 구조 분석(그림 3)에서 두 단일층 모두 간접 반도체입니다: SiAs(E_g = 1.72 eV) VBM은 Y, CBM은 Γ; SiAs₂(E_g = 1.42 eV) 거의 직접 갭으로 VBM은 Y 근처, CBM은 약간 치우쳐 있습니다. HSE06 하이브리드 범함수 계산으로 갭이 SiAs 2.39 eV, SiAs₂ 2.07 eV로 증가했습니다. 오비탈 분해 결과 VBM은 p_y 오비탈 주도, CBM은 s 오비탈 주도로 변형 유도 이동이 간접-직접 전이를 유도할 수 있음을 시사합니다.

캐리어 이동도는 변형 퍼텐셜 이론으로 평가했습니다. 유효 질량(m*)은 이방적으로 변합니다: SiAs에서 m*_e = 0.15 m₀(x), 0.80 m₀(y); m*_h = 0.86 m₀(x), 0.22 m₀(y). SiAs₂에서 m*_e = 0.14 m₀(x), 2.05 m₀(y); m*_h = 0.65 m₀(x), 1.82 m₀(y). 탄성 상수와 변형 퍼텐셜(표 S1)로 전자 이동도: SiAs 0.66×10³ cm² V⁻¹ s⁻¹(x), 0.54×10³ cm² V⁻¹ s⁻¹(y); SiAs₂ 0.26×10³ cm² V⁻¹ s⁻¹(x), 0.11×10³ cm² V⁻¹ s⁻¹(y). 정공 이동도는 SiAs에서 3.90×10³ cm² V⁻¹ s⁻¹(x), 0.30×10³ cm² V⁻¹ s⁻¹(y)에 달해 MoS₂를 능가하며 흑린과 유사한 이방성 수송을 보입니다.

상태 밀도(PDOS) 플롯(그림 4)은 Si와 As 간 강한 s-p 혼성화를 보여주며, 국소화된 p_z 고립 전자쌍 상태가 버클링에 기여합니다. 전하 밀도 맵은 VBM이 Si/As p 오비탈 혼성, CBM은 주로 s 오비탈에서 유래함을 확인합니다.

변형 공학이 전자적 특성을 극적으로 조절합니다. 쌍축 압축 변형은 버클링 높이를 증가시키고 SiAs에서 일련의 간접-직접 전이를 유도합니다(그림 5c). 더 큰 변형에서는 금속 전이가 나타납니다. SiAs₂는 8-10% 인장 변형에서 직접 갭 반도체(E_g = 1.60 eV, HSE)로 전환되며 더 높은 변형에서 금속 상태가 됩니다(그림 5d). 변형에 따른 밴드 구조 진화는 그림 6, 7에 상세히 나와 있습니다. 대표적 직접 갭 구성은 그림 S7에 표시했습니다.

결론

제1원리 DFT 예측에 따르면 SiAs와 SiAs₂ 단일층은 동역학·열역학적으로 안정한 반도체입니다. 간접 밴드 갭(2.39 eV, 2.07 eV)은 변형으로 고도로 조절 가능해 간접-직접 전이와 금속 전이를 유도합니다. 캐리어 이동도는 MoS₂를 능가하며 흑린과 유사한 강한 이방성을 보입니다. 이 특성들은 SiAs와 SiAs₂를 나노스케일 고성능 광전자 소자용 매력적인 플랫폼으로 자리매김시킵니다.

약어

2D

2차원

CASTEP

Cambridge sequential total energy package

CBM

전도 밴드 최소값

DFT

밀도범함수이론

DFPT

밀도범함수 섭동 이론

DP

변형 퍼텐셜

GGA

일반화 기울기 근사

MD

분자동역학

NVT

몰-부피-온도 앙상블

PAW

투영 증강파

PBE

Perdew-Burke-Ernzerhof

PDOS

부분 상태 밀도

TMDs

전이금속 칼코게나이드

VASP

Vienna ab initio simulation package

VBM

원자가 밴드 최대값