ペンタグラフェン：優れた感度を持つ次世代NOxガスセンサー

概要

グラフェンを凌駕する独自の特性を持つ二次元（2D）ペンタグラフェン（PG）は、ナノエレクトロニクスへの有望な応用により広範な注目を集めています。本研究では、単層PGとCO、CO₂、NH₃、NO、NO₂などの代表的な小さなガス分子との電子的・輸送特性を、第一原理計算と非平衡グリーン関数（NEGF）法を用いて調べ、この単層のセンシング能力を探索しました。吸着分子の最適位置とモードを決定し、吸着安定性における電荷移動の重要な役割と、化学結合形成が吸着系の電子構造に及ぼす影響を探りました。単層PGが、適切な吸着強度と明確な電荷移動を伴ってNOx（x = 1, 2）分子に最も選択的に吸着することを実証しました。さらに、NO₂（NO）吸着後のPGの電流-電圧（I-V）曲線では、電流が88%（90%）大幅に減少することを示しました。PGの優れたセンシング性能は、グラフェンやホスホレンなどの他の2D材料に匹敵、あるいは凌駕します。窒素酸化物へのこの超高感度と選択性により、PGは幅広い応用が期待される優れたガスセンサーとなります。

はじめに

グラフェンやホスホレンなどの単層または少層平面結晶で構成される二次元（2D）材料は、材料物理学の新たなパラダイムとして台頭し、独自の構造と物理化学的特性により注目が高まっています。これらは、大きな比表面積と完全に露出した活性サイトに関連しています。これらの特性は、ナノエレクトロニクス、センサー、触媒、太陽エネルギー変換デバイスなどの分野での幅広い潜在的応用を2D材料に与えています。

ペンタグラフェン（PG）は、カイロ五角形タイリングパターンに基づく炭素の新しい2D同素体で、平面シート幾何学において五角形のみ（sp²とsp³配位炭素原子の混合）で構成される個々の原子層を持ちます。ゼロバンドギャップを持つグラフェンとは異なり、その応用を大きく妨げているのに対し、PGは約3.25 eVの準直接固有バンドギャップを持ち、ドーピング、水素化、電界により調整可能です。その異常な原子構造により、PGは1000 Kまでの顕著なエネルギー的、動的、熱的、機械的安定性を持ちます。天然のバンドギャップと頑健な安定性により、PGはナノエレクトロニクス、センサー、触媒への高望ましい特性と大きな可能性を提供するかもしれません。一例として、PGベースの全炭素ヘテロ構造が静電ゲーティングや窒素ドーピングにより調整可能なショットキー障壁を示し、ナノエレクトロニクス応用の可能性を検証しています。興味深いことに、PG上での低温CO酸化のEley-Rideal機構のエネルギー障壁はわずか−0.65 eV（多くの貴金属触媒と同等）であり、BとB/Nドーピングでそれぞれ−0.11 eVと−0.35 eVまで低減可能で、PGが金属フリーかつ低コスト触媒となり得ることを説得力を持って示しています。最近の研究では、PGナノシートがNOの高選択的吸着を示し、ドーピングがH₂、CO、CO₂などのガス分子吸着を改善することも見出されています。グラフェンが優れたセンサー特性を理論的・実験的に示すように、ガス分子の吸着能力は、ガス分子吸着によりPGの電気抵抗率が影響を受けるため、PGもガスセンシング特性を持つことを示唆します。しかし、我々の知る限り、分子吸着がPGの電子特性に及ぼす影響に焦点を当てた先行報告はなく、PGの独自の電子特性を鑑みると、PGベースのガスセンサーの可能性を探ることは極めて望ましいです。

本研究では、密度汎関数理論（DFT）と非平衡グリーン関数（NEGF）計算を用いて、PG単層のガスセンサーとしての可能性を探りました。まず、PG単層へのCOx（x = 1, 2）、NH₃、NOx（x = 1, 2）など代表的分子の吸着挙動を調べました。PG単層へのNOxの優先的吸着（適切な吸着強度を伴う）は、気体NOxへのPGの高い選択性を示します。NO₂吸着前後での電流-電圧（I-V）関係の劇的な変化は、PGの優れた感度を示唆します。ガス分子に対する感度と選択性の両面から、PGは高性能センシング応用の有望な候補となります。

方法

構造緩和と電子計算は、Vienna Ab initio Simulation Package（VASP）に実装された第一原理計算に基づくDFTを用いて実施しました。交換相関相互作用は、Perdew-Burke-Ernzerhof（PBE）汎関数の一般化勾配近似（GGA）で処理しました。PGモデルはxy平面で周期的で、z方向に少なくとも15 Å分離されています。エネルギーカットオフは450 eV、3 × 3超格子に対して9 × 9 × 1 Monkhorst-Packグリッド（TRANSIESTA用は9 × 3 × 9）をブリルアンゾーン積分に用いました。より正確な吸着エネルギーを得るためDFT-D2法を用いました。力収束基準は0.03 eV/Å未満です。NOxは常磁性であるため、その吸着計算にはスピン分極を含めました。輸送特性は、TRANSIESTAパッケージに実装された非平衡グリーン関数（NEGF）法で研究しました。接触領域を通る電流はLandauer-Büttiker公式で計算：\( I(V_b)=G_0 \int_{\mu_L}^{\mu_R} T(E,V_b) dE \)、ここでG₀とTは量子コンダクタンス単位と、ポテンシャルバイアスV_b下でエネルギーEで入射する電子の透過率です。二つの電極間の電気化学ポテンシャル差はeV_b = μ_L − μ_Rです。

結果と考察

吸着系の構造的特性とエネルギーを調べる前に、まず単層PGの格子定数を最適化し、a = b = 3.63 Åを得ました。これは先行報告値と一致します。最も有利な配置を見つけるため、ガス分子吸着のために異なる吸着サイトと配向を調べ、それぞれを3 × 3超格子PG上に配置しました。完全緩和後、NOx分子がPG上で強い化学結合を介して化学吸着するのに対し、他の3分子（COx、NH₃）は物理吸着することがわかりました（図1）。CO、CO₂、NH₃分子はそれぞれ2.40、2.73、2.43 Åの吸着距離でPG上に留まり（表1）、弱いファンデルワールス相互作用を示します。対照的に、双極子NOx分子はC原子のトップ位置に引き付けられ、1.43〜1.56 Åの結合長でC原子と化学結合を形成します。PG/NO₂の場合、NとO原子の両方がPGのC原子と化学結合可能です（図1e）。

吸着構造。a–d CO、CO₂、NH₃、NO吸着時の完全緩和構造モデルの側面図（上）と上面図（下）。e NO₂吸着時の2つの結合モード側面図（結合エネルギーEaを記載）。aにガス分子とPG層間距離、d、eに界面のC-NおよびC-O結合長を示す（オングストローム単位）。簡略化のため、これらの構造モデルをa PG/CO、b PG/CO₂、c PG/NH₃、d PG/NO、e PG/NO₂と略記

分子のPG上での安定性を吸着エネルギー（E_a）で評価しました。E_a = E_pg+gas − E_gas − E_pg、ここでE_pg+gas、E_pg、E_gasはそれぞれガス吸着PG、プリスタインPG、孤立分子の全エネルギーです。表2によると、グラフェンとホスホレンのガスセンサーとしての潜在的利用と同様に、PG/NOとPG/NO₂の吸着エネルギーはそれぞれ分子あたり−0.44 eVと−0.75 eV（ガス捕獲の基準とされる−0.5 eVに接近）で、室温での熱擾乱（k_B Tスケール）に耐える十分な大きさです。一方、PG/COxとPG/NH₃の吸着エネルギーは小さい（−0.05〜約−0.11 eV）、COxとNH₃分子がPGに容易に吸着しないことを示します。これらの結果は、単層PGが有毒NOxガスに高い選択性を持つことを検証します。より重要なのは、PGのNOxセンシング特性が、グラフェン、シリセン、ゲルマネン、ホスホレン、MoS₂など他の2Dナノシートとは異なり、それらがNOxと/or COx（NH₃）を区別できないのに対し、独自のものであることです（表2）。

ガス吸着が、吸着エネルギー決定とホスト層抵抗変化における電荷移動の重要な役割を果たすことは、ほとんどの場合で実証されています。まず、3D電荷密度差Δρ = ρ_tot(r) − ρ_pg(r) − ρ_gas(r)で界面電荷移動を計算し、直観的に可視化しました。ここでρ_tot(r)、ρ_pg(r)、ρ_gas(r)は、ガス吸착あり/なしPGと同一配置の自由ガス分子の電荷密度です。図2に、NOx、COx、NH₃のPG吸着時の計算電子移動を示します。明らかに、界面での電荷密度変化は顕著です。化学吸着NOx系と比較して、PG/COとPG/CO₂界面での電荷再分配は相対的に弱いです。これは共有結合相互作用がファンデルワールス力より強いためです。NH₃吸着では、NH₃分子周辺で電荷再分配が起こります。

電荷密度差プロット。図1とは異なる順序で、各ケースの吸着構造と電荷移動をa–eにプロット。黄色等値面は電子獲得、青は電子損失を示す。等値面値単位はe Å⁻³。共有結合a PG/NOとb PG/NO₂構造での電子移動が他より明らか

Bader法に基づく更なる電荷解析で、これら系での電荷再分配の定量的尺度を得ました（表1）。予想通り、COxとNH₃のPG物理吸着では、PGとガス分子間で少量（< 0.025 e）の電荷移動のみが起こり、弱い結合をさらに裏付けます。対照的に、化学吸着系では電荷移動量が10倍以上多く：PG層からNO₂（NO）分子へ最大0.517 e（0.243 e）が移動します（表1）、これらの強い吸着エネルギーと一致します。吸着強度と電荷移動のこの系統的傾向は、PG上でのガス分子吸着メカニズム理解を助け、また電界によりガス吸着を制御可能であることを示唆します（単層MoS₂上でのNOx吸着の場合と類似）。

次に、ガス吸着がPG電子特性に及ぼす影響を調べました。図3に、ガス分子吸着あり/なしPGの全状態密度（DOS）と、対応する個々からの投影DOSを示します。プリスタインPGから2.10 eVのバンドギャップを得ました。これは純粋PGの先行DFT結果と一致し、PBE/GGA汎関数が半導体バンドギャップを過小評価するためです。これは閾値バイアス（観測可能な電流を生じる電圧）に影響しますが、以下で実証される他の輸送特性には影響しないと期待されます。図3aはプリスタインPGのDOS、図3b、cはCOx吸着でPGの価電子帯（VB）または伝導帯（CB）近傍のDOSが明確には影響されず、これらの小さな吸着エネルギーと弱い電荷再分配と完全に整合します。NH₃吸着でVBトップ近傍に小さな状態が生じます（図3d）が、分子の物理吸着はフェルミ準位近傍のDOSに顕著な変化を与えません。これらの結果は、COxとNH₃吸着がPG電子構造に有意な影響を与えないことを示します。対照的に、NOx吸着PGシートではフェルミ準位近傍に明確な混成状態が観察されます（図3e、f）。この特徴と主要な電荷密度再分配が、NOxとPG単層の強い相互作用を示し、顕著なバンド構造変化をもたらします。これはPG輸送特性に大きな影響を与え、非常に感度の高いガスセンサーとなります。

全電子状態密度。a プリスタインペンタグラフェンDOS。b–f 各ガス分子吸着時のPG全DOS（青線）とガス分子からの部分DOS（赤線）。フェルミ準位をゼロとし、黒破線で表示

一部2D材料がガス分子吸着に極めて敏感で、極めて低密度ガス分子に対応することが実証されています。PGのガス濃度依存感度をシミュレートするため、吸着ガス被覆率がPG特性に及ぼす影響を計算しました。PG/NO系を例に、被覆率5.56%で吸着エネルギーは分子あたり約−0.44 eV。被覆率が3.13〜2.0%へ減少すると、吸着エネルギーは分子あたり約−0.32 eVへ減少します。これはガス濃度変化が主な結論を変えないことを示します。したがって、以下の計算では、5.56%被覆率のPG/NO系モデル（3 × 3超格子使用）を代表として電子的・輸送特性計算に用います。

PG単層のNOxモニタリング感度を定性評価するため、二探針モデルを用いてNOx吸着前後の輸送透過率と電流-電圧（I-V）関係をNEGF法でシミュレートしました（図4a）。物理像を明確にし計算負担を軽減するため、二探針系（疑似「デバイス」構造）を用いました。ここでは「偽電極」がクリーンナノシートの周期的延長から構築され、先行研究で広く用いられています。ここでは、ガス吸着あり/なしの3 × 3 PG超格子を、左右電極と中心散乱領域それぞれに用いました（図4a）。比較のため、ガス吸着なし中心散乱領域でも同一計算を実施。PGのNOx吸着あり/なしI-V曲線を図4b1、4c1に示します。常磁性分子NOxのPG吸着はスピン分極を誘起し、スピン分極電流を生じます。バイアス電圧印加時、左電極フェルミ準位が右電極に対して上方シフトします。したがって、左電極VB最大が右電極CB最小に達した後のみ電流が流れます。結果として、バイアス電圧が3.25 V未満では中心散乱領域を電流が通過しません。これはPGの固有ギャップ（約3.25 eV）に近いです。バイアス電圧が3.25 Vから増加すると、両スピンチャネルの電流が急増します。3.9 Vバイアス下で、ガス吸着なしPGを通過する電流は13.4 μAですが、PGがNO₂分子を吸着すると同バイアス下で電流が急減し1.6 μAとなり、約88%減少します。さらに、PGがNO分子を吸着すると電流が1.34 μAへ減少し、約90%減少です。被覆率効果を探るため、4 × 4と5 × 5超格子へ1分子吸着させた場合を補足図S1でさらに検討。分子とPGシート間相互作用は被覆率で大きく変わらず、類似の吸着エネルギーE_aになります。5 × 5超格子中心領域のPG/NO輸送特性を計算し補足図S2に示します。3.9 Vバイアス下で、1 NO分子吸着の5 × 5超格子中心領域を通過する電流は2.87 μAへ減少（約79%減少）。電流の劇的減少はNOx吸着後の抵抗大幅増加を示し、実験で直接測定可能です。電流の有意変化は、PGセンサーのNOxへの超高感度を意味し、シリセンやホスホレンなど他の2Dナノシートに匹敵または凌駕します（表2で明確に表示）。

二探針系の模式図（a）で、半無限左/右電極領域（赤陰影領域）が中心散乱領域と接触。電極と散乱領域に、それぞれガス吸着あり/なし3 × 3超格子を使用。b1、c1に純粋PGとNO、NO₂吸着PGのI-V曲線を表示。零バイアス下透過スペクトルをc1、c2に表示

NOx吸着PG抵抗増加メカニズムを解明するため、零バイアス下NO₂吸着PG透過スペクトルを計算し図4cに表示。フェルミ準位周辺に幅3.25 Vの零透過領域が観察され、この外側に山岳状特性が透過スペクトルに現れます。DOS（図3f）と同傾向で、PBE汎関数選択が電子構造と輸送特性に大きな影響を与えないことを証明。図3fで、最低空分子軌道（LUMO）状態と最高占有分子軌道（HOMO）状態がギャップエッジに位置し、主にp_z軌道で形成されることを示します。電荷がC p_z軌道からNO₂分子へ移動すると、LUMOとHOMO状態がNO₂吸着により明確に影響を受けます。これは吸着NO₂分子が電荷キャリアの強い散乱中心となり、NO₂分子により誘起されるゾーン中心付近の局在状態により移動度が低下することを示します。つまり、伝導チャネルが阻害され、キャリア寿命または平均自由行程が短くなり、NOx吸着PGで移動度が小さくなります。

ガスセンサーの重要因子の一つである回復時間は、センサーが元の抵抗の80%に戻る時間です。遷移状態理論により、回復時間τはτ = ω⁻¹ exp(E*/K_B T)で計算可能。ここでωは試行周波数（先行報告により約10¹³ s⁻¹）、Tは温度、K_Bはボルツマン定数（8.318 × 10⁻³ kJ/(mol·K)）、K_B Tは室温で約0.026 eV、E*は脱着エネルギー障壁です。回復時間は脱着障壁と密接に関連：同温度で脱着障壁が低いほど、PG表面でのNOx回復時間が短くなります。脱着は吸着の逆プロセスと見なせるため、E_adをポテンシャル障壁(E*)と仮定するのは妥当です。したがって、PG/NOとPG/NO₂のポテンシャル障壁(E*)はそれぞれ0.44 eVと0.75 eV。300 Kでの両系計算回復時間はそれぞれ2.24 × 10⁻⁶ sと0.34 sで、PGセンサーが初期状態へ完全回復可能であることを示します。以上より、PGは高感度と迅速回復時間を持つNOxガス向け潜在材料と結論付けられます。

結論

本研究では、DFT計算を用いてPG単層と代表的ガス分子吸着時の構造的、電子的、輸送特性を体系的に調査しました。結果は、PG単層がシリセンやホスホレンなど他の2D材料と比較して、適切な吸着強度を持つ有毒NOxガスに最も選択的な単層の一つであることを示しました。PG電子抵抗はNO₂吸着で劇的に増加し、超高感度を意味します。一言で言えば、PGは高感度と迅速回復時間を持つNOxガスに優れたセンシング性能を発揮します。このような独自の特徴は、単層PGを優れたガスセンサーとして望ましい候補にします。

データおよび資料の入手可能性

本論文の結論を支持するデータセットは本論文内に含まれ、データと材料に関する更なる情報は、責任著者への動機付けられた要請により関心ある当事者へ提供可能です。

略語

2D:

二次元

CB:

伝導帯

DFT:

密度汎関数理論

GGA:

一般化勾配近似

HOMO:

最高占有分子軌道

LUMO:

最低空分子軌道

NEGF:

非平衡グリーン関数

PBE:

Perdew-Burke-Ernzerhof

PG:

ペンタグラフェン

PG/CO:

CO吸着ペンタグラフェン

PG/CO₂:

CO₂吸着ペンタグラフェン

PG/NH₃:

NH₃吸着ペンタグラフェン

PG/NO:

NO吸着ペンタグラフェン

PG/NO₂:

NO₂吸着ペンタグラフェン

VASP:

Vienna Ab initio Simulation Package

VB:

価電子帯