自从1925年h.s taylor提出催化活性中心的概念以来，人们试图通过各种方法与理论理解催化活性的起源，以期能够快速而准确地预测活性中心的位置，并达到设计高活性催化材料的目的。这其中，基于中间物吸附/脱附以及d带中心等理论的密度泛函理论计算取得了巨大的成功。但是这种方法算量巨大，对算力和人力都提出了极高的要求，此外也很难对材料的所有晶面做出全面的分析，极大限制了其对新材料的预测功能。那么，是否有一种理论，可以绕过复杂的第一性原理计算，实现催化活性中心的快速判断呢？
近年来，拓扑能带理论的发展为高通量预测材料的电子结构提供了新的机遇。以拓扑量子化学理论（topological quantum chemistry）为代表（, 2019, 566, 480），由于晶体能带可由空间群的不可约表示来标记，因此可以从晶体对称性出发，对材料的电子结构和性质做出预测。近期，中国科学院宁波材料技术与工程研究所磁性相变材料团队研究员李国伟在系统研究拓扑表面电子态表面化学效应的基础上（npj computational materials, 2021, 7, 1；applied physics letters, 2020, 116, 070501； advances, 2019, 5, eaaw9867；angewandte chemie, 2019, 131, 13241），与德国马普固体物理化学研究所教授claudia felser及美国普林斯顿大学教授andrei bernevig合作，通过扫描无机晶体数据库中数万个拓扑平庸材料，确定了一种来源于材料体相的本征表面电子态—obstructed surface states（osss），并可以用来快速准确判定无机晶体材料催化活性中心的可能位置。
研究团队从无机晶体数据库中的3万余种非拓扑材料出发（arxiv:2111.02433），计算了这些材料的实空间不变量（real-space invariants），发现存在着这样一类材料，即obstructed wannier电荷中心（owccs）与晶体的原子位置不重合时，会出现一种内禀的金属表面电子态（osss）。分析发现，这些金属电子态完美满足了催化材料活性位点的基本要求：当解离出该晶面时，会有表面态产生；具有良好的导电性；电子态位于费米面附近，易参与电子转移与成键。
为了证实这一猜想，研究人员以多种经典催化制氢材料为例，合成了肉眼可分辨不同晶体面的大尺寸单晶，证实了催化活性位点与osss所在晶面高度重合。以层状范德华材料mox2（x=s, se, te）为例，理论表明，此类材料的owccs位于晶体wyckoff位置2b处，这也预示着表面电子态不在晶体（001）位置，而是位于所有的垂直于（001）晶面的侧面内。这与众多的实验结果相符，即催化活性来源于非（001）晶面的侧晶面。通过导电探针原子力显微镜，他们发现导电性在晶体的边缘位置明显高于无osss存在的惰性面。更重要的是，可以通过考虑晶体对称性及owccs所在位置，判定osss存在的几率并用于预测潜在的高活性催化材料。
借助于这一机理，研究人员发现了3000余种类似的非磁性材料与30种磁性材料。考虑晶体化学稳定性及合成上的可能性及成本，他们从中预测了超过400种潜在的高活性催化材料，并详细给出了osss所在的miller指数面。这证实了osss可以作为寻找催化活性中心的“描述符”，为快速判定活性位点位置、进而针对不同的催化场景设计高效催化材料提供了指导。
这一成果以obstructed surface states as the descriptor for predicting catalytic active sites in inorganic crystalline materials为题发表在advanced materials上，并被遴选为当期卷首插画（frontispiece）论文。
图1 基于拓扑量子化学理论的催化材料搜寻与设计思路
图2 通过本征表面态快速确定催化活性中心所在位置
研究团队单位：宁波材料技术与工程研究所
