一句话总结

DensityTool 是一款 VASP 后处理工具，能够从部分电荷密度和能量本征值中构建空间分辨和自旋分辨的局部电荷密度与自旋密度。与原子球投影相比，该方法提供了更优越的空间分辨率，同时支持对界面和缺陷等非均匀体系进行直接实空间可视化以及平面平均态密度绘图。

核心贡献

• DensityTool 处理 VASP 计算得到的部分电荷密度和能量本征值，用于构建异质固体材料的局部电荷密度与自旋密度。
• 该方法的空间分辨率较原子球投影有显著提升，并支持直接实空间可视化与平面平均局部态密度绘图。
• 该工具能够对界面、缺陷、表面、吸附分子以及混合无机-有机复合材料等非均匀体系中的局域电子结构进行系统性操作与可视化。

引言

电子结构计算是计算材料科学的基础，其中态密度是预测磁性、催化和输运性能的关键描述符。尽管 VASP 仍是此类模拟的标准计算引擎，但提取空间分辨和自旋分辨的电子分布通常依赖于自定义脚本或零散的工具，这导致了可重复性瓶颈并拖慢了分析流程。研究团队提出 DensityTool，这是一个专用的后处理框架，能够直接从 VASP 输出文件中自动化提取和映射空间分辨与自旋分辨的态密度。通过标准化该提取流程，该工具消除了手动编码的开销，使研究人员能够快速表征局部电子环境，从而加速先进材料的设计。

方法

研究团队利用 VASP 计算输出，特别是依赖于能带和波矢的部分电荷密度 $P_{n,\mathbf{k}}(\mathbf{r})$Pn,k​(r) 及对应的能量本征值 $\epsilon_{n,\mathbf{k}}$ϵn,k​，来构建局部态密度（LDOS）和局部自旋态密度（LSDOS）。该方法的核心在于对布里渊区内的这些量进行积分，以获取空间分辨的电子结构信息。在给定能量 $E$E 和位置 $\mathbf{r}$r 处的 LDOS 计算公式为 $L(E, \mathbf{r}) = \frac{N_e}{(2\pi)^3} \sum_n \int_{\text{BZ}} \delta(E - \epsilon_{n,\mathbf{k}}) P_{n,\mathbf{k}}(\mathbf{r}) \, d^3k$L(E,r)=(2π)3Ne​​∑n​∫BZ​δ(E−ϵn,k​)Pn,k​(r)d3k，这实质上代表了能量分辨的部分电荷密度。对于磁性体系，自旋分辨量的计算方式类似，自旋向上和自旋向下分量分别处理，从而得到自旋密度 $s(\mathbf{r}) = \rho^\uparrow(\mathbf{r}) - \rho^\downarrow(\mathbf{r})$s(r)=ρ↑(r)−ρ↓(r) 以及 LSDOS $S(E, \mathbf{r}) = \frac{1}{(2\pi)^3} \sum_n \int_{\text{BZ}} \left[ \delta(E - \epsilon_{n,\mathbf{k}}^\uparrow) P_{n,\mathbf{k}}^\uparrow(\mathbf{r}) - \delta(E - \epsilon_{n,\mathbf{k}}^\downarrow) P_{n,\mathbf{k}}^\downarrow(\mathbf{r}) \right] d^3k$S(E,r)=(2π)31​∑n​∫BZ​[δ(E−ϵn,k↑​)Pn,k↑​(r)−δ(E−ϵn,k↓​)Pn,k↓​(r)]d3k。

框架图展示了计算混合无机-有机体系平面平均 LDOS 和 LSDOS 的流程，其中部分电荷密度首先在与 Si 表面平行的平面上进行平均。该平均步骤通过 PARCHGSPIN 例程执行，降低了数据维度，从而借助 LDOSMAG 和 LSDOSMAG 例程高效计算平面平均 L(S)DOS。所得数据可作为能量和位置的函数进行可视化，有助于深入理解电子态的空间分布。该方法专为兼容 VASP 输出格式而设计，允许将完整的 L(S)DOS 写入 CHGCAR 格式，以便在 VESTA 等工具中直接可视化。

实验

DensityTool 的评估采用了非均匀材料，包括复合体系和钙钛矿平板模型，以验证其从复杂层状系统中提取有意义电子结构信息的能力。通过计算平面平均局部态密度，分析结果表明电子特性紧密遵循交替原子层的排列方式。定性来看，带隙附近的类体态主要源自 PbI₂ 层，而 PbI₂ 终止的表面引入了独特的表面态，这些表面态使整体带隙变窄，且与 MAI 终止的对应表面相比表现出更高的电离势。这些结构与电子特性的关联证实了该工具在表征异质材料方面的有效性，并与公认的理论预期相一致。

研究团队阐述了将 DensityTool 应用于非均匀体系分析的过程，重点聚焦于钙钛矿平板模型。该工具提供了计算电荷和密度分布的例程，并针对磁性体系利用自旋分辨数据进行了特定实现。钙钛矿平板的计算结果显示，电子结构因表面终止方式的不同而呈现显著差异，MAI 终止与 PbI2 终止的表面之间观测到了表面态和带隙的差异。DensityTool 为非磁性和磁性体系均提供了例程，涵盖部分电荷与自旋计算。钙钛矿平板的结果表明，表面终止对电子结构具有显著影响，PbI2 终止的表面显示出更深的表面态和缩小的带隙。平面平均局部态密度表明，带隙附近的类体态主要由 PbI2 层构成，而表面态则受暴露的表面终止方式影响。

研究团队运用 DensityTool 结合电荷、密度及自旋分辨计算的例程，对非均匀钙钛矿平板模型进行了分析。该设置验证了表面终止如何调控局部电子特性，结果表明与 MAI 终止构型相比，PbI2 终止的表面具有更深的表面态和更小的带隙。定性结果表明，带隙附近的类体态源于内部的 PbI2 层，而表面态则直接受暴露的终止化学性质控制。最终，分析证实该工具能够可靠地捕捉依赖于终止方式的电子变化，凸显了表面成分在塑造钙钛矿电子结构中的主导作用。