表面吸附计算要不要保留对称性ISYM？

引言

表面吸附计算里，ISYM 最容易被忽略：干净 slab 看起来有对称性，吸附物一放上去，有时还保留默认 ISYM=2；结果优化后吸附物被拉回某个高对称位置，两个本该不同的原子受力被平均，或者电荷密度被强行恢复成较高对称性。吸附态不是 bulk，小分子、缺陷、覆盖度和偶极方向都会主动破坏表面对称。

VASP 官方说明里，ISYM 控制对称性处理；ISYM=1、2、3 会使用对称性，ISYM=0 或 -1 会关闭对称性。打开对称性时，VASP 会按检测到的空间群对总电荷密度、离子受力和应力张量做对称化，并约化 k 点。这个机制对干净晶体很有用，但对非对称吸附构型，可能把真实的局域畸变平均掉。ISYM 的关键不是省多少时间，而是它是否改变了你允许体系走向的构型空间。

图1 非对称与对称 slab 构型示例。上下表面终止方式、分子修饰和真空方向偶极，都会影响 ISYM 是否适合保留。DOI: 10.1080/14686996.2017.1370962。

ISYM在结构优化里到底做了什么？

它不只是减少 k 点

很多人把 ISYM 理解成加速开关，因为对称性可以减少不可约 k 点数量。但在结构优化中，它还会让电荷密度和力满足检测到的对称操作。若两个表面原子被识别为等价，它们的受力会被对称化；若某个局域电荷分布破坏了对称性，电荷密度也可能被拉回对称形式。当你希望体系自发发生局域畸变、成键偏转或电荷局域时，对称化本身就是一个约束。

F_i^{sym} = (1 / N) Σ_g g F_g^{-1}i 式(1)

式(1)把受力平均的逻辑写成了一个简化图像：被对称操作连接的原子，其力会按对称关系处理。bulk 弛豫中，这种处理往往符合晶体本身的约束；到了单侧吸附、缺陷表面、局域自旋态或有偶极的 slab 里，它也可能削弱真正推动构型降低对称性的那部分力。

SYMPREC 为什么会让结果变得微妙？

SYMPREC 决定位置差异多大时仍被判断为等价。表面模型来自建模软件、坐标精度有限，或优化初期构型还很接近高对称位置时，VASP 可能识别出并不适合吸附态的对称操作。把 SYMPREC 调得很严可以减少误判，但不能替代物理判断：这个吸附态本身是否允许破坏对称性。

图2 对称终止 Si(111) slab 的平面平均电荷密度示例。slab 或吸附层保持对称时，ISYM 能提高效率；上下表面或吸附物不同，模型身份已经变了。DOI: 10.1080/14686996.2017.1370962。

为什么吸附态容易破坏表面对称？

单侧吸附天然降低 slab 对称性

多数表面吸附模型是在 slab 上表面放一个吸附物，下表面固定或钝化。这个模型从几何上已经区分了上、下表面，若再加入 LDIPOL、IDIPOL 或表面电荷转移，真空方向更不应被当作对称等价。单侧吸附优化里，ISYM=0 往往能让构型自由弛豫；若体系含非共线磁性、自旋轨道，或需要完全关闭时间反演相关约化，再把 ISYM=-1 纳入对照。

对称性还会在面内被吸附物打破。例如一个 CO 只吸在某个顶位，一个 OH 倾斜成键，一个 OOH 横跨桥位，这些构型本来就带着方向性，不该被强行拉回高对称方向。吸附物的取向和局域成键，就是优化需要自由寻找的结果。

图3 Al 终止与 O 终止 α-Al2O3(0001) 表面弛豫前后结构。不同终止和表面弛豫会降低对称性，吸附后这种差异通常更明显。DOI: 10.1021/acsomega.1c03771。

覆盖度和取向会让等价位点变成不等价位点

同一个吸附物在低覆盖度、高覆盖度、平躺、直立、桥式配位之间切换时，面内对称性都会变化。若保留 ISYM，VASP 可能把几个本来因吸附而不等价的表面原子继续看成等价，导致能量排序和优化路径偏离真实低对称构型。吸附态的对称性应该由最终构型自己表现出来，而不是由初始 slab 的高对称性预先规定。

*OOH、*COOH、*CHO 等催化中间体常带有明确方向，H 键网络也会选定某一侧。比较不同取向的能量时，这些取向需要独立优化；一旦在对称操作下被折叠成同一种构型，表面上省了计算，实质上丢掉了可比较的构型空间。

图4 EC 分子在 α-Al2O3(0001) 表面的不同取向吸附能。吸附取向会改变局域对称性和能量排序，是关闭 ISYM 做对照的典型场景。DOI: 10.1021/acsomega.1c03771。

什么时候可以保留对称性？

干净 slab、bulk 和明确高对称参考态

计算对象是 bulk、干净对称 slab、表面能或高对称参考结构，且没有单侧吸附、缺陷、外电场、偶极修正和局域磁性破缺时，默认 ISYM 保留下来并不突兀。它能减少 k 点和加速计算，也能让等价原子保持应有关系。对严格高对称路径的能带或声子前处理，对称性还承担着维持模型定义的作用。

干净 slab 能保留对称性，不代表吸附态也能沿用。更清楚的口径是分开处理：clean slab 优化保留 ISYM；吸附态、缺陷态、带局域磁矩态优化采用 ISYM=0；最终静态计算若读取已优化的低对称结构，也沿用同一口径。优化和静态单点的 ISYM 不要无说明地切换。

覆盖度模型若本身就是对称吸附层

完整有序单层是另一种情况：每个表面位点都吸附同一种分子，吸附层与 slab 周期严格匹配，模型本身就定义为有序相。此时保留一部分对称性可以服务于这个有序模型；若目标转向局域缺陷、空位或分子转向，对称性就不再是加速工具，而会变成限制。

E_ads = E_slab+nX − E_slab − nE_X(g) 式(2)

式(2)用于覆盖度吸附能。若 n 个吸附物由对称操作严格连接，ISYM 可帮助维持有序覆盖度模型；若 n 个吸附物可能发生不同倾斜、解离或氢键重排，就应关闭对称性，让它们独立弛豫。

哪些体系不该被对称性束住？

单侧吸附、缺陷、局域电荷和磁性

单侧 slab 吸附、含氧空位、掺杂、台阶、边界或低覆盖度吸附，都已经把原来的表面对称性拆开了。吸附物还可能倾斜、旋转、转移质子或解离；体系也可能出现局域磁矩、电荷极化或 LDIPOL 修正。迁移、NEB 和 AIMD 更依赖低对称路径展开，VASP 官方也把 ISYM=0 放在分子动力学相关场景中使用，因为动力学本身不应被静态对称操作限制。

对表面反应，关闭对称性通常不是为了追求“更低能”这个单一目标，而是为了允许真实低对称路径出现。若关闭后结构大幅偏离，先检查是否合理：吸附物是否滑到相邻位点、是否发生解离、固定层是否过少、真空方向是否设置正确。ISYM=0 给体系自由度，但不替代结构合理性检查。

图5 EC 在 α-Al2O3(0001) 不同覆盖度下的电荷流动。覆盖度和取向改变会引入方向性电荷重排，保留对称性前应确认模型本身确实有序。DOI: 10.1021/acsomega.1c03771。

ISYM=0 和 ISYM=-1 怎么选？

一般表面吸附优化中，ISYM=0 已经能关闭常规空间对称性；它仍会利用 Ψ_k=Ψ*_-k 这样的关系来约化布里渊区采样。ISYM=-1 则完全关闭对称性。非共线磁性、自旋轨道、外场，或任何不希望出现对称相关约化的任务，都值得把 ISYM=-1 单独拿出来测试。普通共线自旋吸附优化则可以先从 ISYM=0 建立一致口径。

VASP里怎样做对照？

先用同一结构做静态对照

同一个已优化结构分别做 ISYM 默认值和 ISYM=0 静态单点，能把问题拆得比较清楚：能量、力、磁矩和电荷如果几乎不变，对称性没有明显改写该结构的能量口径；差异一旦放大，后续数据就应回到关闭对称性的统一口径。

ΔE_sym = E_ads(ISYM=0) − E_ads(ISYM=2) 式(3)

式(3)诊断的是 ISYM 对吸附能的影响。两个计算若走向不同优化路径，ΔE_sym 同时包含“不同终态结构”和“不同电子口径”两部分；因此先做同一结构静态对照，再比较各自充分优化后的最终构型，结论才不容易混在一起。

SYSTEM = surface_adsorption_isym_check ENCUT  = 520 EDIFF  = 1E-5 EDIFFG = -0.03 IBRION = 2 NSW    = 160 ISIF   = 2 ISYM   = 0 LDIPOL = .TRUE. IDIPOL = 3 # clean slab 单独测试默认 ISYM；吸附态优化从 ISYM=0 起步更清楚

结果部分只需把模型身份和计算口径交代清楚：干净 slab 采用默认对称性优化；吸附态因为单侧吸附破坏表面对称，结构优化和静态能量均采用 ISYM=0。若比较过默认 ISYM 与 ISYM=0，再给出关键吸附能差异和构型趋势。对称性设置要和模型物理身份一致，不能只沿用 bulk 的默认习惯。

本章要点总结

要点：ISYM 不只减少 k 点，还会按检测到的对称性处理电荷密度、离子受力和应力。

要点：单侧吸附、缺陷、低覆盖度、局域磁性、偶极修正、NEB 和 AIMD 通常更适合关闭对称性。

要点：干净 slab 或明确有序覆盖度模型可保留对称性，但吸附态应重新判断，不应无条件继承默认值。

要点：ISYM=0 与默认值的差别，最好用同一结构静态单点和充分优化结果两类对照分别判断。