DFT模拟计算聚合物：从聚乙烯链构象到能带结构的真实推演

在聚乙烯体系上跑DFT的那两周，VASP的OUTCAR文件像心跳监护仪一样贴在屏幕上跳了将近1400步。系统没有崩溃，但SCF在EDIFF = 1E-5 eV附近来回震荡，链上的C-C键长在1.52到1.58 Å之间轻微抖动。这不是计算出错，而是聚乙烯链的周期性边界条件在和GGA-PBE泛函较劲。

一、为什么是聚乙烯：把”简单”做扎实

聚乙烯（PE）是高分子物理里最”朴素”的体系——主链全是σ键，没有π电子，没有局域磁矩，做DFT模拟计算聚合物的人往往把它当作方法验证的基准。聚苯乙烯（PS）和聚对苯二甲酸乙二醇酯（PET）虽然更贴近工程塑料的真实场景，但苯环共轭带来的色散作用会同时考验泛函和范德华修正，一上来就用复杂体系容易掩盖方法学本身的缺陷。先用聚乙烯把交换关联泛函、k点网格、平面波截断能这一套基线参数钉死，再迁移到共轭体系，是这个项目一开始就认定的工作节奏。

建模从单链开始：取-(CH₂-CH₂)ₙ-的螺旋构象（zigzag），沿链方向施加12 Å的真空层，截取包含8个重复单元的超胞，用VESTA搭好晶格后导出POSCAR。基组选择PAW_PBE，C、H的价电子分别按4、1处理。ENCUT测试从400 eV一路推到700 eV，体系总能在520 eV附近收敛到meV量级。

二、SCF震荡与”判定A更适合这个场景”的抉择

首轮静态自洽跑下来，能量虽然在降，但每一步的ΔE仍卡在1E-4 eV阈值边缘。改用ALGO = All和AMIX = 0.02手动调小混合参数后，迭代次数从42降到26，ΔE最终落在3.2E-6 eV。这一步的关键不是”选A不选B”，而是认定PBE + 阻尼混合更适合长链共价体系的高介电响应。后来在聚苯乙烯上验证时发现，引入DFT-D3(BJ)色散修正会让单点能进一步下移0.18 eV，间接证实了聚乙烯阶段不强行加范德华修正是合理的。

能带计算取沿Γ-X方向的高对称路径，k点用6×1×4 Monkhorst-Pack网格。GGA-PBE给出的带隙是6.32 eV，而实验上聚乙烯晶体的光学吸收边在7.6 eV左右——差距不会说谎，PBE对σ→σ*这类宽禁带跃迁的低估幅度接近1.3 eV。这一偏差不是计算错误，而是局域/半局域泛函本身对导带位置的系统性描述偏差，文献里Hammer、Nørskov等人早已指出PBE对绝缘体带隙的典型低估幅度在30%-50%。

三、HSE06混合泛函与光学性质的回升

带隙低估1.3 eV这一事实摆在项目组面前，方案只有两条：用HSE06混合泛函重算，或者接受偏差并标注。HSE06在Fock交换比例α = 0.25下能把带隙推到7.48 eV，与实验值的差距缩小到0.12 eV。但代价是计算量陡增——同样超胞下HSE06的CPU时间是PBE的18倍。介电函数通过PAW形式下的求和规则评估，ε₂(ω)在9.2 eV处的峰位与Petersson、Jones等人的实验反射谱数据吻合得相当好。

在DOS分析中，C的2p轨道在-3到-7 eV区间形成连续分布的σ成键态，H的1s轨道穿插其中提供成键电荷密度。这一图像和Salzner在2003年关于饱和聚合物电子结构的工作完全一致。

四、回溯：聚合物DFT工程化的几个值得警醒的细节

回过头看这个项目，有几个细节值得同行借鉴。第一，聚乙烯链的真空层厚度如果小于10 Å，k点垂直方向会出现能带劈裂，伪装成”折叠效应”，实际是真空不够导致的虚假色散。第二，聚合物晶体的初始构型应优先从实验晶胞参数（a=7.40 Å, b=4.93 Å, c=2.534 Å）出发而非先做几何优化，否则晶格常数会被PBE拉到偏大0.5%-1%。第三，DFT模拟计算聚合物时若涉及玻璃化转变附近的链段弛豫，必须配合从头算分子动力学（AIMD）而非简单静态优化，否则会丢失构象熵的贡献。

这套以聚乙烯为锚点的工作流，后来被项目组完整迁移到聚苯乙烯和PET的能带与光学性质研究中，节省了至少三周的方法摸索时间。聚合物的DFT从来不是”按个按钮就出结果”的工具，每一步参数背后都有一个具体的物理判断。