ms理论计算

我刚接触 ms理论计算 时，最不习惯的是它的”任务流”思维——所有计算都挂在 Study 树上，参数散落在不同模块里。第一次算一个钙钛矿的能带，我在 CASTEP 里忘了开自旋，band gap 算成金属，排查半天才发现是 Setup 页的 spin 没勾。那次让我明白，ms理论计算的坑大多在模块间的参数一致性，而不是方程本身。

一、为什么 ms理论计算首选 CASTEP

MS 里做周期体系我默认 CASTEP（平面波赝势 DFT），做孤立分子则换 DMol³（原子轨道基）。CASTEP 的好处是工作流完整：Build 晶体 → 切表面 → 几何优化 → 能带/DOS/光学/弹性一站式。我算二维材料剥离能、表面吸附都用它。做 ms理论计算 时，赝势选超软（ultrasoft）还是 PAW 要心里有数：PAW 精度高但慢，超软快但极硬体系要谨慎。

二、核心原理：平面波 DFT 在 MS 里怎么跑

ms理论计算 的 backbone 是 Kohn-Sham 方程在平面波基组下求解。关键参数三件套：截断能（Energy cutoff，我常用 400~500 eV）、k 点网格（按体系尺寸取 0.03 Å⁻¹ 间距）、收敛容差（力 0.01 eV/Å、能量 1e-5 eV/atom）。曾经一个氧化物因 cutoff 只用 300 eV，晶格常数偏小 1.5%，调到 450 eV 才稳。这三件套不做收敛测试，后面所有性质都悬空。

三、关键技术要点：几何优化先于一切

很多人急着算性质，跳过了充分几何优化。我做 ms理论计算 时，优化阶段就把自旋、DFT+U、色散都设好，让结构先松弛到位，再在同一精度下算能带和 DOS。一个金属有机框架，优化未收敛就跑性质，孔径算错导致吸附位点判断全反。优化收敛标准我也卡得比默认紧——位移 1e-4 Å、力 0.01 eV/Å，避免”看着收敛其实没到”。

四、实操流程：从建模到性质

我的标准动作：导入/搭建晶体 → 加对称性检查 → 切表面设真空层 → CASTEP 几何优化（开自旋/色散按需）→ 提高 k 点做性质计算（Bandstructure + DOS + 需要的 Optical/Elastic）→ 用 Analysis 工具出图。做ms理论计算最容易被省的是性质计算的 k 路径，能带图若路径不全，会漏掉真实带边位置，我习惯按 Brillouin 区高对称点点满。

五、常见踩坑：任务参数没继承

ms理论计算 里最常见的事故是优化用一套参数、性质计算又改了另一套，导致结构与分析不自洽。我一律在优化后再”计算性质”而非新建任务，保证继承结构。另一个坑是 DMol³ 算分子忘了设电荷和多重度，开壳层当闭壳层，轨道能全错。还有溶剂化（COSMO）忘了开，液相体系直接气相算。

六、复盘：MS 是”平台”不是”黑箱”

回过头看，ms理论计算 真正的门槛不是懂 DFT，而是把 MS 的工作流参数串成一条自洽链。我接项目时默认先要体系类型（晶体/表面/分子）和目标性质，再选模块定参数。把赝势、截断能/k点、几何优化三件事做规整，CASTEP 和 DMol³ 的输出才敢写进报告。被证明有用的，是愿意为一次收敛测试多等半小时的耐心。下次有人催”能不能快点”，我会说——省掉的那次收敛测试，后面要用三天排查来还。

一个 TiO₂ 的例子：默认 k 路径只走了 Γ−X−M，漏了 Γ−Z 方向，带边被误判在 X 附近，补上 Z 之后才发现真实导带底在 Γ−Z 中间，带隙大了 0.3 eV。这种漏路径的错最容易发生在低对称表面，我现在出带隙前一定先画完整 BZ 路径核对。MS 的 k 点路径是可以手动编辑的，别全信自动生成的——自动路径有时会跳过你最关心的那个高对称点。