グラフェン仕事関数の制御：サイズ、官能基化、欠陥、ドーピングの第一原理計算

概要

電子デバイス応用に向け、第一原理密度汎関数理論を用いてグラフェンの仕事関数（WF）を体系的に設計・評価しました。無修飾、水酸基修飾、欠陥導入、ドーピングされたグラフェン構造を探索し、WFがリボン幅、表面官能基化、欠陥位置、ドーパント種に強く依存することを明らかにしました。無修飾ジグザググラフェンはアームチェア型より高いWFを示し、リボン幅拡大は一貫してWFを低下させます。水酸基はWFを上昇させ、その効果は基数とクラスタリングで増幅されます。欠陥は一般にWFを低下させ、中心部空孔の影響が最も顕著です。BやAlによる置換ドーピング（p型）はWFを大幅に上昇、NやP（n型）は急激に低下させ、Siは穏やかで安定した効果を示します。また形成エネルギーを算出しドーパント導入の容易さを評価：BとNは容易、Alは困難。これらの知見は、デバイス設計におけるグラフェン電子特性制御の堅実な理論的基盤を提供します。

背景

グラフェンの卓越した電気的、機械的、化学的特性は、センサ、電界効果トランジスタ（FET）、光起電力電極、ショットキーダイオード、真空管、LED接合など多様な応用へと導き、従来材料を凌駕する性能を発揮しています。しかし、デバイス性能はグラフェン仕事関数（WF）にかかっており、電荷注入、ショットキー障壁高さ、全体的なバンドアライメントを決定します。典型的な実験WF値は4.2–4.8 eVの範囲です。化学ドーピング、表面官能基化、欠陥工学、歪み、静電ゲーティング、レーザー照射など、WF調整戦略は多数報告されています。しかし、グラフェンのカイラリティ、サイズ、欠陥位置を横断的に比較した研究は限られ、ドーパントの熱力学的実現可能性が見落とされることも多々あります。本研究は、包括的な第一原理解析によりこれらのギャップを埋めます。

計算手法

全計算をCASTEPコードと密度汎関数理論（DFT）で実施。交換相関には一般化勾配近似（GGA）を採用し、不均一電子密度の優れた取り扱いを確保。真空層を15 Åに設定しスラブ間相互作用を排除、340 eVカットオフのウルトラソフト擬ポテンシャルを使用。9×9×1モンクホルスト–パックk点メッシュとメトフェッセル–パクソン緩和（0.05 eV）で収束を確保。仕事関数は真空準位とフェルミエネルギーの差（WF = V₀ – E_f）から取得。

結果と考察

リボン幅とカイラリティのWFへの影響

幅1–7単位格子のジグザグおよびアームチェアグラフェンナノリボン（GNR）をモデル化。両カイラリティで幅拡大に伴いWFが単調に低下；幅広リボンはバルクグラフェンWF（~4.48 eV）に収束。ジグザグGNRは一貫してアームチェア型より高いWFを示し、これは異なる結晶対称性に起因。バンドギャップ解析（表1）も幅拡大に伴うギャップ縮小を示し、WF傾向を支持。

水酸基化、欠陥、およびそれらの位置の影響

水酸基官能基化はWFを上昇：単一OHで4.50 eV、複数OHで5.10 eVまで到達。対称的・クラスタ化したOH配置が最高WFを与え、分散配置ではやや低下（4.66–4.83 eV）。欠陥は単一炭素空孔としてモデル化し、一般にWFを低下。中心空孔で最低WF（4.34 eV）、中心外れ空孔では無修飾GNR比でわずかに上昇（4.36 eV）。これらは欠陥誘起WF低下の先行実験観測と整合。

B、N、Al、Si、Pドーピングの効果

ドーピング濃度を2.4–14.6%（4×4超格子中1–6原子）で変化。p型ドーパントBとAlはWF上昇、Bは14.6%で5.15 eV達成。n型ドーパントNとPはWFを劇的に低下、Pは最高濃度で3.23 eV。Siドーピングは穏やかで安定したWF効果。形成エネルギー計算から、BとNは熱力学的に有利、Al、Si、Pは高い形成エネルギーを示し導入が困難。

結論

本系統的第一原理研究により、グラフェンWFはリボン幅、水酸基官能基化、欠陥工学、選択的ドーピングにより精密に調整可能であることを実証。これらの知見は、次世代デバイスアーキテクチャにおけるグラフェン電子特性最適化の明確な指針を提供します。