歪み調整可能なバンドギャップと優れたキャリア移動度を持つSiAsおよびSiAs₂単層：第一原理的知見

概要

サイズableなバンドギャップと高キャリア移動度を兼ね備えた原子層厚二次元（2D）半導体の発見は材料科学の主要課題である。最近のSiAs単結晶合成により、少層SiAsが機械的剥離可能であることが示され、単層作製の実現可能性が示唆された。密度汎関数理論（DFT）を用いて、SiAsとSiAs₂単層が動力学的・熱力学的に安定な半導体であり、それぞれ2.39eVと2.07eVの間接バンドギャップを持つことを実証。二軸歪み印加でこれら間接ギャップが直接ギャップ（SiAs 1.75eV、SiAs₂ 1.60eV）に転換し、更大歪みで金属転移誘起。キャリア移動度計算でSiAs正孔移動度最大3.9×10³cm²V⁻¹s⁻¹を示し、MoS₂を上回り、黒リンと類似の異方性を示す。これら知見はSiAsとSiAs₂を次世代光電子デバイス有望候補として位置づける。

背景

二次元材料は現代凝縮系研究の基石となり、調整可能な電子的・光学的・機械的特性を提供。グラフェン発見後、シリセン、窒化ホウ素、遷移金属ダイカルコゲナイド（TMDs）、黒リン、ボロフェン、アルセニン、テルルレン等多様な2D結晶ファミリーが実験的に実現され、各々独自利点と課題を持つ[1-23]。

グラフェンは卓越キャリア移動度持つがゼロバンドギャップでデジタル応用制限。TMDsは固有ギャップ持つが低移動度[26-28]、ホスホレンはバンドギャップと高移動度兼備但空気中化学不安定[13,29]。最近のSiAsとSiAs₂層状結晶成長成功は、安定で剥離可能な2D半導体への新経路を開き、潜在的優位特性を秘める[30-32]。

本研究では第一原理DFTでSiAsとSiAs₂単層の構造的・電子的・輸送特性を予測。機械的安定性確認、バンド構造計算、歪み誘起バンドギャップエンジニアリングとキャリア移動度評価を行い、これら材料を高性能光電子デバイス有力候補として位置づけた。

計算手法

全計算はVienna Ab-initio Simulation Package（VASP）で、一般化勾配近似（GGA）内Perdew-Burke-Ernzerhof（PBE）交換相関汎関数使用。プロジェクタ拡張波（PAW）法でコア-価電子相互作用記述。20Å真空ギャップで層間結合防止、500eV平面波カットオフで収束確保。ブリルアンゾンサンプリングに15×5×1 Monkhorst-Packグリッド使用。構造緩和は総エネルギー10⁻⁴eV、力0.01eV/Åで収束。

フォノン分散はPHONOPY超胞アプローチで、VASP密度汎関数摂動理論（DFPT）由来力定数使用。熱安定性確認のため、3×3×1超胞（タイムステップ2fs）で300K、6ps NVT分子動力学シミュレーション実施。ラマンスペクトルはCASTEPでPBEレベル計算[39-41]。

結果と考察

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とSi-As結合長2.35Åと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.13eV/原子、SiAs₂ 4.98eV/原子、アルセニン（2.99eV/原子）とシリセン（3.71eV/原子）超える[18,43]。フォノン計算（図2a）はΓ点近く横音響分枝軟化除き虚数モードなし、動力学的安定性確認。300K分子動力学（図2b）は最小エネルギー変動と構造保持示し、熱的頑健性示唆。

バンド構造解析（図3）で両単層は間接半導体：SiAs（E_g＝1.72eV）VBM＝Y、CBM＝Γ；SiAs₂（E_g＝1.42eV）ほぼ直接、VBM Y付近、CBMわずかずらし。HSE06ハイブリッド汎関数計算でギャップ増大：SiAs 2.39eV、SiAs₂ 2.07eV。軌道分解でVBMはp_y軌道主導、CBMはs軌道主導、歪み誘起シフトで間接-直接転移駆動示唆。

キャリア移動度は変形ポテンシャル理論で評価。有効質量（m*）異方的変化：SiAsでm*_e＝0.15m₀（x）、0.80m₀（y）；m*_h＝0.86m₀（x）、0.22m₀（y）。SiAs₂でm*_e＝0.14m₀（x）、2.05m₀（y）；m*_h＝0.65m₀（x）、1.82m₀（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.60eV、HSE）転換、更高歪みで金属状態出現（図5d）。歪み下バンド構造進化詳細図6、7。代表的直接ギャップ構成図S7。

結論

第一原理DFT予測でSiAsとSiAs₂単層は動力学的・熱力学的安定半導体。間接バンドギャップ（2.39eV、2.07eV）は歪みで高度調整可能、間接-直接転移と金属転移可能。キャリア移動度はMoS₂超え、黒リン類似顕著異方性示す。これら特性はSiAsとSiAs₂をナノスケール高性能光電子デバイス有力プラットフォームとして位置づける。

略語

2D

二次元

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

価電子帯最大値