物理吸附和化学吸附在DFT里怎么区分？构型搜索、实空间函数与轨道响应分析

说明：本文华算科技主要介绍物理吸附和化学吸附在 DFT 中的状态定义、候选构型搜索、实空间函数、轨道响应、振动信号与反应坐标分析。

物理吸附和化学吸附的计算定义？

吸附类型描述界面相互作用的性质。物理吸附主要由色散、静电和诱导极化稳定，吸附物原有化学身份与分子骨架大体保留；化学吸附出现表面—吸附物新键、配位变化或分子内键断裂。DFT 输出总能、电子密度与 Kohn–Sham 轨道，physisorption 或 chemisorption 需要研究者依据同一模型的一组信号作出归属。

结构优化找到的是给定初始构型附近的局域极小值。单一初始高度或单一取向可能只落入一个势阱，另一个吸附态仍留在势能面上。相同吸附物应布置多种高度、朝向和位点，分别优化完整分子态、配位态及可能的解离态；自旋、电荷、覆盖度、晶面和超胞在比较中保持一致。

TiO₂(001) 上的 CO₂ 给出直观对照：线性分子远离表面时保留 O=C=O 骨架；C 原子或 O 原子进入表面配位后，构型转为单齿、双齿或桥式碳酸盐。化学身份与配位拓扑发生改变，比任意吸附能阈值更接近两类吸附的定义。

弱吸附构型也可能具有可观的总稳定能。大分子与孔壁接触面积增大后，许多微弱色散项累加，E_ads 的绝对值可超过小分子成键体系。能量大小反映选定参考态之间的差值，相互作用类型则由电子密度分布、分子内键响应和表面配位共同限定。

TQBQ-COF 孔腔中的三种药物分子保持完整骨架，框架也未形成新的局域配位。此类构型给出非共价主客体吸附的候选状态，归属仍由后续实空间分析确认，避免把“距离较远”直接写成物理吸附。

非共价区域和表面新键留下哪些信号？

NCI/RDG 直接分析低密度区域的梯度特征。在 sign(λ₂)ρ 接近零的位置，绿色等值面常对应色散主导的弱相互作用；负值区的蓝色区域可来自氢键或较强静电吸引；正值区的红色区域对应位阻排斥。颜色定义与等值面阈值取自论文图注，跨论文比较时记录相同绘图参数。

药物/TQBQ-COF 复合体在孔壁与吸附分子之间出现大面积绿色区域，弱色散接触覆盖了较宽空间。蓝色斑块对应局域吸引，红色尖峰来自芳香骨架内部或拥挤部位的排斥。大面积绿色区域与完整分子构型相互对应，物理吸附归属由结构和密度梯度两类信息共同支撑。

QTAIM 把电子密度划分为原子盆，并在键临界点 BCP处读取 ρ(r)、∇²ρ(r) 和总能量密度 H(r)。低 ρ(r)、正拉普拉斯与 H(r)≥0 常见于闭壳层非共价接触；H(r)

巯基吡啶/Al₁₂P₁₂ 的键路径与 NCI 区域同时存在，部分共价和静电成分可以共存。此类状态处在连续谱中，命名时写明“配位化学吸附”“部分共价吸附”或“非共价作用占主导”，比强行压成两个互斥标签保留更多物理信息。

轨道/振动/分子内键怎样约束吸附归属

化学吸附形成新配位时，吸附物前线轨道与表面态会产生能量对齐、杂化展宽和占据变化。PDOS 中的峰重叠只能定位可能参与相互作用的轨道；成键/反键贡献还由 COHP/ICOHP、键长与占据共同约束。差分电荷密度适合定位积累和耗尽，Bader 电荷给出原子盆积分，两者不承担化学键级的直接测量。

CO₂/TiO₂ 的化学吸附态中，Ti d 电子进入 CO₂ 的 2π^* 反键轨道，谱峰位置和占据随之改变；C–O 键由 1.18 Å 拉长到约 1.26 Å。反键态占据增加与分子内键拉长构成同向证据，单纯的电荷转移数值则无法区分极化与成键。

局域势能面的曲率决定振动频率，界面作用改变曲率后便会移动相应峰位。物理吸附常保留吸附物原有振动模，只出现静电场、氢键或色散环境造成的位移；表面新键、配位改变和反键态占据会带来新振动模、峰分裂或较大的频移。频率计算还负责确认优化构型没有非目标虚频，过渡态则要求沿反应坐标的一条虚频。

Al₁₂P₁₂ 吸附巯基吡啶后出现构型相关的 IR/Raman 峰形与强度变化。频移方向依赖具体振动模和电子回馈，红移并非化学吸附的通用标签。可靠归属需要追踪振动本征矢，确认变化来自吸附物内部键、表面—吸附物键或整体平移/转动受限。

怎样处理阈值/局域极小/解离路径的冲突

固定的 0.5 eV、0.8 eV 或某个吸附距离阈值只适合特定数据集的经验筛选。参考态改变会平移 E_ads：完整分子、自由基、解离片段和凝聚相化学势对应不同零点；大分子色散累积、离子吸附和氢键网络也会跨越经验能区。分类时保留能量数值及其定义，标签由相互作用机制承担。

半局域泛函常低估长程色散，未加色散修正的优化可能丢失物理吸附势阱；某些色散方案又会改变平衡高度和吸附能。强关联表面的 U 值、自旋态、表面电荷与溶剂模型会移动化学吸附态能量。方法敏感性应作为分类结果的一部分，例如记录“PBE 无稳定弱吸附态，PBE-D3 出现浅势阱”。

完整分子吸附态与解离产物之间若存在势垒，二者是势能面上的两个状态。CI-NEB、dimer 或 LST/QST 可追踪最小能量路径，反应能与活化能分别描述端点差和路径最高点。几何优化中直接断键只说明该初始构型落入解离产物盆，无法给出从气相分子到产物的动力学难度。

CO₂/TiO₂ 的六种初始吸附构型通向 HCOO^• 时呈现不同曲线：部分路径出现最高点，部分路径持续降能。吸附态的局域配位决定后续成键是否还要跨越势垒，物理吸附态也可能是反应前驱体。DFT 归属记录应写出候选构型、E_ads 定义、NCI/QTAIM 或 COHP 信号、分子内键和振动响应，以及 NEB 或过渡态频率验证。