LOBSTER专题

LOBSTER键长分布、DOE与DOV

用Si两层近邻演示BWDF与BWDFCOHP，并说明辅助能量和轨道指标的读取方法。

更新于 2026-11-07约 3 分钟阅读按任务与环境选择适用配置

阅读本篇目录

除了能量分辨曲线，还可以按原子间距离整理成键权重。BWDF与BWDFCOHP是键权重相关分布；它们不是普通RDF，也不能在没有确认原生定义时直接改成“每原子键数”或“总键能”。

1 用同一Si波函数扩展近邻范围#

配套Si作业自动生成distance-distribution/。与最近邻例相比，lobsterin改用：

cohpGenerator from 0.1 to 4.1
BWDF 0.02
BWDFCOHP

其余能量窗、基组和VASP波函数保持不变。0.02是距离分箱步长，单位Å。BWDFCOHP需要同时启用BWDF，不能只写一行BWDFCOHP 0.02。在本版本中，缺少主开关的错误组合会在部分计算后终止；应修正输入后在新目录重新做后处理，不需要重跑已经成功的VASP静态计算。

head work/p100/si/distance-distribution/BWDF.lobster
head work/p100/si/distance-distribution/BWDFCOHP.lobster
head work/p100/si/distance-distribution/data/BWDFCOHP.csv

2 结果怎么读#

Si两层近邻的原生成键权重分布

实际几何壳层约2.35126和3.83959 Å，对应非零分箱中心约2.36126和3.84959 Å；二者有分箱中心偏移，不能把中心直接当成优化后的精确键长。两处BWDFCOHP原生权重约−4.56160246、−0.01884239，BWDF约**+0.37944866、−0.10437774**。

本例有28条原生周期原子对记录，包括第二层同编号原子的不同周期镜像。图保留原生权重，没有自行按原子数、重复次数或距离重新归一；原始CSV还保留了零距离占位行，图仅显示物理距离范围。

只选2.6 Å内单一近邻距离时，本组会得到仅一个零行的分布文件，虽然其他COHP文件正常。这样的空分布不能用来证明已经得到有效的键长谱。本篇采用两个实际几何壳层，能直接检查峰位置与结构对应关系。

3 DensityOfEnergy#

默认Si输入同时启用：

DensityOfEnergy

查看：

head work/p100/si/reference/DensityOfEnergy.lobster
head work/p100/si/reference/data/DensityOfEnergy.csv

这是局域投影框架下能量分布相关的辅助输出。先保留原生列和能量轴，并结合所用版本手册解释；不能直接把它当VASP总能，也不要把任意曲线面积称为形成能或反应能。形成能仍需相应结构及参照态的总能计算。

4 Density of Values：DOV#

启用：

writeDOV

原生文件名是DOSCAR.DOV.lobster。查看和导出：

head work/p100/si/reference/DOSCAR.DOV.lobster
head work/p100/si/reference/data/DOV.csv

可在配套环境中用Doscar(doscar="DOSCAR.DOV.lobster", structure_file="CONTCAR")读取。它是程序定义的另一种投影诊断量，不能看到DOSCAR字样就当作通常用于占据电子数检查的DOS；电子数核对仍使用DOSCAR.lobster。

本篇示例对应LOBSTER 6.0.0。扩展到其他元素、键长窗口和轨道组合时，先看原生实际输出、非零数据和几何对应关系，再决定如何作统计比较。

输入准备 · COHP与ICOHP · COOP、COBI与轨道 · DOS与电荷 · 基组与投影质量 · 键长分布与辅助输出 · 缓存与结果导出

需要根据你的体系选择资源？联系讨论计算需求。实际可用软件与权限，以当前账户和运行环境为准。