手机版

DFT模拟计算聚合物：链状体系电子结构计算的实践与局限

发布时间：2026-06-15 来源：科研学术网

小中大

字号：

DFT模拟计算聚合物这件事，跟算晶体和算小分子都不一样。小分子体系没有周期性问题，晶体体系有明确的布拉伐格子和k点采样规则，而高分子链处于两者之间——有周期性但又不是三维的。处理这类体系时，周期性边界条件的设置方式直接决定了计算的物理合理性。

一维周期性：链方向怎么设

用VASP算聚合物，最自然的做法是把链方向设为周期性方向（比如z轴），在x和y方向留足够大的真空层。真空层的厚度需要保证相邻镜像链之间的相互作用可以忽略——对共轭聚合物来说，链间π-π堆叠的衰减距离大约在10-15 Å，真空层至少要设15 Å（x和y方向各15 Å）。

但真空层越大，平面波基组的计算效率越低。为什么？因为VASP的FFT网格是均匀划分的，真空区域里的波函数振幅接近零，但FFT仍然要在整个网格上操作。真空层从15 Å增加到20 Å，计算时间可能翻倍，但对物理结果的改善微乎其微。团队在做聚噻吩衍生物的计算时，测试了12 Å、15 Å和18 Å三种真空层厚度，发现HOMO-LUMO能隙的差别在0.02 eV以内——远小于泛函本身的系统误差。最终选了15 Å作为标准设置。

链长收敛性：重复单元够了没

聚合物计算中最核心的收敛性测试是链长——用多少个重复单元才能代表无限长链的性质？这个问题没有标准答案，取决于你关注的物理量。能隙对链长的依赖很强（短的寡聚物能隙大，随链长指数衰减趋近于聚合物极限），而几何参数（键长、二面角）在3-4个重复单元后就基本收敛了。

团队在聚3-己基噻吩（P3HT）的计算中，从1个到6个重复单元（nT到6T）系统地算了HOMO-LUMO能隙的变化。1T的能隙是3.8 eV，4T降到2.1 eV，6T是1.9 eV——从4T到6T的变化只有0.2 eV，外推到无限链的极限值大约是1.7 eV。实验上P3HT的光学吸收边在1.9 eV附近，考虑到DFT-PBE本身会低估能隙约30-40%，4T的计算值已经可以给出合理的趋势判断了。

这个链长收敛测试的计算量不小——每个重复单元对应一个独立的结构优化+电子步计算，6个重复单元就是6组完整的任务。但在不做收敛测试的情况下直接取一个”看起来够”的链长去算，风险是结果可能在一个没有收敛的链长上做分析，结论的可靠性无从谈起。

van der Waals力的处理：PBE的盲区

GGA-PBE泛函有一个已知的系统性缺陷：它不能描述van der Waals（色散）相互作用。对无机晶体来说这个问题不那么致命（主要是离子键和共价键），但对聚合物来说，侧链之间的堆积、链间的相互作用都依赖vdW力。如果不加修正，算出来的最优构型里侧链间距会偏大，密度偏低。

常用的修正方案有三种：DFT-D2/D3经验色散修正（Grimme）、vdW-DF（非局域相关泛函）、optB88-vdW（优化交换泛函）。对聚噻吩体系，团队试了PBE-D3和optB88-vdW两种。侧链间C-H…π距离的差异在0.1-0.2 Å，对总能量有约0.3 eV/重复单元的差别。不算大，但对研究堆积构型（face-on vs edge-on）的相对稳定性来说，这个差异足以改变结论。

选PBE-D3还是optB88-vdW？PBE-D3更轻量（计算量几乎不增加），optB88-vdW对能量描述更自洽但计算成本增加约30%。对于只需要定性判断堆积模式的场景，PBE-D3够了；如果要算绝对堆积能跟实验数据对标，optB88-vdW更合适。

k点采样和计算成本

一维周期性链的k点采样只在周期性方向上设置。一般来说，一个重复单元在链方向上取一个k点（Gamma点）就够了——因为重复单元本身的长度已经大于10 Å，Brillouin zone很小。但如果重复单元很短（比如聚乙炔的CH单元只有约1.4 Å），就需要在链方向上取更多k点。

计算成本方面，聚合物体系的成本主要取决于真空层和截断能。一个4重复单元的P3HT，真空层15 Å，ENCUT 400 eV，Gamma点，在32核上大约需要6-8小时完成结构优化。这个计算量对VASP来说不算大，但如果要做多组不同侧链取代基的对比（比如换烷基链长度、加极性基团），总计算量会快速叠加。