CO与H<sub>2</sub>在应变Fe(110)表面的竞争吸附

李守英; 赵卫民; 乔建华; 王勇

doi:10.7498/aps.68.20190660

摘要

Abstract

作者及机构信息

Authors and contacts

文章全文

1. 引　言
2. 计算模型和方法
3. 计算结果与讨论
3.1 CO分子的吸附
3.2 H2分子的吸附
3.3 势能随距离的变化
3.4 应变对CO与H2吸附的影响
4. 结　论

参考文献

施引文献

摘要

为研究CO降低临氢管线钢的氢脆机制, 采用基于自旋极化密度泛函理论的第一性原理方法, 研究H₂和CO在Fe(110)表面吸附过程中的势能变化以及不同应变时的吸附. 研究发现Fe(110)表面对CO的吸引力大于H₂, 且预先吸附的CO能阻碍H₂的解离, 减弱H与Fe之间的作用力. 态密度分析结果表明CO中的C原子与Fe原子有多个共轭峰, 有强烈的共轭杂化作用. 不同应变Fe(110)表面的吸附结果表明CO在Fe(110)表面的吸附能比H₂更负, CO与表面的结合强度更大, CO优先吸附. 结合热力学定量计算分析CO分压对氢覆盖度影响, 结果表明随着CO的分压升高, 氢覆盖度降低. 表面拉应变越大, 需要的CO分压越高. 拉应变使得H₂, CO吸附能差减小, CO阻碍氢吸附的能力降低是拉应变表面需要更高CO分压的原因.

关键词:

氢脆 /
竞争吸附 /
氢吸附

Abstract

In this work, the competitive adsorption behavior of H₂ and CO on strained Fe(110) are investigated by the first-principles method based on the spin-polarized density functional theory to study the hydrogen embrittlement of steels. The results show that the most stable adsorption site for CO is top site, and the orbital of CO molecule hybridizing with Fe 3p and 4s states illustrates a strong electronic interaction between them. The adsorption energy values of CO at the four calculated adsorption sites are more negative than those of H₂, which favors the binding with Fe(110) surface. The potential energy variations for CO and H₂ molecules close to the surface are calculated. The attractive force of the Fe(110) surface acting on CO in 1.5–3 Å is greater than that acting on H₂. The pre-adsorbed CO increases the dissociation energy barrier of H₂ from 0.08 eV to 0.13 eV but reduces the force between H₂ and surface. The surface tensile strain enhances the interaction between hydrogen and Fe(110), which, however, is reduced by the compressive strain. The opposite tendency is found in the adsorption of CO. The binding strength of CO is stronger than that of H₂ on the strained Fe(110) surface. The difference in adsorption energy between CO and H₂ decreases with tensile strain increasing. The effect of surface strain and partial pressure of CO gas phase on the surface coverage ratio of H atom are also calculated quantitatively based on thermodynamics at 298 K, with the partial pressure of H₂ set to be 10 MPa. The surface ratio of the H atom decreases with partial pressure of CO increasing. The hydrogen coverage drops nearly to zero when the partial pressure of CO reaches a certain value. This result reveals that CO can inhibit hydrogen adsorption on Fe surface. In the case where the surface ratio of hydrogen decreases to 1%, the corresponding CO partial pressures are 105 Pa, 1.1 × 10³ Pa, 2.4 × 10⁵ Pa on –2%, 0, 2% strained Fe(110) surface, respectively. High CO partial pressure is needed to suppress the hydrogen adsorption since the binding strength of CO is close to that of H₂ on the expanded surface.

Keywords:

PACS:

74.70.Xa	(Pnictides and chalcogenides)
74.25.F-	(Transport properties)
75.50.Ee	(Antiferromagnetics)
76.50.+g	(Ferromagnetic, antiferromagnetic, and ferrimagnetic resonances; spin-wave resonance)

作者及机构信息

1.
中国石油大学(华东)材料科学与工程学院, 青岛　266580

2.
青岛滨海学院机电工程学院, 青岛　266555

通信作者: 赵卫民, zhaowm@upc.edu.cn

基金项目: 山东省自然科学基金(批准号: ZR2017MEE005)资助的课题

Authors and contacts

1.
Schools of Materials Science and Engineering, China University of Petroleum (East China), Qingdao 266580, China

2.
Schools of Mechanical and Electronic Engineering, Qingdao Binhai University, Qingdao 266555, China

Corresponding author: Zhao Wei-Min, zhaowm@upc.edu.cn

Funds: Project supported by the Natural Science Foundation of Shandong Province, China (Grant No. ZR2017MEE005)

文章全文

1. 引　言

近年来, 环境保护、地球温室效应、能源枯竭已经成为世界各国共同努力研究的重要课题, 各国都在致力于开发清洁能源和可再生能源^[1]. 太阳能、风能具有丰富、清洁、可再生的优点而受到广泛关注. 但是这些可再生能源具有间歇性、地域性、且不易储存和运输的特点^[2,3]. 相对而言, 氢具有高效、可存储和运输的优点, 被视为最理想的能源载体. 世界各国致力于研究将太阳能、风能产生的电力转换为氢能^[4,5]. 大量、长距离氢能输送, 管道是最高效、低成本的方式^[6,7]. 管线钢在输送氢气的过程中管道内的高压氢环境会引起氢的吸附、渗透, 使得氢进入管线钢内部, 进而导致管线钢的氢致脆化^[8,9]. 氢脆是影响输氢管线工作最突出的因素. 焊接接头区域由于表面存在焊接残余应力及工作应力是输氢管线的薄弱环节^[10]. 表面应力改变材料的表面电子结构, 影响吸附^[11—13]. 氢的解离吸附是氢进入管线钢内部的前提. 目前铁表面的氢吸附主要围绕无应力表面的吸附^[14—16], 关于应力加载铁表面的氢吸附未见报道.

当前, CO和H₂在金属Ni, Fe等表面混合吸附研究目的主要是催化合成碳氢化合物^[17,18], 降低钢的氢脆的研究较少. Bernasek等^[19]通过TEAS (thermal energy atom scattering)试验研究高真空度下CO和H₂在Fe(111)面的吸附, 发现预先饱和吸附CO后, 表面不再吸附H₂. Wang等^[20]研究发现在Fe(111)面, 每个CO分子能阻碍2个H原子的吸附. Huo和Liao^[18]研究发现CO与H₂的稳定吸附构型与表面CO/H₂覆盖度比值有关. 通过比较吸附能, 多数学者们认为CO分子优先于H₂在铁表面吸附^[21,22]. 但是关于CO对H₂的解离过程以及CO分压对吸附氢覆盖度的影响鲜有报道.

Fe(110)面是α-Fe的密排面, 具有最低的表面能, 是塑性变形过程中的新增台阶面^[23,24]. 本文基于密度泛函理论, 研究在应变Fe(110)表面, CO对氢解离、吸附以及覆盖度的影响. 本文研究结果对于CO定量掺杂减少氢吸附、降低临氢管线氢脆具有重要的理论指导意义.

2. 计算模型和方法

本文计算应用Material Studio软件中基于密度泛函理论的CASTEP^[25]模块完成, 电子波函数采用平面波展开, 布里渊区采用Monkhorst-Pack方法均匀k点取样, 应用广义梯度近似的rPBE^[26]泛函(修正后的Perdew-Burke-Ernzerhof泛函)描述电子之间的交换关联能, 利用超软赝势^[27]描述价电子与离子实之间的关系, 选取平面波截断能为425 eV, 自洽场循环收敛于1.0 × 10^–5, 力收敛于0.03 eV/atom, 几何优化采用Broyden-Fletcher-Goldfarb-Shanno算法. α-Fe单胞优化k点取19 × 19 × 19, 优化后的晶胞参数为a = 2.845 Å, 与实验结果2.866 Å基本一致. 基于优化后的单胞模型, 构建7层2 × 2的Fe(110)层模型, 设置真空层厚度为12 Å. 表面模型以及吸附模型几何结构优化时, k点取值为10 × 10 × 1, 底下4层Fe原子固定, 上面3层Fe原子以及吸附原子进行弛豫. 文中所选参数均进行了收敛性测试. 先后将CO与H₂分子放入12 Å × 12 Å × 12 Å的立方盒子进行优化, k点取1 × 1 × 1, 优化后CO分子的键长为1.154 Å, H₂键长为0.748 Å, 分别与实验值1.128 Å, 0.741 Å符合得较好.

Fe(110)表面对称点, 如图1所示. top表示顶位, lb表示长桥位, sb表示短桥位, tf表示三重洞位.

$图 1 Fe(110)面及其对称点\r\nFig. 1. Fe(110) surface and high symmetry sites.$

图 1 Fe(110)面及其对称点

Fig. 1. Fe(110) surface and high symmetry sites.

下载: 全尺寸图片幻灯片

在Fe(110)面的 $[1\bar 11]$ , $\left[ {1\bar 1\bar 1} \right]$ 方向, 以施加等轴应变的方式进行应力加载. 应变ε应用(1)式计算:

$\varepsilon = \Delta l/{l_0},$

(1)

式中l₀表示无应变的晶胞模型尺寸, Δl表示该尺寸的变化量. 应变施加范围为–2%—+2%.

自由气体分子i在表面吸附过程中吉布斯自由能的变化 $\Delta G_{{\rm{ads}}}^i$ 通过(2)式求解:

$\begin{split} & \Delta G_{{\rm{ads}}}^i \\ =\; & E_{{\rm{ads}}}^i - \left( {\Delta H_{{\rm{gas}}}^ \circ \left( {{P^ \circ },T} \right) - T{S^ \circ }\left( {{P^ \circ },T} \right) + {k_{\rm{B}}}T{\rm{ln}}\frac{{{p_{\rm{i}}}}}{{{p^ \circ }}}} \right),\end{split}$

(2)

式中 $E_{\rm ads}^{i}$ 为通过密度泛函理论计算出气体分子i的吸附能, P °为标准大气压, T表示温度, k_B为玻尔兹曼常数, ${\Delta H_{{\rm{gas}}}^ \circ}$ (P °,T)表示从(P °,0 K)到(P °,T)状态的焓变, S °(P °,T)表示状态(P °,T)的熵. ${\Delta H_{{\rm{gas}}}^ \circ }$ (P °,T)与S °(P °,T)的数据可通过查热力学表^[28]计算获得.

H₂与CO的混合吸附过程中, 氢的覆盖度θ_H通过(3)式^[29]计算:

${\theta _{\rm{H}}} = \frac{{{\rm{exp}}\left( { - \dfrac{{\Delta G_{{\rm{ads}}}^{{{\rm{H}}_2}}}}{{{k_{\rm{B}}}T}}} \right)}}{{1 + {\rm{exp}}\left( { - \dfrac{{\Delta G_{{\rm{ads}}}^{{{\rm{H}}_2}}}}{{{k_{\rm{B}}}T}}} \right) + {\rm{exp}}\left( { - \dfrac{{\Delta G_{{\rm{ads}}}^{{\rm{CO}}}}}{{{k_{\rm{B}}}T}}} \right)}}.$

(3)

式中, ${{\Delta G_{{\rm{ads}}}^{{{\rm{H}}_2}}}}$ 表示H₂吸附的吉布斯自由能变化, ${{\Delta G_{{\rm{ads}}}^{{\rm{CO}}}}}$ 表示CO吸附的吉布斯自由能变化.

3. 计算结果与讨论

3.1 CO分子的吸附

CO分子的C原子容易与金属原子作用^[30]. 在CO吸附计算中, C原子靠近表面的对称点top, sb, lb, tf. 几何结构优化完成后, 发现CO垂直吸附于上述对称点. 用(4)式定义CO的吸附能:

$E_{{\rm{ads}}}^{{\rm{CO}}} = {E_{{\rm{slab}} + {\rm{CO}}}} - {E_{{\rm{slab}}}} - {E_{{\rm{CO}}}},$

(4)

式中, E_slab, E_slab+CO分别表示Fe(110)吸附CO分子前后表面的能量; E_CO表示CO分子能量; $E_{{\rm{ads}}}^{{\rm{CO}}}$ 表示CO分子的吸附能.

由吸附能的定义式可知, 负值表示吸附过程放出热量, 正值表示吸收热量. 吸附能的值越负, 表示该吸附构型越稳定. CO以及H₂在不同位置的吸附能如表1所列. 从表1数据可以看出本文计算结果与文献[31]数据符合得很好, 误差仅在0.5%, 说明本文计算结果是可信的. 微小的数据差别可能来自于rPBE与PBE近似的不同. 由表1数据可以看出, 4个位置CO的吸附能都为负值, 且top位的吸附能值最负, 表明该位置是CO分子最优最稳定吸附位置. 低能电子衍射试验也表明CO分子垂直吸附在top位^[32].

表 1 CO与H₂在Fe(110)表面不同位置的吸附能(eV)

Table 1. Adsorption energies of CO and H₂ on high symmetry sites of Fe(110).

吸附物	top	sb	lb	tf
CO	–1.89	–1.64	–1.82	–1.83
CO^[31] (PBE)	–1.88	–1.63	–1.80	—
H₂	–0.47	–0.82	–1.04	–1.33

下载: 导出CSV

| 显示表格

图2为CO分子吸附在top位前后的分波态密度(PDOS)图. 可以看出吸附前, 孤立的CO分子中C原子和O原子有强烈键合作用, 从低能到高能的分子轨道分别是3σ, 4σ, 1π, 5σ, 2π. 其中5σ是最高占据分子轨道, 2π是最低未占据分子轨道. 吸附前, Fe原子的态密度分布在–8—22 eV能量区间. CO吸附使得Fe原子自旋向上态密度和自旋向下态密度差异减小. 吸附后, Fe原子的4s, 3p在–22.6 eV, –8.9 eV处, 3d在–5.7 eV出现较大的共轭峰. 表明表面Fe原子与C原子存在强烈的杂化耦合作用. 该杂化耦合作用使得C原子的2p态展宽, 形成成键态和反键态. 成键态使得C原子与Fe原子的作用增强, 是CO在Fe表面吸附能较大的原因.

$图 2 CO, Fe(110)面及其top位吸附前后分波态密度图　(a)吸附前; (b)吸附后\r\nFig. 2. Projected local density of states from Fe on clean Fe surface and a CO molecular in vacuum (a) and with a CO adsorbed surface (b).$

图 2 CO, Fe(110)面及其top位吸附前后分波态密度图　(a)吸附前; (b)吸附后

Fig. 2. Projected local density of states from Fe on clean Fe surface and a CO molecular in vacuum (a) and with a CO adsorbed surface (b).

下载: 全尺寸图片幻灯片

3.2 H₂分子的吸附

计算过程中发现, 在靠近Fe(110)表面时, 氢分子容易解离成为氢原子吸附在Fe(110)表面. 为提高计算效率, 通过氢原子的吸附根据(5)式计算氢分子的吸附能 ${{ E_{{\rm{ads}}}^{{{\rm{H}}_2}}}}$ ^[29]:

$E_{{\rm{ads}}}^{{{\rm{H}}_2}} = 2{E_{{\rm{slab}} + {\rm{H}}}} - {E_{{{\rm{H}}_2}}} - 2{E_{{\rm{slab}}}},$

(5)

式中E_slab, E_slab+H分别是表面吸附氢原子前后体系的能量; $E_{\rm H_2}$ 表示氢分子的能量. 表1列出了氢分子的吸附能.

氢分子在tf位置的吸附能最负, 为–1.33 eV, 表明tf位的吸附最稳定, 与Wang等^[20]的计算数据–1.36 eV基本一致. 比较CO和H₂的吸附能数据可以看出, 4个吸附位CO的吸附能均比H₂的吸附能更负, 表明CO与Fe(110)表面的结合强度大于H₂与Fe(110)面的结合强度. 其中CO吸附在sb位与H₂吸附在tf之间结合强度差值最小, 为0.31 eV.

计算了H原子在表面tf位吸附前后的态密度, 如图3所示. 从图3可以看出, 吸附后H的1s电子向左移动到–7.3— –5.2 eV, 存在明显的波峰. Fe原子的4s电子在–7.8— –5.8 eV之间出现新的峰、3p电子与3d电子在–6.2 eV出现弱小的峰, 表明Fe的4s, 3p, 3d电子与H的1s电子存在共轭杂化作用, 共轭峰值较小, 说明杂化作用较弱. H原子的吸附对Fe原子自旋向上和自旋向下的态密度差异影响较小.

$图 3 H与Fe(110)面tf位吸附前后分波态密度图　(a)吸附前; (b)吸附后\r\nFig. 3. Projected local density of states from Fe on clean Fe surface and a H atom in vacuum (a) and with a H adsorbed surface (b).$

图 3 H与Fe(110)面tf位吸附前后分波态密度图　(a)吸附前; (b)吸附后

Fig. 3. Projected local density of states from Fe on clean Fe surface and a H atom in vacuum (a) and with a H adsorbed surface (b).

下载: 全尺寸图片幻灯片

3.3 势能随距离的变化

为分析CO对H₂的吸附以及解离过程的影响, 计算H₂, CO分子逐渐靠近Fe(110)表面的势能变化. 计算过程中限定吸附物的Z向坐标, 仅对X, Y坐标进行弛豫.

Wang 等^[33]研究发现垂直构型靠近铁表面的H₂难以解离. 且前面分析表明CO垂直吸附于Fe(110)表面. 因此本文分别计算与表面垂直的CO、平行H₂以及预先吸附CO平行的H₂在靠近表面过程中的势能变化, 计算结果如图4所示. 图中横坐标各距离分别为CO中的C原子到表面的距离、H₂两个H原子的中心到表面的距离以及预先吸附CO时H₂的中心到表面的距离. H₂在靠近表面的过程中, 能量先降低, 后又逐渐升高, 到2.4 Å能量到达最高点. 结构优化中发现距离小于2.4 Å时, 如不固定H₂的Z坐标, H₂直接解离为H原子并吸附在tf位. 距离减小过程中, CO的吸附能不断降低, 到达1.5 Å时, 能量最低, 形成CO的分子吸附态. 能量随距离变化的斜率表示作用力的大小. 1.5—3 Å范围内, 表面对CO的吸引力明显大于H₂, 这是表面优先吸附CO的原因之一.

$图 4 CO与H2吸附过程中的势能变化\r\nFig. 4. Potential energy variations of CO and H2 moving towards Fe(110).$

图 4 CO与H₂吸附过程中的势能变化

Fig. 4. Potential energy variations of CO and H₂ moving towards Fe(110).

下载: 全尺寸图片幻灯片

Gholizadeh和Yu^[34]研究认为小分子的解离过程需要考虑零点能校正. 为准确计算H₂的解离能, 同时提高计算效率, 通过频率计算对吸附H₂, CO+H₂的初始结构以及其在2.4 Å的结构进行零点能校正. 零点能校正公式如(6)式:

${E_{{\rm{ZPE}}}} = \frac{1}{2}\mathop \sum \limits_{j = 1}^{3n} h{\nu _j},$

(6)

式中, ν_j表示振动频率, h表示普朗克常数.

经零点能校正后, H₂的解离能为0.08 eV, 与文献[35]的0.1 eV基本一致, 表明H₂分子易解离. 预先吸附CO后, H₂的解离能升高到0.13 eV, 说明预先吸附CO阻碍H₂的解离, 并且使得最终吸附在tf位的H₂的吸附能由–1.33 eV升高–1.10 eV, 表明预先吸附的CO使得H和Fe(110)面的作用减弱.

3.4 应变对CO与H₂吸附的影响

计算了应变–2%—2%之间Fe(110)面的CO吸附. 发现应变表面sb, lb, tf位置的CO均向top位弛豫, 最终吸附在Fe(110)的top位置. top位CO的吸附能以及tf位H₂的吸附能和应变之间的关系如图5所示, 两者之间的吸附能差如图6所示. 可以看出H₂的吸附能与应变之间存在近似负的线性关系. 拉应变使得H₂的吸附能更负, 压应变使得吸附能值升高. 这表明拉应力使得氢与Fe(110)面的作用力增强, 压应力减小了氢与Fe(110)面的相互作用. 与H₂的吸附相反, 拉应变使得CO吸附能升高, 吸附作用减弱; 压应变使得吸附能更低, 吸附作用增强. 在计算应变范围以内, CO的吸附能均比H₂的吸附能更负, 说明CO与Fe的作用力更强. 通过图6可以看出, 拉应变使得H₂与CO的吸附能差ΔE减小, 压应变使得ΔE增大. 压应变使得CO阻碍氢吸附的效果更为明显.

$图 5 CO吸附能、H2吸附能和应变之间的关系\r\nFig. 5. Relationship between adsorption energy of CO and H2 and strain.$

图 5 CO吸附能、H₂吸附能和应变之间的关系

Fig. 5. Relationship between adsorption energy of CO and H₂ and strain.

下载: 全尺寸图片幻灯片

$图 6 H2和CO吸附能的差与应变之间的关系\r\nFig. 6. Relationship between strain and the difference of CO and H2 adsorption energy.$

图 6 H₂和CO吸附能的差与应变之间的关系

Fig. 6. Relationship between strain and the difference of CO and H₂ adsorption energy.

下载: 全尺寸图片幻灯片

为定量分析CO对应变Fe(110)面氢覆盖度θ_H的影响, 将CO, H₂的吸附能代入(2)式, 利用(2)式计算出的吉布斯自由能的变化 $\Delta G_{{\rm{ads}}}^i$ , 代入(3)式计算298 K时的θ_H. 图7表示298 K, 10 MPa的H₂分压条件下, CO的分压与θ_H之间的关系. 从图7可以看出, 无论表面有无应力, 随着CO的分压升高, 氢覆盖度降低. CO分压升高到一定值时, 氢覆盖度接近0, 表明CO掺杂能抑制氢吸附. –2%, 0, 2%平面应变Fe(110)表面, θ_H降低到1%时, CO分压分别为105 Pa, 1.1 × 10³ Pa, 2.4 × 10⁵ Pa. 拉应变使得H₂和CO吸附能差值减小, CO阻碍氢吸附的能力降低, 需要更高的CO分压降低氢的覆盖度.

$图 7 CO的分压与θH之间的关系\r\nFig. 7. Relationship between CO pressure and coverage of H$

图 7 CO的分压与θ_H之间的关系

Fig. 7. Relationship between CO pressure and coverage of H

下载: 全尺寸图片幻灯片

4. 结　论

本文采用密度泛函理论研究H₂, CO在应变Fe(110)表面的竞争吸附特性. 主要得出如下结论:

1)预先吸附CO阻碍H₂的解离, 使得H₂与Fe(110)面的作用减弱;

2)态密度分析结果表明CO中的C原子与Fe原子存在强烈的杂化耦合作用, 而H原子与Fe原子的杂化作用较弱;

3)无论表面有无应变, CO与表面的结合强度均高于H₂; 拉应变使得两者吸附能差减小, CO阻碍氢吸附的能力降低;

4) CO掺杂对氢覆盖度的计算发现, 随着CO的分压升高, 氢覆盖度降低; 拉应变越大, CO阻碍氢吸附的作用越弱, 需要的CO分压越高; 2%平面应变Fe(110)表面, θ_H降低到1%时, CO分压为2.4 × 10⁵ Pa.

感谢中国石油大学张宏玉老师在计算过程中给予的指导和帮助.

参考文献

[1]	Liu Z, Han H, Xiang C 2018 Energ. Policy 115 92 Google Scholar
[2]	Gorji T B, Ranjbar A A, Mirzababaei S N 2015 Sol. Energy 119 332 Google Scholar
[3]	Mahian O, Kianifar A, Kalogirou S A, Pop I, Wongwises S 2013 Int. J. Heat Mass Transfer 57 582 Google Scholar
[4]	Sherif S A, Barbir F, Veziroglu T N 2005 Sol. Energy 78 647 Google Scholar
[5]	舟丹 2017 中外能源 22 16 Zhou D 2017 Sin. Glo Energ. 22 16
[6]	Dodds P E, Mcdowall W 2013 Energ. Policy 60 305 Google Scholar
[7]	Briottet L, Batisse R, Dinechin A 2012 Int. J. Hydrogen Energy 37 9423 Google Scholar
[8]	Nanninga N, Slifka A, Levy Y 2010 J. Res. Nat. Inst. Stand. Technol. 115 437 Google Scholar
[9]	Kim C M, Kim Y P, Kim W S 2017 J. Mech. Sci. Technol. 31 3691 Google Scholar
[10]	Zhao W, Min Y, Zhang T, Deng Q, Jiang W, Jiang W 2018 Corros. Sci. 133 251 Google Scholar
[11]	Goikoetxea I, Juaristi J, Muino R D 2014 Phys. Rev. Lett. 113 066103 Google Scholar
[12]	Jeon J, Yu B D, Hyun S 2016 J. Kor. Phy. Soc. 69 1776 Google Scholar
[13]	Yang L, Shu D J, Li S C 2016 Phys. Chem. Chem. Phys. 18 14833
[14]	Huo C, Li Y W, Wang J G 2005 J. Phys. Chem. B 109 14160 Google Scholar
[15]	Dan C S 2005 Catal. Today 105 44 Google Scholar
[16]	Kunisada Y, Sakaguchi N 2015 J. Jap. Ins. Met. A 79 447 Google Scholar
[17]	Amaya R S, Linares D H, Duarte H A 2016 J. Phys. Chem. C 120 10830
[18]	Huo C F, Liao X Y 2007 J. Phys. Chem. C 111 4305 Google Scholar
[19]	Bernasek S L, Zappone M, Jiang P 1992 Surf. Sci. 272 53 Google Scholar
[20]	Wang T, Tian X X, Yang Y, Li Y W, Wang J G, Beller M, Jiao H J 2016 Catal. Today 261 82 Google Scholar
[21]	Huo C F, Ren J, Li Y W 2007 J. Catal. 249 174 Google Scholar
[22]	Burke M L, Madix R J 1990 Surf. Sci. 237 20 Google Scholar
[23]	Xie W, Peng L, Peng D 2014 Appl. Surf. Sci. 296 47 Google Scholar
[24]	潘金生, 田民波 2011 材料科学基础 (北京: 清华大学出版社出版) 第156页 Pan J S, Tian M B 2011 Fundamentals of Material Science (Beijing: Tsinghua University Press) p156 (in Chinese)
[25]	Clark S J, Segall M D, Pickard C J,Hasnip P J, Probert M I J, Refson K, Payne M C 2005 Z. Kristallogr. 220 567 Google Scholar
[26]	张凤春, 李春福, 文平 2014 物理学报 63 197101 Google Scholar Zhang F C, Li C F, Wen P 2014 Acta Phys. Sin. 63 197101 Google Scholar
[27]	王明军, 李春福, 文平, 张凤春, 王垚, 刘恩佐 2016 物理学报 65 037101 Google Scholar Wang M J, Li C F, Wen P, Zhang F C, Wang Y, Liu E Z 2016 Acta Phys. Sin. 65 037101 Google Scholar
[28]	Chase M W 1998 NIST-JANAF Thermochemical Tables (4th Ed.) (New York: The American Institute of Physics for The National Institute of Standards and Technology) p641
[29]	Kuwabara A, Saito Y, Koyama Y 2008 Mater. Trans. 49 2484 Google Scholar
[30]	陈运红 2005 硕士学位论文 (广州: 暨南大学) Chen Y H 2005 M. S. Thesis (Guangzhou: Ji'nan University) (in Chinese)
[31]	Jiang D E, Carter E A 2004 Surf. Sci. 570 167 Google Scholar
[32]	Moon D W, Cameron S, Zaera F, Eberhardt W, Carr R, Bernasek L, Gland J L, Dwyer J 1987 Surf. Sci. 180 L123
[33]	Wang H, Nie X, Guo X, Song C S 2016 J. CO₂ Uti. 15 107 Google Scholar
[34]	Gholizadeh R, Yu Y X 2015 Appl. Surf. Sci. 357 1187 Google Scholar
[35]	Staykov A, Yamabe J, Somerday B P 2014 Int. J. Quantum Chem. 114 626 Google Scholar

施引文献

期刊类型引用(1)

李泽众，洪文山，谢涛，刘畅，罗会仟. 铁砷化物超导体的自旋激发谱. 物理学报. 2025(01): 219-233 .

百度学术

其他类型引用(4)

图 1 Fe(110)面及其对称点

Fig. 1. Fe(110) surface and high symmetry sites.

下载: 全尺寸图片幻灯片

图 2 CO, Fe(110)面及其top位吸附前后分波态密度图　(a)吸附前; (b)吸附后

Fig. 2. Projected local density of states from Fe on clean Fe surface and a CO molecular in vacuum (a) and with a CO adsorbed surface (b).

下载: 全尺寸图片幻灯片

图 3 H与Fe(110)面tf位吸附前后分波态密度图　(a)吸附前; (b)吸附后

Fig. 3. Projected local density of states from Fe on clean Fe surface and a H atom in vacuum (a) and with a H adsorbed surface (b).

下载: 全尺寸图片幻灯片

图 4 CO与H₂吸附过程中的势能变化

Fig. 4. Potential energy variations of CO and H₂ moving towards Fe(110).

下载: 全尺寸图片幻灯片

图 5 CO吸附能、H₂吸附能和应变之间的关系

Fig. 5. Relationship between adsorption energy of CO and H₂ and strain.

下载: 全尺寸图片幻灯片

图 6 H₂和CO吸附能的差与应变之间的关系

Fig. 6. Relationship between strain and the difference of CO and H₂ adsorption energy.

下载: 全尺寸图片幻灯片

图 7 CO的分压与θ_H之间的关系

Fig. 7. Relationship between CO pressure and coverage of H

下载: 全尺寸图片幻灯片

表 1 CO与H₂在Fe(110)表面不同位置的吸附能(eV)

Table 1. Adsorption energies of CO and H₂ on high symmetry sites of Fe(110).

吸附物 top sb lb tf

CO –1.89 –1.64 –1.82 –1.83
CO^[31] (PBE) –1.88 –1.63 –1.80 —
H₂ –0.47 –0.82 –1.04 –1.33

下载: 导出CSV

[1]	Liu Z, Han H, Xiang C 2018 Energ. Policy 115 92 Google Scholar
[2]	Gorji T B, Ranjbar A A, Mirzababaei S N 2015 Sol. Energy 119 332 Google Scholar
[3]	Mahian O, Kianifar A, Kalogirou S A, Pop I, Wongwises S 2013 Int. J. Heat Mass Transfer 57 582 Google Scholar
[4]	Sherif S A, Barbir F, Veziroglu T N 2005 Sol. Energy 78 647 Google Scholar
[5]	舟丹 2017 中外能源 22 16 Zhou D 2017 Sin. Glo Energ. 22 16
[6]	Dodds P E, Mcdowall W 2013 Energ. Policy 60 305 Google Scholar
[7]	Briottet L, Batisse R, Dinechin A 2012 Int. J. Hydrogen Energy 37 9423 Google Scholar
[8]	Nanninga N, Slifka A, Levy Y 2010 J. Res. Nat. Inst. Stand. Technol. 115 437 Google Scholar
[9]	Kim C M, Kim Y P, Kim W S 2017 J. Mech. Sci. Technol. 31 3691 Google Scholar
[10]	Zhao W, Min Y, Zhang T, Deng Q, Jiang W, Jiang W 2018 Corros. Sci. 133 251 Google Scholar
[11]	Goikoetxea I, Juaristi J, Muino R D 2014 Phys. Rev. Lett. 113 066103 Google Scholar
[12]	Jeon J, Yu B D, Hyun S 2016 J. Kor. Phy. Soc. 69 1776 Google Scholar
[13]	Yang L, Shu D J, Li S C 2016 Phys. Chem. Chem. Phys. 18 14833
[14]	Huo C, Li Y W, Wang J G 2005 J. Phys. Chem. B 109 14160 Google Scholar
[15]	Dan C S 2005 Catal. Today 105 44 Google Scholar
[16]	Kunisada Y, Sakaguchi N 2015 J. Jap. Ins. Met. A 79 447 Google Scholar
[17]	Amaya R S, Linares D H, Duarte H A 2016 J. Phys. Chem. C 120 10830
[18]	Huo C F, Liao X Y 2007 J. Phys. Chem. C 111 4305 Google Scholar
[19]	Bernasek S L, Zappone M, Jiang P 1992 Surf. Sci. 272 53 Google Scholar
[20]	Wang T, Tian X X, Yang Y, Li Y W, Wang J G, Beller M, Jiao H J 2016 Catal. Today 261 82 Google Scholar
[21]	Huo C F, Ren J, Li Y W 2007 J. Catal. 249 174 Google Scholar
[22]	Burke M L, Madix R J 1990 Surf. Sci. 237 20 Google Scholar
[23]	Xie W, Peng L, Peng D 2014 Appl. Surf. Sci. 296 47 Google Scholar
[24]	潘金生, 田民波 2011 材料科学基础 (北京: 清华大学出版社出版) 第156页 Pan J S, Tian M B 2011 Fundamentals of Material Science (Beijing: Tsinghua University Press) p156 (in Chinese)
[25]	Clark S J, Segall M D, Pickard C J,Hasnip P J, Probert M I J, Refson K, Payne M C 2005 Z. Kristallogr. 220 567 Google Scholar
[26]	张凤春, 李春福, 文平 2014 物理学报 63 197101 Google Scholar Zhang F C, Li C F, Wen P 2014 Acta Phys. Sin. 63 197101 Google Scholar
[27]	王明军, 李春福, 文平, 张凤春, 王垚, 刘恩佐 2016 物理学报 65 037101 Google Scholar Wang M J, Li C F, Wen P, Zhang F C, Wang Y, Liu E Z 2016 Acta Phys. Sin. 65 037101 Google Scholar
[28]	Chase M W 1998 NIST-JANAF Thermochemical Tables (4th Ed.) (New York: The American Institute of Physics for The National Institute of Standards and Technology) p641
[29]	Kuwabara A, Saito Y, Koyama Y 2008 Mater. Trans. 49 2484 Google Scholar
[30]	陈运红 2005 硕士学位论文 (广州: 暨南大学) Chen Y H 2005 M. S. Thesis (Guangzhou: Ji'nan University) (in Chinese)
[31]	Jiang D E, Carter E A 2004 Surf. Sci. 570 167 Google Scholar
[32]	Moon D W, Cameron S, Zaera F, Eberhardt W, Carr R, Bernasek L, Gland J L, Dwyer J 1987 Surf. Sci. 180 L123
[33]	Wang H, Nie X, Guo X, Song C S 2016 J. CO₂ Uti. 15 107 Google Scholar
[34]	Gholizadeh R, Yu Y X 2015 Appl. Surf. Sci. 357 1187 Google Scholar
[35]	Staykov A, Yamabe J, Somerday B P 2014 Int. J. Quantum Chem. 114 626 Google Scholar

[1]	赵有鹏, 刘晓勇, 刘辉, 房坤, 王佳, 罗先甫, 徐宁, 孙绪鲁, 刘煜, 高宇昊, 吴泽鹏, 李雪峰, 张欣耀. Ti-2.5Al-2Zr-1Fe在慢应变速率下的氢脆行为与机理研究. 物理学报, 2024, 73(21): 216103. doi: 10.7498/aps.73.20240896
[2]	盛喆, 戴显英, 苗东铭, 吴淑静, 赵天龙, 郝跃. 各Li吸附组分下硅烯氢存储性能的第一性原理研究. 物理学报, 2018, 67(10): 107103. doi: 10.7498/aps.67.20172720
[3]	徐紫巍, 石常帅, 赵光辉, 王明渊, 刘桂武, 乔冠军. 电化学析氢反应中单层MoSe2氢吸附机理第一性原理研究. 物理学报, 2018, 67(21): 217102. doi: 10.7498/aps.67.20180882
[4]	姜平国, 汪正兵, 闫永播, 刘文杰. W20O58(010)表面氢吸附机理的第一性原理研究. 物理学报, 2017, 66(24): 246801. doi: 10.7498/aps.66.246801
[5]	姜平国, 汪正兵, 闫永播. 三氧化钨表面氢吸附机理的第一性原理研究. 物理学报, 2017, 66(8): 086801. doi: 10.7498/aps.66.086801
[6]	王小霞, 刘鑫, 张琼, 陈宏善. 吸附氢分子的振动态及熵的计算. 物理学报, 2017, 66(10): 103601. doi: 10.7498/aps.66.103601
[7]	刘峰斌, 陈文彬, 崔岩, 屈敏, 曹雷刚, 杨越. 活性质吸附氢修饰金刚石表面的第一性原理研究. 物理学报, 2016, 65(23): 236802. doi: 10.7498/aps.65.236802
[8]	刘秀英, 李晓凤, 于景新, 李晓东. Pd负载共价有机骨架COF-108上氢溢流机理的密度泛函理论研究. 物理学报, 2016, 65(15): 157302. doi: 10.7498/aps.65.157302
[9]	黄向前, 林陈昉, 尹秀丽, 赵汝光, 王恩哥, 胡宗海. 一维石墨烯超晶格上的氢吸附. 物理学报, 2014, 63(19): 197301. doi: 10.7498/aps.63.197301
[10]	卢其亮, 罗其全, 陈莉莉. C@Al12团簇吸附H的密度泛函理论研究. 物理学报, 2010, 59(1): 234-238. doi: 10.7498/aps.59.234
[11]	刘峰斌, 汪家道, 陈大融, 赵明, 何广平. 不同密度氢吸附金刚石(100)表面的微观结构. 物理学报, 2010, 59(9): 6556-6562. doi: 10.7498/aps.59.6556
[12]	戴伟, 唐永建, 王朝阳, 孙卫国. 自制吸附仪储氢性能测试研究. 物理学报, 2009, 58(10): 7313-7316. doi: 10.7498/aps.58.7313
[13]	宋红州, 张平, 赵宪庚. Be(0001)薄膜中的量子尺寸效应与吸附氢的第一性原理计算. 物理学报, 2007, 56(1): 465-473. doi: 10.7498/aps.56.465
[14]	宋红州, 张平, 赵宪庚. 原子氢在Be(1010)薄膜上吸附的第一性原理计算. 物理学报, 2006, 55(11): 6025-6031. doi: 10.7498/aps.55.6025
[15]	于洋, 徐力方, 顾长志. 氢吸附金刚石（001）表面的第一性原理研究. 物理学报, 2004, 53(8): 2710-2714. doi: 10.7498/aps.53.2710
[16]	胡晓明, 林彰达. 用低能电子衍射研究氢在Si(100)表面吸附引起的相变. 物理学报, 1996, 45(6): 985-989. doi: 10.7498/aps.45.985
[17]	谢剑钧, 张涛, 路文昌. 氢在担载金属表面的吸附研究. 物理学报, 1993, 42(11): 1815-1821. doi: 10.7498/aps.42.1815
[18]	侯晓远；杨曙；董国胜；丁训民；王迅. 用HREELS研究氢在GaAs和InP的（111），（III）表面上的吸附. 物理学报, 1987, 36(8): 1070-1074. doi: 10.7498/aps.36.1070
[19]	胡际璜, 刘国辉, 王迅. Si(111)表面吸附氢会形成双氢化相吗?. 物理学报, 1986, 35(9): 1192-1198. doi: 10.7498/aps.35.1192
[20]	金晓峰, 丰意青, 庄承群, 王迅. 用热脱附谱研究氢在Si(100)清洁表面的吸附. 物理学报, 1984, 33(6): 747-754. doi: 10.7498/aps.33.747

期刊类型引用(1)
1. 李泽众，洪文山，谢涛，刘畅，罗会仟. 铁砷化物超导体的自旋激发谱. 物理学报. 2025(01): 219-233 . 百度学术
其他类型引用(4)

计量

文章访问数: 13167
PDF下载量: 245
被引次数: 5

姓名
邮箱
手机号码
标题
留言内容
验证码

搜索

留言板

CO与H₂在应变Fe(110)表面的竞争吸附