ZnO修飾C₃NナノシートにおけるSF₆分解生成物吸着の理論的知見

概要

優れた絶縁媒質であるSF₆は、高電圧絶縁装置に広く用いられ、電力システムの安全運転を支えています。しかし、長期間運用される装置では避けられない部分放電によりSF₆が分解し、絶縁性能が劣化します。本研究では、第一原理計算（DFT）を用いて、ZnO修飾C₃N（ZnO-C₃N）ナノシートにおけるSF₆分解生成物の吸着およびセンシング特性を調べ、SF₆絶縁装置の運転状態評価に向けた新規ナノ材料候補の提案を目指しました。まずZnO-C₃N単層の構造を検討し、続いて代表的な4種のSF₆分解生成物（H₂S、SO₂、SOF₂、SO₂F₂）に対する吸着挙動をシミュレーションしました。その結果、ZnO-C₃N単層はSO₂、SOF₂、SO₂F₂に対して望ましい反応性と感度を示し、これらのガス分子に激しい変形と大きな吸着エネルギーをもたらすことが判明しました。この知見は、ZnO-C₃N単層ベースのガス吸着材がSF₆絶縁装置から不純物ガスを除去する用途に応用可能であることを示唆しています。解析によれば、ZnO-C₃N単層は絶縁強度の維持と電力システムの安全運転確保に資する材料として有望です。

はじめに

ナノテクノロジーの急速な発展に伴い、新規ナノ材料ベースのセンサ応用が近年増加しています。迅速な応答、低消費電力、低コスト、高感度といった利点から、ナノセンサは医療、軍事、ガス検知、環境モニタリングなどの分野で盛んに研究されています[1-4]。中でも抵抗変化型センサは、構造と作動機構がシンプルなため広く採用されています。初期にはグラフェンが、大きな比表面積、高いキャリア移動度、優れた熱伝導性など優れた物理化学的特性によりガス検知材料として注目されました[5-8]。しかし、ゼロバンドギャップ特性のため[9,10]、CO、CO₂、CH₄、N₂、NO₂、NH₃、H₂、H₂Oなどの一般的ガス識別で性能が劣ります[11]。その後、研究者たちの努力により、グラフェンと同等の特性を持ちながらゼロバンドギャップを持たない多数の新規ナノ材料がガスセンシング分野に登場しました。遷移金属ダイカルコゲナイド（TMDs）[12-14]、金属炭化物/窒化物[15]、層状III-VI族半導体[16,17]、III-V族窒化物[18-20]などが代表例です。

新たに現れたグラフェン類似材料の中でも、C₃Nはヘキサアミノベンゼン三塩酸塩単結晶の直接熱分解や2,3-ジアミノフェナジンの重合により合成され[21,22]、ガス吸着材として大きな注目を集めています[23-25]。C₃Nは本質的に0.39 eVの間接バンドギャップを持つ半導体であり、異なる直径の量子ドット作製によりバンドギャップを制御可能です[22]。微視的には、C₃Nは2×2グラフェン超格子の2つの炭素原子が窒素で置換された構造と見なせ、6つの炭素と2つの窒素からなる平面ハニカム格子を持ちます。添加されたN原子により、固有C₃Nはグラフェンと同等の構造安定性を保ちながら、より強い化学活性と高いキャリア移動度を示し、ガス検知候補として競争力があります[26]。吸着能力の面では、固有C₃NがNO₂とSO₂検知で優れた選択性と感度を持つ一方、他のガスでは明確な吸着効果がないと報告されています[27]。しかし、不純物粒子による修飾でC₃Nの表面反応性を大幅に向上できることが明らかになっています。例えば、PashangpourとPeyghan[28]は、固有およびドープC₃NナノシートのCO吸着能を比較し、Alドーパントが固有C₃Nよりはるかに強い結合相互作用をもたらすことを示しました。その後、Zargham Bagheri[29]はSiドープC₃Nのアセトン吸着を理論研究し、C原子をSiで置換することで吸着エネルギーが-9.7から-67.4 kcal/molに増大し、より多くのC原子を置換すると感度が向上することを発見しました。

金属酸化物は、ガス相互作用のための化学反応性を高める表面修飾材として一般的です。金属酸化物半導体の一つであるZnOは、約3.37 eVのバンドギャップ、約60 meVの励起結合エネルギー、化学的安定性、優れた光触媒特性、特定ガスへの高い活性を持ちます[30]。ZnOは0次元（0D）、1次元（1D）、2次元（2D）のナノ構造形態で成長可能で、それぞれナノクラスタ、ナノワイヤ/ナノチューブ、ナノシート/ナノリボンが例として挙げられます[31]。サイズと形態が容易に制御できる特性[32,33]から、ZnOナノ粒子はナノ表面のセンシング性能向上ドーパントとして有望視されています[34-36]。最近、単一分子金属酸化物ドーパントによるナノ材料表面活性向上の理論研究がいくつか提案されています。E. Mohammadi-Maneshら[37]はCuおよびCuO修飾グラフェンのH₂S吸着能を理論的に調べ、H₂S吸着後に修飾グラフェンの伝導度が固有グラフェンと比較して有意に変化することを発見しました。AsadiとVaezzadeh[38]は、BおよびCuO修飾グラフェンシートを設計し、DFTに基づきH₂SとCO検知を研究しました。これらの研究ではDFTを用い、定義した超格子に基づき原子/分子構造を周期的に拡張し、系全体の物理特性を計算しました。この手法では、センサとしての基板上への原子/分子吸着を用いてナノ構造を研究します。これらの報告に刺激を受け、本研究では単一金属酸化物分子ドーピングについて更なる関連研究を行い、ZnOナノ粒子の代わりにZnO分子でC₃Nをドープすることで、ZnOのガスセンシングへの影響を簡略化して探索しました。

SF₆は、優れた絶縁およびアーク消弧特性を持つため、ガス絶縁開閉装置（GIS）の媒質として広く用いられています[39]。GIS内部の避けられない欠陥、例えば部分放電（PD）によりSF₆が分解し、SF₄、SF₃、SF₂などの低フッ化硫化物が生成されます[40]。これらの副生成物は微量の水分や酸素とさらに反応し、H₂S、SO₂、SOF₂、SO₂F₂などの安定化学物質を生成します[41]。これらの副生成物の絶縁信頼性はSF₆よりはるかに低く、放置するとPD進行を加速します。したがって、GISの安全運転を保証するため、これらのガスを検知または除去することが不可欠です。本論文では、ZnOをドーパントとして選び、ZnO修飾C₃N（ZnO-C₃N）単層モデルを構築し、典型的なSF₆分解生成物（H₂S、SO₂、SOF₂、SO₂F₂）に対する吸着性能を理論的に研究しました。構造変化、電子移動挙動、バンド構造、状態密度（DOS）を分析することで、ZnOドーパントがC₃N表面とガス分子の相互作用に与える影響を包括的に調べました。本研究の目的は、ZnO-C₃N単層の詳細な吸着およびセンシング機構を明らかにし、SF₆絶縁装置内の不純物ガス検知または除去への潜在的応用に資することです。

計算詳細

全てのDFT計算はDmol³パッケージ[42]で実施しました。現実系の非均一電子密度をより良く記述するため、一般化勾配近似（GGA[43]）におけるPerdew-Burke-Ernzerhof（PBE）汎関数とTS分散補正[44,45]を用いて電子交換相関項を扱いました。コア処理にはDFT半コア擬ポテンシャル（DSSP）を採用し、原子軌道基底セットには分極機能付き二重数値基底（DNP）を選択しました[46]。モンクホルスト-パックk点メッシュ6×6×1を幾何学最適化と電子構造計算の両方で用いました[47]。幾何学最適化におけるエネルギー収束許容値、最大力、最大変位はそれぞれ1.0×10⁻⁵ Ha、0.002 Ha/Å、0.005 Åに設定しました[48]。また、単層とその周期像の間に15 Åの真空層を設け、それらの相互作用を回避しました。

図1に示すように、C₃N単層の2×2×1超格子とガス分子を構築し、ドーピングと吸着プロセスの前に最適化しました。最適化単層におけるC-N結合長（1.422 Å）は、N原子の半径がC原子より大きいため、C-C結合長（1.418 Å）よりわずかに長くなります。本研究で得られた格子定数は4.92 Åで、文献[25,49]の報告値と類似しています。分子と単層間の電荷移動はヒルシュフェルト解析で計算し、ガス分子の電荷変化をQ_Tで表しました。正のQ_Tはガス分子の電子供与挙動を、負の値は電子受容挙動を示します。

結果と考察

ZnO修飾C₃N単層の解析

幾何学最適化後、ZnOをC₃N単層表面に異なる向きと位置で配置し、ZnO-C₃Nの最も合理的な構成を探索しました。図2aによると、ZnO粒子は六員環構造の中心（P_H1, P_H2）、C-CおよびC-N結合の中点（P_B1, P_B2）、C原子（P_C）およびN原子の真上において、面に垂直（O₁, O₂）および平行（O₃）な向きでC₃N単層に接近します。ZnO-C₃N単層の安定性を評価するために生成エネルギー（E_form）を定義し、以下のように計算しました：

$$ E_{\mathrm{form}} = E_{\mathrm{ZnO-C_3N}} - E_{\mathrm{ZnO}} - E_{\mathrm{C_3N}} \tag{1} $$

ここで、E_ZnOとE_C₃Nはドーピング前のZnO分子とC₃N単層のエネルギー、E_ZnO-C₃NはZnO-C₃N構造のエネルギーです。拡張原子構造において近接「結合」が生じると全エネルギーが低下し、負のE_formとなります[50]；最も大きなE_form（最も負の値）を示す構造を吸着および更なる解析用に選択しました。

全てのZnO-C₃N単層構成は補足図S1, S2, S3に示されています。これらの結果は、大きなE_formを持つ構造の多くがO₂向きであり、ZnO粒子がO原子を向けてC₃N表面に接近し、C原子によって捕捉されることを好むことを示しています。さらに、本研究の各修飾プロセスはE_formが負値であるため自発的であり、ZnO粒子をO₂向きでS_C位置に配置したときに最大のE_formが得られます。図2bに見られるように、二原子分子は40°傾いてC₃Nに付着しています。Zn-O結合は1.733 Åから2.182 Åに伸長し、C-O結合は1.338 Åと測定されました。ZnOの影響でC₃N表面はもはや平坦ではなく、ある程度の歪みが生じ、O原子に最も近いC原子が表面から引き出されます。ZnO-C₃N単層の電子的振る舞いをさらに議論するため、変形電荷密度（DCD）と状態密度（DOS）を図3に示します。図3aでは、赤領域が電荷密度の増加、青領域が減少を表します。ZnO分子が吸着すると、C₃N単層から0.255 eを引き抜き、O原子周辺に明確な赤領域が認められます。一方、Zn原子は青領域に囲まれ、O原子とZn原子の電気陰性度差を反映しています。さらに、O原子とC原子間の電荷密度の有意な上昇はC-O結合の形成を示唆し、これは図3cに示すO 2p軌道とC 2p軌道の状態間の強い混成によっても支持されます。図3bのDOS曲線から、ZnO導入により系のDOSが増加し、いくつかの新しいピークが出現することが明らかです。新たに現れた小さなピークは約-2.5 eVでO原子に由来し、-5.6 eVに位置する最大のピークは明らかにZn 3d軌道に起因します。DOSの変化と軌道混成の確認は、ZnO粒子がC₃N表面に強固に吸着し、系全体の電子構造に大きな影響を与える事実を裏付けています。

ZnO-C₃N単層の吸着挙動

各ガス分子をZnO-C₃N単層表面上の様々な向きで配置し、吸着パラメータを十分に比較して最も望ましい構成を選択しました。例えば、H₂SとSO₂のような三原子分子では、3原子で構成される平面を表面に対して平行または垂直にし、S原子を上向きまたは下向きにしました。吸着エネルギー（E_ads）を用いて異なる吸着構造のエネルギー変化を記述し、以下のように計算しました：

$$ E_{\mathrm{ads}} = E_{\mathrm{ZnO-C_3N/gas}} - E_{\mathrm{ZnO-C_3N}} - E_{\mathrm{gas}} \tag{2} $$

ここで、E_gasとE_ZnO-C₃Nは吸着前の孤立ガス分子とZnO-C₃N単層のエネルギー、E_ZnO-C₃N/gasはガス吸着後の系のエネルギーを表します。各状況の局所的最小全エネルギーを得た後、最大のE_adsを持つ構造のみを選択し、図4に示して更なる議論を行いました。また、電荷移動機構をより良く理解するため、電子密度差（EDD）を図5に描きました。

H₂S吸着系（図4a）では、H₂S分子が平行位置で吸着し、H₂S分子のH原子とZnOドーパントのO原子間の最近接原子距離は2.042 Åと測定されました。捕捉されたH原子のH-S結合は、孤立状態の1.357 Åと比較して1.374 Åに伸長しましたが、もう一方のH-S結合は吸着過程で変化しませんでした。幾何学構造のわずかな変形は、H₂S分子とZnO-C₃N単層の相互作用が弱いことを示唆します。E_ads（-0.444 eV）とQ_T（-0.077 e）を合わせると、H₂S分子がZnO-C₃N単層に安定して付着できないことが明らかであり、ZnO-C₃N単層はH₂S検知には不向きと判断されます。SO₂分子（図4b）では、両方のO原子がZn原子によって捕捉され、距離はそれぞれ2.020および2.031 Åでした。SO₂の存在によりZnO粒子は表面に近づき、Zn-O-Cで形成される角度は129°から118°に減少しました。同時にヒルシュフェルト解析から、S原子が0.164 eを失って電子供与体として振る舞い、O原子への吸着効果が弱まるため、S-O結合が1.481から1.619 Åに伸長することがわかりました[51]。逆に、ZnOのO原子の電子増加（0.292 e）による原子間力の強化で、Zn-O結合が2.182から1.869 Åに短縮しました。幾何学的変化と電子移動に加え、E_adsが-1.222 eVまで達することは、吸着過程における強い相互作用のさらなる証拠であり、ZnO-C₃N単層のSO₂検知への潜在的応用を示唆します。SOF₂およびSO₂F₂吸着系（図4c, d）では、標的分子が形態の劇的な変化を経験していることがわかります。SOF₂のF原子はS-F結合の拘束から解放され、1.830 Åの距離でZn原子に吸着されます。SO₂F₂ではさらに劇的な変化が観察され、両方のS-F結合が切断され、1.802および1.883 ÅのZn-F結合が形成されます。Zn-F結合形成に加え、SOF₂系では残ったSOF基もZn-O結合を通じてZn原子によって捕捉されます。一方、SO₂F₂系では興味深いことに、SO₂F₂分解で生成したSO₂基は捕捉されずZn原子から距離を保ち、そのS-O結合長は孤立SO₂と同一です。ガス分子の顕著な変形は、SOF₂系で-1.683 eV、SO₂F₂系で-3.497 eVという大きなE_adsと関連しています。これに基づき、SOF₂およびSO₂F₂のZnO-C₃Nへの吸着は、この過程で大量の電荷移動を伴う強い化学吸着と判定され、有意な電子軌道混成の存在が示唆されます。大きなE_adsとQ_Tにより、ZnO-C₃N単層はSO₂、SOF₂、SO₂F₂に対して他のナノ材料より安定した吸着を提供できます（表1参照）。ZnO-C₃N単層の吸着構造は、リストアップされたナノ材料より0.358-3.281 eVおよび0.038-0.811 e大きく、ガス検知用途におけるこの材料の吸着性能を保証します。しかし、大きなE_adsの結果として、これらのガス分子（SO₂、SOF₂、SO₂F₂）が強い相互作用力から逃れるのは困難と推測されます。したがって、センサポイズニングによる性能劣化を防ぐため、高温アニーリングや紫外線照射などの対策を講じてZnO-C₃N単層の脱着性能を改善すべきです。具体的な脱着性能解析は「ガスセンシング性能評価」節で提供します。

図5のEDDでは、青部分が電子蓄積領域、紫部分が電子枯渇領域を示します。H₂S系では、H原子とO原子間に小さな蓄積領域が見られるものの、大部分の蓄積と枯渇領域はH₂S分子周辺に位置し、H₂S分子内の少量の電荷移動と分子軌道の再分配を示唆します。SO₂吸着系では、S原子とZn原子周辺に明確な枯渇領域があり、蓄積領域は主にO原子周辺とS原子上方に分布します。この現象はSO₂分子の電子受容体としての役割を確認し、ヒルシュフェルト解析から得られたQ_T（-0.426 e）と一致します。SOF₂およびSO₂F₂系では、蓄積と枯渇領域の大規模さから、両系で顕著な電荷移動と電子混成が起きていることがわかります。蓄積領域は標的ガスの原子間に分布し、枯渇領域は主にZn原子周辺に局在化しており、Zn原子の明確な電子供与特性を反映しています。これらの電子的振る舞いは、ZnO-C₃N単層がガス分子に強い吸着を持つという仮定をより説得力あるものにします。

ガス吸着時のZnO-C₃N単層の電子的特性

グラフェン[56]、SWCNT[57]、MoSe₂[58]ベースのガスセンシング研究で報告されているように、DOSはガスとナノ構造間の電子的振る舞いを調べるもう一つの重要なパラメータです。図6aから、H₂Sの分子軌道再分配が図5aのEDDから導かれた結論と一致することがわかります。H 1sとO 2p軌道の混成は-4 eVおよび-6 eV付近で見られますが程度は低く、弱い相互作用と新しいH-O結合形成の微小な可能性を示します。SO₂系（図6c）では、反結合軌道がフェルミ準位付近にわずかに移動し、軌道の一部が分離から連結へ変化し、SO₂分子の電子構造の顕著な再分配を表します。原子間相互作用（図6d）では、O 2p、Zn 4s、Zn 3d軌道が-6、-4、-2 eVなど複数のエネルギー準位で混成していることがわかります。この混成はO原子とZn原子間の強い化学的相互作用を意味し、最適化構造で計算されたZn-O結合形成を効果的に支持します。SO₂FおよびSO₂F₂系では、構造の劇的な変形により分子軌道が強く活性化・再分配され、多くの新しい軌道が形成されます。SOF₂のF 2pとO 2p軌道がZn 4sおよび3d軌道と-8、-7、+3 eVで強く混成します。F原子とZn原子間の混成は-7、-5、+3.5 eV付近で確認できます。Zn原子と捕捉されたF、O原子間の顕著な混成は、安定な化学結合（Zn-FおよびZn-O結合）形成の証拠であり、ZnO-C₃N単層とこれら2ガス間の強い吸着効果の説明になります。4つの吸着系から得られた結果を合わせると、H₂S分子を除く他の3分子（SO₂、SO₂F、SO₂F₂）がZnO-C₃N単層に露出したときに強固に吸着されることがわかります。この結論は、本物質が高電圧装置でのガス除去用途の可能性を持つことを証明します。

ガスセンシング性能評価

抵抗変化型デバイスでのガス検知には、吸着後の適度な伝導度変化が必要です。系の伝導度は以下の式に従いバンドギャップに関連します[59]：

$$ \sigma = A \times e^{-E_g / 2k_B T} \tag{3} $$

ここでAは定数、k_Bはボルツマン定数（8.62×10⁻⁵ eV K⁻¹）、Tは温度です。伝導度とバンドギャップの間に反比例関係が認められ、バンドギャップが広いほど電子が禁制帯を越えるのが困難になります。図7aは、ZnO-C₃N単層のバンドギャップが0.168 eVと非常に小さく、固有C₃Nのバンドギャップ（0.39 eV）の半分以下であることを示しています。一方、伝導帯底（M点）と価電子帯頂（Γ点）の位置が異なることから、半導体特性と間接バンドギャップの性質は維持されています。吸着系のバンドギャップについては、図7b-dに異なる変化が見られます。H₂S系では、伝導帯底の下方移動によりバンドギャップが0.125 eVに減少します。他の系では、価電子帯頂付近に出現した新しい不純物準位がM点付近でフェルミ準位と接し、これらの系のバンドギャップがゼロになります。これはZnO-C₃N単層への強いp型ドーピングと見なせます[43,60]。吸着構造の半導体的性質はゼロバンドギャップによる金属的特性で覆い隠される可能性があります[61]が、ゼロバンドギャップは伝導度の可視的な向上をもたらし、ZnO-C₃N単層ベースのデバイス応答性能向上に極めて重要です。応答性能（R）をさらに分析するため、以下の式で計算しました[62]：

$$ R = \frac{\left| \frac{1}{\sigma_{\mathrm{gas}}} - \frac{1}{\sigma_{\mathrm{gas}}} \right|}{\frac{1}{\sigma_{\mathrm{pure}}}} = \left| \frac{\sigma_{\mathrm{pure}} - \sigma_{\mathrm{gas}}}{\sigma_{\mathrm{gas}}} \right| \tag{4} $$

ここでσ_pureとσ_gasは、吸着前と吸着後のZnO-C₃N系の伝導度を表します。計算によれば、H₂S系のRは0.567、他の3系は0.962であり、すなわちZnO-C₃N単層表面で吸着が起こると伝導度がそれぞれ56.7%および96.2%増加することになり、これらのガス存在を検知可能です。

回復時間（τ）は、ガス検知用センサの特性評価におけるもう一つの重要パラメータで、吸着ガス分子を除去するのに要する時間を示します。文献をレビューすると、τはファントホフ-アレニウスの式で計算可能です[63]：

$$ \tau = F^{-1} e^{-E_a / k_B T} \tag{5} $$

ここでFは試行頻度で、本研究では10¹² s⁻¹と定義しました。E_aは脱着エネルギー障壁で、ここではE_adsの値と同一と仮定し、k_BとTは式(3)と同様に定義します。式(5)から、E_adsが大きくなるにつれて吸着ガスの脱着が困難になりますが、作動温度を上げることで制御可能です。

表2に、4種のガスがZnO-C₃N単層表面から除去されるのに必要な回復時間を示します。H₂S分子は小さなE_adsを持つため、脱着エネルギー障壁が低く、マイクロ秒オーダーの極めて短い回復時間となります。一方、大きなE_adsを持つ系では、作動温度ではガス分子を表面から分離することは不可能で、脱着には数日を要します。498 Kおよび598 K（電気機器では稀にしか発生しない温度）に達すると、SO₂とSOF₂それぞれについて脱着プロセスが分オーダーまで加速されます。ガス分子（SO₂、SOF₂、SO₂F₂）と表面間の極めて強い吸着は、ZnO-C₃N単層がSF₆分解生成物を除去し、電力システム内の良好な絶縁状態を維持するガススカベンジャーとしての潜在的応用を示唆します。また、実際の構造では、C₃N上に大量のZnOナノ結晶が存在するため、効果は大幅に向上すると期待されます。さらに、元の構造と比較して、ZnO-C₃N単層から放出されるガスの活性は、分子構造の激しい変形（SOF₂およびSO₂F₂）により大きく損なわれ、システムへの影響はほとんど及ぼせません。H₂Sについては、ZnO-C₃N単層の不安定な相互作用と極めて短い回復時間は、吸着密度が小さいため効果的な検知を提供できないと推測されます。

結論

本論文では、ZnO修飾C₃Nモデルを確立し、ZnO粒子をC₃N表面に様々な向きと位置で配置することで最適構造を探索しました。これにより、ZnO-C₃N単層の4種のSF₆分解生成物（H₂S、SO₂、SOF₂、SO₂F₂）に対する吸着パラメータを、吸着前後のE_ads、DOS、Q_T、バンド構造の分析により取得しました。H₂S分子はナノ構造上で安定吸着しにくい一方、他のガスはZnO粒子に強く捕捉されることが判明しました。これらの結果は、ZnO-C₃N単層の吸着性能が、高電圧装置からSO₂、SOF₂、SO₂F₂を除去するガススカベンジャーとしての潜在的応用を可能にし、絶縁強度と電力システムの安全運転を維持することを確認しました。さらに、フロンティア分子軌道理論は、上記ガス吸着による伝導度の明確な変化を考慮すると、ZnO-C₃N単層がSF₆絶縁装置の誘電状態を評価する指標としての可能性を持つことを示唆しています。

データと材料の入手可能性

全てのデータと材料は本文および補足ファイルに提供されています。