VASP计算吸附能

一、项目背景与科学问题

水在金属氧化物表面的吸附行为是光催化、湿敏传感器、表面润湿与电化学界面研究的基础问题。二氧化钛（TiO₂）作为最经典的宽禁带光催化材料，其（110）面是锐钛矿相热力学最稳定的暴露面，理解水分子在 TiO₂(110) 表面的吸附构型、吸附强度与解离倾向，对揭示光催化产氢、光生载流子分离以及表面羟基化过程至关重要。

本项目的核心科学问题在于：水分子在 TiO₂(110) 表面究竟是分子吸附还是解离吸附？两种吸附模式的能量差是多少？不同表面位点（五配位钛位 Ti₅c、桥氧位 O_br）的吸附能力差异如何？这些问题的定量答案需要第一性原理计算给出。

二、计算方法与理论背景

2.1 吸附能定义

吸附能定义为 E_ads = E(slab+H₂O) − E(slab) − E(H₂O)，其中 E(slab+H₂O) 为吸附体系总能量，E(slab) 为干净表面能量，E(H₂O) 为孤立水分子能量。吸附能为负表示放热吸附，数值越大吸附越强。该定义忽略了基组重叠误差与零点能，因此项目中对关键构型进一步做了零点能修正。

2.2 DFT 与 Hubbard U

常规 GGA-PBE 泛函会严重低估 TiO₂ 中 Ti 3d 态的带隙与局域性，导致吸附能偏差。本项目采用 PBE+U 方法（对 Ti 3d 施加 U=4.2 eV），并用 DFT-D3 色散修正捕捉范德华相互作用，使水分子与表面的弱相互作用得到合理描述。

三、计算设置与模型构建

• 软件与方法：VASP，PAW 赝势，PBE+U（U=4.2 eV，Ti 3d）+ DFT-D3。
• 表面模型：锐钛矿 TiO₂(110) 四层原子层板，采用 2×2 表面超胞，真空层厚度 15 Å，底部两层原子固定模拟体相约束。
• 平面波截断能：450 eV；k 点网格：3×2×1（含 Gamma 点）。
• 收敛标准：电子自洽 10⁻⁵ eV，离子弛豫力收敛 0.02 eV/Å。
• 计算内容：分子吸附（水通过 O 端配位 Ti₅c）、解离吸附（OH+H 分别吸附于 Ti₅c 与 O_br）、O_br 位 H 供体吸附三种构型，并进行 Bader 电荷分析。

四、结果分析与案例图解读

图 1 的左面板（a）给出三种吸附构型的吸附能柱状图。分子吸附于 Ti₅c 位最强，E_ads = −0.75 eV，表明水分子更倾向于以完整分子形式通过氧孤对电子配位表面钛原子；O_br 位 H 供体吸附为 −0.51 eV，弱于 Ti₅c 位；解离吸附为 −0.42 eV。这一位点依赖关系说明，表面阳离子位点（Ti₅c）是吸附的活性中心，而桥氧位仅通过氢键贡献较弱的吸附。

图 1 的右面板（b）给出分子吸附与解离吸附之间的能量关系示意。以气相水和干净表面为零点，分子吸附态能量为 −0.75 eV，解离吸附态为 −0.42 eV，两者之差 ΔE = 0.33 eV，说明在零温零压下分子吸附是热力学稳定构型，解离吸附需要跨越 0.33 eV 的能垒。结合 Bader 电荷分析，分子吸附时水分子向表面转移约 0.12 e，而解离吸附后 OH 与表面成键更共价，电荷转移更大。这一结论与实验上低温下 TiO₂(110) 表面水以分子吸附为主、升温后部分解离的现象一致。

两面板结合，左面板定量比较了不同位点的吸附强度，右面板揭示了分子与解离吸附的热力学关系，共同回答了水在 TiO₂(110) 表面吸附模式与吸附能的核心问题。

五、结论与展望

本项目用 VASP 的 PBE+U 方法计算了水在锐钛矿 TiO₂(110) 表面的吸附能，得到 Ti₅c 位分子吸附 −0.75 eV、解离吸附 −0.42 eV 的结果，明确了分子吸附为热力学稳定构型。案例图以位点吸附能柱状图与吸附模式能量图双面板呈现。

后续可拓展：(1) 计算分子吸附到解离吸附的完整过渡态，得到活化能与动力学信息；(2) 研究水覆盖度对吸附模式转变的影响；(3) 考察掺杂（如 N、Fe）对表面亲水性与吸附能的调控机制。