吸附动力学模拟分析表面吸附模型的化学键 _吸附

所属部门：研之成理理论化学研习社

文章插图
量子化学计算从薛定谔波动方程出发。在分子体系的计算中，我们通过波函数分析可以得到一系列信息。理论化学最主要的目的是归纳化学反应的特征、预测化学反应。因此而发展出了各种概念来描述电子结构，比如分子轨道（MOs）、局域分子轨道（LMOs）、键级、AIM（atoms in molecules）、费米空穴（Fermi Hole）相关的概念，以及各种电荷、各种能量分解方法。
而在表面科学中，分子和载体的化学相互作用同样重要，是分析表面催化反应，表面重构等一系列表面科学重要方法。今天就讲讲两种常用化学键分析方法在周期性体系中的应用：（1）SSNBO（Solid State NATURAL BOND ORBITAL ANALYSIS）和（2）SSAdNDP（Solid State Adaptive Natural Density Partitioning）
（1）自然键轨道理论（NBO）是一种对密度矩阵部分对角化，从而将分子轨道部分定域化的量子化学理论。广义地说，根据对角化和定域化的程度不同，该理论中研究的轨道包括自然原子轨道，自然杂化轨道（缩写：NHO），自然键轨道（缩写：NBO）和自然半定域化分子轨道（缩写：NLMO）。这些自然轨道可以视作从原子轨道线性组合得到分子轨道的中间步骤，按照定域化的程度由低到高，它们的关系如下：原子轨道→ NAO → NHO →NBO → NLMO → 分子轨道
在计算化学中，自然（定域化）轨道用于计算电子密度在原子上与在分子间的键上的分布。这些轨道在相应的单原子或双原子区域内具有“最大占据数”的特点。也就是说，以自然轨道为基表示一阶约化密度算符时，矩阵的对角元能够尽可能大，通常可以非常接近或达到2 。于是自然键轨道就给出了波函数对应的最主要的自然路易斯结构。自然 Lewis 结构上对应的自然键轨道占据数通常包括了绝大部分的电子密度，对于常见的有机分子，能达到 99% 以上。自然轨道的概念最早由 Per-Olov L?wdin 在1955年引入，指的是对应于一个多电子波函数的唯一的一组具有最大占据数的单电子波函数。
在分子体系中，NBO已经成为被广泛接收的波函数分析手段。在2012年，J. R. Schmidt发布了可用于处理周期性体系的Periodic NBO software（http:// schmidt.chem. wisc.edu / nbosoftware）下面我们就以CO分子在Pd（111）表面上吸附的例子来说明NBO在催化体系中的具体应用。如图一（A）Pd-C的成键轨道占据1.98符合一个双中心键特征，其中C原子发生sp杂化，在成键过程中，s贡献63%，p贡献37%；而Pd则主要由dz2（78%）参与成键。如图一（D）C-O的pi反键空轨道和Pd的dxz轨道发生部分重叠，则Pd上的部分电子转移到CO的pi*反键上（占据数0.202）形成弱反馈pi键。同样我们可以分析其它几个化学键的成键方式和成分。

文章插图
图一，CO在Pd（111）表面上吸附的NBO轨道。（a）Pd-C sigma轨道，(b) C-Osigma轨道，（c）C-O pi成键轨道，（d）C-O pi反键轨道，Pd上d轨道孤对电子。

文章插图
表一，CO在Pd(111)上的化学键占据数，轨道贡献和轨道成分。
补充一点。在投影到高斯基组上的波函数也可以用于Solid State Adaptive Natural Density Partitioning(SSAdNDP) 分析周期性材料中的多中心键，尤其是对于体系中存在大pi键的二维材料可以给出直观地化学图象，如图二。

文章插图
图二，ssAdNDP 分析Be5C2体系中的都中心化学键
【吸附动力学模拟分析表面吸附模型的化学键】NBO对化学键分析中强大的作用，但是要描述两个化学片段（ex：CO/Pd）之间的电子流向和作用能量分解就需要ETS-NOCV和pEDA方法来实现，请看理论化学研习社下期文章(*^_^*) 。
相关内容阅读：
从量子力学到密度泛函理论（一）
全硼富勒烯的发现 —— 实验与计算的完美结合
J. K. N?rskov：利用理论计算来设计新的固体催化剂