一维纳米线模型怎么建？

引言

一维纳米线和二维 slab 的边界条件完全不同。slab 通常在一个方向留真空、两个方向保持周期；纳米线只沿轴向周期重复，横向两个方向都要隔开周期镜像。模型如果没有把这件事处理清楚，后续带结构、磁矩、吸附能和电势都会混入相邻镜像线之间的相互作用。一维纳米线建模的核心，是把轴向真实周期和横向真空隔离分开检查。

本文围绕三个实操问题展开：怎样定义纳米线轴向晶胞，横向真空怎样收敛，轴向 k 点和横向偶极怎样检查。示例参数只作为起点，金属线、半导体线、磁性线和带吸附物的线状模型都需要单独测试。

图1 M₆X₆ 一维原子线的结构、声子谱和不同过渡金属硫族原子线模型。图中给出了一维周期结构的轴向重复单元和横向孤立特征。DOI: 10.1021/jacsau.0c00049。

纳米线晶胞先确定哪个方向？

建一维模型时，第一步不是加真空，而是确定真实周期方向。若纳米线沿 z 方向延伸，c 轴应对应线方向的重复长度；a 和 b 方向主要用于横向隔离。这样设置后，KPOINTS 才能自然写成 1 × 1 × N_k，结构图、能带路径和电势检查也能对应到同一个方向。

Δk_z ≈ 2π / (N_k c) 式(1)

式(1)说明轴向 k 点间隔由周期长度 c 和采样数 N_k 共同决定。c 越短，轴向布里渊区越大，需要更多 k 点；若把纳米线方向放错，或者在横向也给出多个 k 点，就会把真空方向的无意义采样引入计算。一维模型的 KPOINTS 重点在轴向，横向通常保持 1 个点。

对从体相切出的纳米线，还要确认端面是否已经被周期连接。若结构本身沿轴向有二聚化、螺旋、磁序或 Peierls 畸变，最小重复单元可能不是体相晶胞的一个周期。此时应比较 1 倍、2 倍甚至更长轴向单元的总能、键长和磁矩，而不是只沿用最短晶胞。

横向真空怎样判断够不够？

横向真空的任务是削弱相邻镜像线之间的静电、范德华和波函数尾部相互作用。若纳米线直径为 D，晶胞横向长度 a 和 b 应给出足够的空区，使相邻镜像线的最近距离 L_vac 达到收敛要求。对极性线、带电线或强偶极吸附态，所需横向距离往往比普通非极性线更大。

L_vac = min(a,b) − D 式(2)

式(2)只是几何估算，真正的判断要看目标量。可从 12 Å、15 Å、18 Å、22 Å 等横向尺寸做单点测试，比较总能每原子、带隙、磁矩、吸附能和 LOCPOT 横向平台。若横向真空改变后目标能量差仍明显漂移，说明镜像线之间仍存在可观相互作用。

图2 Fe₆S₆ 原子线的 AIMD 能量、线间距相互作用能和分子反应能量图。图中线间距能量曲线可用于理解横向镜像相互作用。DOI: 10.1021/jacsau.0c00049。

横向真空也会影响计算成本。a、b 过大时，平面波格点数量和 FFT 成本会升高，但它不会像增加轴向超胞那样直接改变物理周期。因此更合理的做法是：先用较低精度找出不会发生镜像干扰的范围，再在最终模型上提高 ENCUT、电子收敛和轴向 k 点。

轴向 k 点怎样收敛？

纳米线的能带沿轴向展开，k 点不足会直接影响费米能级、带隙、磁矩和吸附能。半导体纳米线可以从 1 × 1 × 5、1 × 1 × 7、1 × 1 × 9 开始；金属线或窄带磁性线应使用更密轴向网格，并保持 ISMEAR 和 SIGMA 口径一致。

ΔE_ads(N_k) = E_ads(N_k) − E_ads(N_k + ΔN) 式(3)

式(3)提醒 k 点收敛要看吸附能差，而不是只看纳米线本身总能。若 H、O、CO 或 OH 吸附改变了费米能级附近态密度，干净线通过的 k 点不一定适合吸附态。干净纳米线、吸附态和过渡态应采用可比较的轴向 k 点口径。

图3 ZnO 纳米线和 ZnO@CNT 复合结构的能带与态密度。图中轴向能带说明一维模型的 k 点采样会直接影响电子结构判读。DOI: 10.1038/s41598-018-33991-x。

如果纳米线轴向超胞被拉长，例如为了放置低覆盖度吸附物或缺陷，k 点数可以相应减少，但等效 k 点密度要保持可比。一个常用做法是先用原始轴向周期收敛 k 点，再按轴向长度比例折算到更长超胞，并用静态单点确认目标量差异。

横向偶极和电势怎样检查？

非对称吸附、极性表面或包覆结构可能让纳米线横向出现净偶极。此时即使真空足够大，LOCPOT 横向平均电势也可能显示缓慢漂移。VASP 的偶极修正可作为诊断工具，但不能代替横向真空收敛；尤其是一维线状体系，横向两个方向都可能有分量，单一 IDIPOL 方向不一定能描述全部情况。

更可靠的检查方式是同时看结构对称性、横向电势平台、差分电荷密度和真空尺寸收敛。若只在 x 方向偏心吸附，可把纳米线居中、增大 a 方向真空，并测试 IDIPOL = 1；若同时存在 x 和 y 方向不对称，则应优先扩大横向超胞并检查不同放置方向是否改变能量。偶极修正只能减少特定方向的周期电势误差，不能修复过小真空或不合理构型。

图4 ZnO 纳米线封装于碳纳米管中的俯视和侧视结构。图中展示了线状模型的横向截面、轴向方向和包覆环境。DOI: 10.1038/s41598-018-33991-x。

带电纳米线更难处理，因为周期性背景电荷和横向真空会共同影响电势零点。若不是研究明确的带电态，优先使用中性模型；若必须改变电荷数，应报告真空收敛、电势参照和电荷修正口径，并避免把绝对总能差直接解释为吸附热力学。

VASP 输入可以怎样起步？

一个中性半导体纳米线的起始设置可以采用 1 × 1 × N_k 的 Gamma-centered 网格，横向 a、b 给出 15 Å 以上空区，并按总能差、带隙、磁矩和吸附能继续收敛。若体系含过渡金属或开壳层吸附物，应开启自旋并测试 MAGMOM 初值。

KPOINTS 参考示例 Gamma 0 Gamma 1 1 9 0 0 0  INCAR 参考示例 ENCUT = 500 EDIFF = 1E-5 EDIFFG = -0.02 IBRION = 2 NSW = 150 ISPIN = 2 ISMEAR = 0 SIGMA = 0.05 LREAL = Auto

如果纳米线是金属性，ISMEAR 可改用适合金属的展宽设置，并单独测试 SIGMA；如果是强局域 3d 体系，还要考虑 DFT+U 与磁序。纳米线模型最终能否用于吸附和催化判断，取决于横向真空、轴向 k 点、磁性初值和结构周期是否同时通过检查。

图5 TiI₃ 纳米线中 1、3、7 链模型的结构示意。图中体现了横向截面大小会改变一维模型的物理身份和周期镜像距离。DOI: 10.1039/d2na00877g。

对吸附计算，可先在干净纳米线上确认横向真空和轴向 k 点，再放入吸附物测试不同轴向覆盖度。若吸附物沿线方向相隔太近，得到的是高覆盖度周期阵列；若横向真空不足，得到的是相邻纳米线共同作用后的假模型。把这两类误差分开，后续吸附能和能垒才有可比性。

本章要点总结

要点：一维纳米线只有轴向是真实周期，横向两个方向主要用于隔离周期镜像。

要点：横向真空应通过能量、磁矩、带隙、吸附能和电势平台共同检查。

要点：轴向 k 点直接影响一维能带、费米能级和吸附能差，不能照搬 slab 网格。

要点：横向偶极修正只能作为辅助诊断，不能替代足够大的真空和合理构型。