做表面量子化学计算的人应该都经历过这种抓狂的时刻模型建立得很规矩结构优化也收敛得很好但最后算出来的吸附能、功函数或者表面态跟实验数据怎么都对不上。我最早做Pt表面的吸附能计算时就在这个地方卡了好几天。参数调了一圈泛函也换了最后才明白问题根本不在数值设置上而是我搭的那个理想截断表面在真实世界里压根就不存在——真实表面的原子会为了降低能量自发重新排列形成一种全新的表面几何和周期结构。这个现象就是量子化学和表面科学里常说的表面重构surface reconstruction。这篇内容想把这几个问题掰开揉碎讲清楚表面重构到底是怎么发生的、有哪些载入教科书的经典例子、在量子化学计算里怎么建模型才能把重构考虑进去、以及实操中怎么躲开那些容易踩的坑。无论你是刚开始接触表面科学的学生还是已经在跑DFT的研究人员这些内容应该都能用得上。1. 切断晶体的那一刀表面重构的物理本质与驱动力想象一块晶体内部的每一个原子都待在周期性势能的能量最低点上上下左右前后都有邻居配位环境完全满足电子轨道也各得其所。现在你沿着某个晶面把它一刀切开两侧露出来的原子立刻变得不一样了它们的头顶少了大半个配位球有些键被硬生生截断指向真空的方向好端端多出一个没有配对的悬挂键。这种状态不是能量最低的稳定态原子必然会做出调整问题只在调整到什么程度。1.1 弛豫和重构根本不是一回事我在帮学生改报告或者审稿的时候经常看到有人把表面弛豫relaxation和表面重构reconstruction混为一谈。这两个词的差别非常本质直接决定了你要不要换超胞、要不要重新建模型。如果表面只发生弛豫表面原子会相对于理想位置在垂直方向移动比如第一层和第二层的间距收缩或者扩张但表面在水平方向的平移对称性完全没有改变。原本是(1×1)周期优化完之后还是(1×1)周期。打个生活化的比方这就像一队人原地调整了前后间距队伍的基本格子还在只是站得或紧或松。金属的很多密排表面比如Cu(111)、Al(111)主要就是这种弛豫行为第一层间距收缩几个百分点往下的层间距交替收缩或扩张但水平对称性始终不变。如果发生重构事情就严重多了表面原子不仅仅在垂直方向动一动还会在平行方向发生大规模重新排布。原来的(1×1)周期变成(2×1)、(1×2)、(√3×√3)R30°甚至(7×7)这样的大周期。简单来说表面自己在重新组队原本整齐的队形被打散形成了一个全新的表面相。这里要记住一个最根本的判别标准重构的本质不是看原子动了多少而是看表面平移对称性是否发生改变。之所以反复强调这个区分是因为它直接决定计算模型只弛豫的表面用一个小超胞就能描述而重构的表面必须用包含重构周期的大超胞建模成本和方法完全不同。很多算出来的结构和实验对不上的案例追根溯源都是从一开始就把这两个概念弄混了。1.2 原子凭什么要重新排队悬挂键、表面能与表面应力那么问题来了原子明明待得好好的为什么会选择离开平衡位置去重新排队归根到底一句话为了降低体系的表面能。先看最直观的悬挂键问题。以硅这种典型的共价晶体为例体相里的每个硅原子都与周围四个硅原子形成四根共价键配位数是4。沿(111)面切开后最表层每个硅原子只剩三根键连向体相朝真空的方向多出一根悬挂键。悬挂键对应的电子态落在禁带里能量高且不稳定会想方设法找伴配对。最直接的办法就是让表面原子动起来邻近的悬挂键相互靠近形成二聚体表面上方补充一些原子形成吸附原子或者把某些区域整体调成另一种堆垛方式。这些大动作可以重新分布电子、压低表面总能量代价是原子偏离原来的理想晶格位置——这就是重构的原始动力。悬挂键只是问题的一部分真实表面还存在表面应力。表面层原子想要的自然晶格常数往往和体相不完全一样因为上方邻居消失后原子间的平衡距离会相应调整。如果应力足够大表面干脆放弃继续呼应体相的晶格通过重构释放一部分应力。所以你会看到表面重构经常伴随着明显的原子面内位移重构后的表面应力通常显著低于理想截断表面。从热力学角度看所谓重构就是表面原子在给定温度和化学环境下找到的那条更省能量的自组织路径。对共价半导体来说悬挂键的能量惩罚太高所以倾向于剧烈地重新成键对金属来说自由电子可以比较灵活地重新分布重构幅度相对温和但像贵金属(110)这类表面依然会形成非常有特色的重构。明白了这一点你就能理解为什么不能想当然地拿一个理想截断表面当真实表面去算——理想表面在能量上就不是一个稳定相。2. 三种载入教科书的经典重构一次看懂表面世界的多样性上面把重构的物理图像和驱动力讲清楚了但光有概念还不够得落到具体结构上才容易立起直观印象。下面我选三个已经被实验反复验证过的经典案例元素半导体Si(111)、化合物半导体GaAs(110)和过渡金属Pt(110)。这三个体系刚好代表了三种截然不同的表面行为——剧烈的大周期重构、只弛豫不重构、温和的中周期重构。把它们放在一起看你对什么表面可能重构、重构长什么样、为什么会这样的判断会有一个比较完整的参照系。2.1 Si(111)-7×7藏了三十年的原子拼图做表面科学的人大概都听过Si(111)-7×7。这个重构在实验上很早就被低能电子衍射LEED看到了干净的Si(111)表面不是原来的(1×1)周期而是出现了一个每个方向都放大7倍的衍射图案。如果每个方向都扩到7倍一个表面单胞里就要装下49个表面原子位。这么大规模的长程重构整个表面学界花了三十多年才把原子细节彻底搞清楚。直到1985年Takayanagi等人在实验衍射数据的基础上提出了著名的DAS模型也就是dimer-adatom-stacking fault model二聚体、吸附原子、堆垛层错三者的组合。大致描述是重构后的表面单胞分成两块堆垛方式不同的三角形区域边界上形成二聚体墙表面上方多出12个吸附原子四个角上留下所谓的corner hole空位。这套精巧结构最直接的收益是表面悬挂键数量从理想截断时的49个大幅降到19个表面能被实实在在压下去一大截。这个例子对做计算的人有两点提醒。第一重构的尺度可以非常大像Si(111)这种看起来简单的表面最小合理模型也要覆盖7×7的单胞相当于表面一层就要摆几十个原子起步。第二这样复杂的结构绝不是靠简单的结构优化就能从理想截断表面里自己长出来的必须先把DAS模型的原子坐标手工搭好再放进DFT里弛豫。不然你优化一万步大概率还是停在某个和(1×1)差别不大的局域极小上。2.2 GaAs(110)为什么它就是懒得重构讲完一个拼命重构的再讲一个反例GaAs的(110)解理面。GaAs是闪锌矿结构沿着(110)面很容易解理而且解理出来的表面在实验上始终是干净的(1×1)周期没有出现Si那样的巨大重构。这是不是意味着表面原子纹丝没动当然不是——它只是没改变周期但原子位置确实发生了明显偏移表面的Ga原子向体相方向沉下去As原子向真空方向顶出来形成大约30度的表面翘曲。这种只发生在垂直方向的调整就是前面说的弛豫。为什么GaAs不需要重构关键在于电荷补偿。切出(110)面后表面同时出现Ga和As两类原子各自带有一个悬挂键一个容易失电子一个容易得电子。电子会自发地从Ga的悬挂键转移到As的悬挂键重新极化之后形成稳定的表面成键轨道原先那些让体系不安的高能悬挂键基本被中和掉了。电子结构稳住了原子自然不需要大动干戈去重新排布。这个案例值得反复体会不是所有半导体表面都必须重构。判断一个表面需不需要重构不能只看材料类型要看表面能不能通过电荷重新分布和几何弛豫把表面能降到足够低。这个道理同样适用于你在计算里去判断某个陌生表面先想清楚化学机制再决定建模方案而不是盲目地往大超胞里塞原子。2.3 Pt(110)-1×2 missing row金属表面的沟槽美学金属表面也有重构而且对催化计算特别重要。拿Pt(110)来说这个表面在真空条件下会形成所谓的missing row结构每两行表面原子中缺一行形成规则的沟槽表面周期从(1×1)变成(1×2)。直觉上可能觉得缺一行原子好像表面原子更少了这怎么可能更稳定但金属和共价晶体不一样它的结合主要靠自由电子的离域和d电子轨道的重新杂化。去掉一行原子之后剩余原子的配位环境发生改变低配位台阶、拐角原子的电子密度重新分布反而让整个表面的能量降了下来。类似的重构在Au(110)、Ir(110)等贵金属表面也都能观察到这说明它是普遍规律而不是个别表面的巧合。对计算化学家来说这个重构的实战意义很大如果你真的想在Pt(110)上研究CO吸附、析氢反应或者加氢反应缺失行边缘上的低配位原子往往就是催化活性最高的位点。你如果用(1×1)理想表面建模型相当于把整个活性位点直接忽略掉了。常见的做法是至少构建2×1的超胞并把missing row的初始几何手动摆进去再优化。材料与表面重构类型核心驱动机制建模注意Si(111)7×7大周期重构悬挂键数量大幅缩减超胞必须覆盖7×7周期GaAs(110)无重构只弛豫电荷转移补偿悬挂键注意面间距弛豫即可Pt(110)missing row 1×2重构d电子重杂化与低配位稳定化至少用2×1超胞还原沟槽结构3. 量子化学计算里如何正确对待表面重构讲了这么多物理图像接下来进入操作层面。在DFT这类量子化学计算里处理表面重构三件事最重要别跳过文献直接闷头优化、超胞和收敛设置要足够容纳重构、用表面能来判断不同重构的热力学稳定性。3.1 先查文献别让结构优化替你回答一切很多同学拿到一个新表面第一反应是搭一个理想截断的slab直接优化指望结构优化自己发现真实的表面几何。这个想法可以理解但在表面重构问题上风险很大。原因在于DFT结构优化实际上是在求解势能面的一个局域极小点而势能面上往往同时存在很多极小值。理想截断表面本身对应一个局部极小比如Si(111)的7×7重构则对应另一个更深的全局极小。从理想表面出发做共轭梯度或准牛顿优化系统很难自行翻越能垒到达更深的极小最终大概率只是做了一次轻微弛豫就停下来跟真实的重构结构毫无关系。所以我现在的工作习惯是拿到一个新表面第一步永远是查文献包括LEED、STM、AFM这些实验表征结果确认目标表面在常温常压或者目标条件下到底呈现什么表面周期。在Materials Project这类数据库和文献里通常能找到线索。只有在确认了重构模式之后才去手工搭建对应超胞把重构结构的初始坐标放进去再让DFT精细弛豫。这样结构优化只是修正细节而不是猜全局结构成功率会高很多。3.2 slab建模、超胞选择与收敛标准具体到slab模型有几个参数直接影响能否正确描述重构。首先是超胞尺寸。重构的本质是平移对称性改变新周期如果在某个方向是原来的2倍你的超胞在这个方向就必须至少是原来的2倍如果是7×7那就得是7倍。否则无论如何优化这个方向都不可能演化出重构模式因为周期性边界条件不允许出现比超胞更大的重构波矢。这是一个硬约束不随算法优化而改变。其次是slab厚度和真空层。厚度不足表层的弛豫和重构行为会被底层拽住得到不真实的几何。我一般从5层左右开始测试逐步加到8层以上直到表面能和原子位移都收敛。真空层至少15埃起步最好做到20埃甚至更高避免上下两个表面通过周期性镜像互相作用。对半导体和绝缘体如果slab两侧都有悬挂键可以考虑用假氢终止或者干脆做两个面都是重构面的对称slab。再就是收敛标准。重构前后的能量差常常很小对金属表面而言不同重构模式的相对稳定性可能只有每个原子几十毫电子伏特的量级。因此结构优化的受力收敛至少要设到0.02电子伏特/埃甚至更紧能量收敛标准同步提高配合足够密的k点网格和平面波截断能收敛测试。别小看这些设置很多重构算不出来或者两种重构能量一模一样的结论最后查下来都是收敛标准太松导致的假象。3.3 表面能对比怎么判断哪种重构是热力学赢家判断一个重构在热力学上是不是更稳定标准做法是比较表面能surface energy。对单质表面定义很直接γ (E_slab − N × E_bulk) / (2A)。其中E_slab是slab的总能量E_bulk是体相中每个原子的能量A是slab一个面的面积分母中的2来自slab的上、下两个面。实际操作有几个细节要当心。第一如果slab上下两个面不对称常见做法是上面自由弛豫、下面固定模拟体相这种情况下分母里的2就不太严格因为上下两个面的结构和能量贡献不一样。解决的办法是使用对称slab方案或者在非对称slab上加入偶极修正dipole correction否则体系存在净偶极矩总能量和功函数都会受影响。第二比较1×1、2×1、7×7这些不同重构时N的取值必须一致而且slab要做到足够厚确认表面能已经收敛到不随厚度变化。第三如果是化合物表面还要考虑化学势窗口比较不同化学计量比表面时要在元素化学势允许的范围内讨论。有了表面能这把标尺就可以把不同重构模式放在一起比表面能越低的结构在给定热力学条件下越可能出现在真实表面。不过也要记住DFT算的是0 K的静态能量不含有限温度下的振动自由能。如果实验观察到的高温重构在0 K表面能比较里并不是最低的别急着怀疑算错了很可能是温度带来的熵效应在起作用。4. 实用判断框架与避坑指南前面讨论的表面重构基本假设是干净的表面在真空中达到热力学平衡。但在真实应用场景里表面几乎不可能一尘不染催化反应中有反应分子吸附电化学界面有电解质和偏压外延生长中有源源不断的沉积原子。所以这一节把目光放到更实用的层面。4.1 吸附可以诱导重构表面从来不是静止的吸附诱导重构是被很多人忽略但又非常重要的现象。某些重构结构只有在吸附特定分子之后才会稳定存在。比如氧在Cu(110)表面会诱导出一系列基于Cu-O链的表面重构其中最典型的就是(2×1)结构CO在某些金属表面吸附后也能稳定或改变原本的missing row重构半导体表面吸附氢、氧之后原来的大周期重构同样可能被解锁成另一种结构。表面和吸附物是一个耦合体系不是简单的刚性衬底加分子。对做催化计算的人来说这个问题很现实你选了一个(1×1)或者(1×2)表面模型吸附质一放上去最上层原子真的还停在原来的位置上吗我现在的习惯是每做完一轮吸附优化都检查一下固定层以上的原子位移如果吸附后最上层原子整体偏离了原始重构位移达到百分之几埃的量级那大概率说明你需要考虑吸附诱导重构至少要把这种重构纳入过渡态的初始结构设计。表面是活的它会对外界环境做出响应。4.2 我踩过的四个坑希望你别再踩这些年我在表面计算里反复踩过的坑归结起来有四个个个都是血泪教训。第一个是把(1×1)理想表面当成万能起点。就是我开头说的那次经历把Pt(110)当成完美(1×1)面算吸附能怎么都不对。后来想清楚了真实条件下的Pt(110)存在missing row重构(1×1)模型从根上就错了后面所有努力都是在错误的地基上盖楼。第二个是超胞不够大重构模式根本装不下。如果一个表面在文献里是(2×1)重构你却用(1×1)超胞去跑优化程序在数学上就不可能给出(2×1)的原子排布周期性边界条件锁死了可能解的空间。这个坑很隐蔽因为整个优化过程看起来一切正常能量也一路在降你以为自己在找重构其实早就在一开始就被超胞框死了。第三个是优化完只看总能量不看对称性。我见过不止一次跑完结构优化发现能量降低了就以为找到了重构。但把终态原子坐标和初态一对比对称性根本没变只是层间距收缩了一点。这只能叫弛豫不叫重构。判断重构一定要回到那个根本标准表面平移对称性是否改变。第四个是完全不顾温度条件。有些重构是低温或高温稳定相比如某些半导体表面在高温下会从(1×1)变成(2×1)到低温又是另一种结构。DFT算的是0 K的平衡结构实验是在特定温度下测的直接拿0 K结论去比实验对不上很正常。至少要心里有数你算的到底是哪个热力学极限下的稳定性。4.3 我的表面建模检查清单最后把我每次建表面模型都会过一遍的清单放在这里希望能帮你少走弯路查文献或数据库确认目标表面在相关条件下是(1×1)还是某种重构明确重构周期比如2×1、7×7、√3×√3R30°。根据重构周期确定最小超胞超胞尺寸必须能容纳该周期研究吸附时最好在重构单胞基础上再扩展一到两倍降低吸附质之间的镜像相互作用。建slab厚度5层起测真空层15到20埃底层固定一到两层模拟体相顶层自由弛豫。手动按重构结构摆放初始原子坐标然后交给DFT做结构优化。优化收敛后用两个标准验证是否真的得到了重构一是对比初末态原子坐标和表面平移对称性二是计算表面能和(1×1)或其他重构模式做比较。如果是化合物表面明确化学势区间和终止面不要拿非化学计量比的模型直接比较表面能。做吸附或催化计算时记录吸附后最上层原子位移判断是否需要考虑吸附诱导重构。这份清单不是什么金科玉律但它是在无数个白跑通宵的计算里提炼出来的。当你新建的表面模型又一次跟实验对不上时大概率能在清单的某一项里找到原因。说回文章最开始那个让我纠结了好几天的Pt(110)吸附能。后来把模型换成含missing row重构的(1×2)超胞同一套泛函、同一套k点吸附能一下子就对了。那次经历给我的收获是表面重构不是一个可以事后补救的细节而是决定整个模型真实性的前置条件。现在每当要研究一个新表面我第一个问题永远是——这个表面在真实环境里到底长什么样想清楚了再动手比埋头优化几圈参数值钱得多。