γ-Graphyne对CH4、CO、CS2和Cl2吸附行为的第一性原理研究：位点依赖的电子响应与Cl2优先敏感性

《physica status solidi (b)–– basic solid state physics》：CH4, CO, CS2, and Cl2 Adsorption on γ-Graphyne: A First-Principles Study

　　

采用色散校正密度泛函理论（DFT），在四种高对称位点上比较了CH4、CO、CS2和Cl2在原始γ-Graphyne上的吸附行为，分析内容包括吸附能、优化几何构型、Kohn–Sham能带隙、投影态密度（PDOS）及L?wdin布居。CH4以物理吸附为主，吸附能为?0.25至?0.27 eV，几乎不改变0.51 eV的能带隙。CO和CS2的最低能量态（分别为?0.26和?0.33 eV）同样引起可忽略的能带隙变化，而扰动更强的构型在能量上不利。相比之下，Cl2在炔碳位点附近初始化后形成两个1.78 ?的C─Cl键，Cl─Cl键长从2.01 ?伸长至2.88 ?，吸附能为?0.42 eV，能带隙降至0.47 eV。L?wdin布居表明对称的轨道重排，未分辨出净电荷转移，与以共价化学吸附为主的机制一致。在sp²碳附近初始化的构型产生最大的能带隙缩减（27.18%）以及约0.16–0.19个电子从γ-Graphyne转移的估算值，但该构型在能量上不利（吸附能 = +0.11 eV）。上述结果在所考察分子中识别出对Cl2的优先计算响应，同时界定了与传感相关的构型范围。

**研究背景与问题**

二维材料因其大的可及表面积而成为气体传感的重要平台，纳米结构化和化学功能化可进一步调控灵敏度和选择性。石墨烯虽具有高载流子迁移率和优异导电性，但其半金属电子结构和零带隙限制了需要可控电导调制的应用。γ-Graphyne由Baughman等理论提出、Li等实验报道，含33.33%炔基基团，兼具sp²和sp杂化碳原子，具有本征带隙和高载流子迁移率，被视为石墨炔家族中最稳定的成员。CS2是毒性前体，用于纤维素、粘胶纤维和橡胶生产；Cl2是毒性腐蚀性气体，用于水处理、消毒和聚氯乙烯生产；CH4具有可燃性和窒息危险；CO是不完全燃烧的高毒产物。此前研究已考察石墨炔与多种气体的相互作用，但CH4、CO、CS2和Cl2在原始γ-Graphyne上的吸附尚未在统一计算框架下系统比较，尤其对Cl2和CS2的电子响应研究不足。因此，研究人员采用范德华（vdW）校正的密度泛函理论（DFT），在H1、H2、T1和T2四个位点上考察四种分子的吸附，分析吸附能、优化几何、能带隙变化、PDOS和L?wdin布居，以识别优先构型并比较吸附物诱导的电子响应。

**关键技术方法**

研究人员采用Quantum ESPRESSO软件包进行第一性原理DFT计算，使用GGA-PBE泛函加Barone等参数化的vdW校正，辅以LDA-PZ仅作结构基准。采用周期性边界条件、平面波基组和超软赝势，动能截断45 Ry、电荷密度截断450 Ry，结构优化使用4×4×1 Monkhorst–Pack网格、DOS/PDOS计算使用8×8×1网格。吸附模型采用含48个碳原子的2×2超胞，晶格参数13.4764 ?，真空层15–18 ?。分子在四个高对称位点以约3.0和1.5 ?两种初始距离放置，进行无约束全结构弛豫。L?wdin布居分析通过projwfc.x模块将Kohn–Sham态投影到正交化赝原子轨道上获得。

**研究结果**

*结构基准与内聚能*：GGA-PBE+vdW给出的晶格参数a = b = 6.88 ?与实验值6.9 ?最为接近（偏差0.26%），Csp─Csp、Csp─Csp²和Csp²─Csp²键长与文献模拟吻合良好。内聚能为?8.623 eV/atom，与Kenarsari等报道的?8.564 eV/atom接近（差异0.7%）。原始γ-Graphyne的Kohn–Sham带隙为0.51 eV，与先前半局域DFT结果一致。

*CH4吸附*：CH4在H1、T1和T2位点为弱放热吸附（Eads = ?0.25至?0.27 eV），H2位点不利（+0.32 eV）。分子–表面距离3.46–3.80 ?，C─H键保持1.10 ?，带隙基本不变，PDOS显示分子态仅出现在带边以下深处（约?4.2和?11.8 eV），证实为物理吸附且电子响应有限。

*CS2吸附*：CS2存在两种吸附 regime。初始距离约3.0 ?的"_far"构型在四个位点均呈弱放热物理吸附（?0.31至?0.33 eV），H1和H2为最低能量位点，分子内C─S键保持1.56 ?。初始距离约1.5 ?时，H1和H2位点弛豫回物理吸附构型；而T1和T2位点形成S─C共价键，C─C距离降至约1.5 ?，C─S距离约1.85 ?，分子内C─S键变为1.61和1.79 ?，但吸附能强正（+0.54 eV），能量上不可行。T2位点PDOS显示分子态更离域并接近导带边，伴随13.22%的带隙缩减，但该构型能量不利。

*CO吸附*：CO以氧朝下、碳朝下和平行三种取向在四个位点采样。氧朝下取向在各构型中均表现为物理吸附（?0.21至?0.24 eV），C─O键保持1.14 ?，带隙约0.50–0.51 eV。碳朝下取向多数构型仍为物理吸附，但CO_vertical_C_near_T2显示C─C距离1.44 ?、C─O键长1.21 ?，吸附能近零（?0.01 eV），指示部分轨道杂化。平行取向中CO_horizontal_near_H2为最低能量构型（Eads = ?0.26 eV，带隙0.51 eV）；CO_horizontal_near_T2虽形成短接触（C─C 1.44 ?、C─O 1.42 ?），但吸附能强正（+2.88 eV），能量不利。研究人员指出H2位点对CO物理吸附在能量上更有利，而T2位点仅在强制接近条件下促进化学键形成且能量代价显著。

*Cl2吸附*：Cl2表现出最强的位点依赖性。"_far"构型在四个位点均为物理吸附（?0.27至?0.30 eV），Cl─Cl键保持2.00–2.01 ?。初始距离约1.5 ?时，H1和H2位点仍保持物理吸附（?0.30 eV）；T1位点Cl─Cl距离增至3.88 ?，形成两个不等价C─Cl接触（1.93和1.80 ?），产生最大带隙缩减（27.18%），L?wdin分析表明约0.16–0.19个电子从γ-Graphyne转移至Cl2，但Eads = +0.11 eV能量不利。T2位点形成两个近等价C─Cl键（1.779和1.780 ?），Cl─Cl距离增至2.882 ?，Eads = ?0.42 eV为最放热构型，带隙降至0.47 eV，L?wdin布居显示强对称轨道重排但无净电荷转移，指示以共价C─Cl键合为主并伴随强Cl─Cl活化。PDOS分析表明T2位点Cl态与薄片广泛杂化，与短C─Cl键和有利化学吸附一致；T1位点Cl态直接贡献于两个带边，产生最大带隙调制但能量不利。

**讨论与结论**

研究人员指出，吸附物诱导的带隙变化仅作为计算电子响应描述符，不能直接等同于器件灵敏度。对于CS2和CO，化学吸附样构型虽产生约13%和11%的带隙缩减，但正吸附能使其相对分离参考体系不利；最低能量物理吸附构型的带隙变化低于2%。CH4同样弱相互作用且电子扰动小。这些结果表明增强对CH4、CO和CS2的响应可能需要功能化或杂原子掺杂。Cl2在所考察四种气体中产生最显著的电子响应，能量有利的T2态结合共价C─Cl键形成与可测带隙调制，而T1态显示更强电荷转移和前沿态扰动但吸附能为正。计算局限性包括：描述的是自由悬浮原始γ-Graphyne在静态0 K总能量框架下、2×2超胞施加的周期性吸附物覆盖度下的行为；衬底可能引入应变、屏蔽静电相互作用、转移电荷、与γ-Graphyne态杂化或阻塞吸附位点；有限温度将增加振动和构型熵、覆盖度效应及吸附/脱附动力学。因此能量排序和DOS趋势应视为本征0 K参考，需显式衬底、有限温度自由能或从头算分子动力学及反应路径计算来确立操作温度稳定性和恢复性。

研究结论部分翻译如下：色散校正DFT计算揭示了原始γ-Graphyne对所考察气体的吸附物依赖响应，定性顺序为Cl2 >> CS2 ≈ CO >> CH4。CH4物理吸附（Eads = ?0.25至?0.27 eV）并产生可忽略的带隙变化。CO和CS2的最低能量态同样为物理吸附（Eads = ?0.26至?0.33 eV），带隙修饰低于2%。化学吸附样的CO和CS2构型产生更大带隙缩减，但正吸附能使其能量不利。Cl2呈现两种 distinct 化学吸附 regime：T2位点形成两个近等价C─Cl键（1.779和1.780 ?），Cl─Cl距离从2.01增至2.882 ?，该构型能量有利（Eads = ?0.42 eV）并将带隙降至0.47 eV，L?wdin布居显示强对称轨道重排但无分辨净电荷转移，指示以共价C─Cl键合为主和强Cl─Cl活化；T1位点Cl─Cl分离达3.88 ?，形成两个不等价表面接触（1.93和1.80 ?），转移约0.16–0.19个电子从薄片至Cl2并产生最大带隙缩减（27.18%），但Eads = +0.11 eV使其能量不利。综合吸附能、优化几何、PDOS、带隙变化和L?wdin布居，研究人员将γ-Graphyne鉴定为在所考察四种分子中对Cl2具有优先电子响应的候选材料，但这是比较性第一性原理预测而非已证明的器件选择性，需实验输运测量及包含衬底、有限温度自由能、吸附/脱附势垒和真实覆盖度的计算来量化灵敏度和恢复性，并可探索功能化或杂原子掺杂以增强对CH4、CO和CS2的有利态响应。