Posts

本文总结了Abinit结构优化的计算流程，以及代码的实现方法。

在利用 ABINIT + LOBSTER 进行 $\mathrm{Al}(111)$ slab 吸附体系的 COHP 计算时，发现所有 24 个任务在 SCF 收敛后进入长达约 12 小时的静默阻塞状态，无任何输出写入日志文件。本文记录该问题的诊断过程与修复方案。

将赝势从 norm-conserving（psp8，ta015）切换到 PAW（xml，PseudoDojo pbe/paw-sr-11，ta031）后，所有 COHP 计算的 SCF 迭代出现严重的收敛困难。能量在前几步快速下降后开始振荡上升，残差（vres2）在 $\sim1\times10^{-2}$ 处停滞不前。本文记录该问题的诊断与修复过程。

本文档记录使用 ABINIT 8.6.3 计算波函数，配合 LOBSTER 5.1.1 进行 COHP（Crystal Orbital Hamilton Population）分析的完整流程，包括版本选择理由、关键参数设置、输入文件模板、资源估算及常见问题排查。

等离子球（plasma globe）内部的发光"电弧"严格说更接近低压气体放电丝状通道，不是空气中高温金属电弧那种连续强电流电弧。它的形成涉及高频高压电场、气体击穿、电离、电子雪崩、等离子体鞘层和介质阻挡放电等过程。本文将逐一解释其中每一个关键公式的推导过程。

在使用 ABINIT 进行单原子参考能计算时，O 和 S 两种开壳层原子出现了 SCF 自洽迭代不收敛的问题。本文记录从问题发现到最终解决的完整排查过程，涵盖输入文件设置、根因分析和修正方案。

在使用 Python 绘制 DFT 吸附能曲线的过程中，发现某些元素的曲线图在特定 $z$ 坐标附近出现不符合物理预期的异常振荡。排查后定位到问题根源：CubicSpline 三次样条插值在面对尖峰数据时会产生严重的过冲（overshoot）。本文将记录这一问题的分析过程和修复方案。

在密度泛函理论（DFT）计算化学研究中，绘制 adatom（吸附原子）沿表面法线方向的吸附能曲线是分析表面催化性质的重要手段。本文介绍一个完整的 Python 脚本，利用 Matplotlib 将六种元素（N、H、C、Si、O、S）在四种吸附位点（Top、Bridge、fcc、hcp）上的 $\Delta E_{\mathrm{ads}}$ 曲线以 $3 \times 2$ 网格布局呈现，每个面板同时包含主图和局部放大图。