VASP 多尺度与实验接口操作导航
串联ATAT、CALPHAD、热力学、局域探针、衍射与连续介质的真实例子,说明哪些量可以直接传递。
阅读本篇目录
下载实际DFT数据、导出脚本和参考图,保存为 example.zip:
python3 -m zipfile -e example.zip .
cd vasp-multiscale-reference-examples
bash results.sh
bash plot.sh
本包复用本站已经完成的VASP数据,只做参照、单位和张量方向的整理,没有新增DFT或有限元计算。分析使用NumPy和Matplotlib,集群 ase/3.29.0 环境已提供所需依赖。
1. 根据下游任务选数据#
| 下游目标 | VASP或后处理提供什么 | 实际操作教程 |
|---|---|---|
| 合金簇展开、Monte Carlo | 一致的构型能量、晶格与ECI模型 | VASP到ATAT |
| CALPHAD数据准备 | 成分、构型、参照态、静态形成能 | Cu–Au数据导出 |
| 有限温度自由能 | 静态能、声子自由能及必要电子项 | QHA / Al电子与总自由能组合 |
| NMR | 绝对屏蔽及张量约定 | 水分子屏蔽 |
| EFG、核四极相关输入 | 电场梯度张量、核素参数 | Fe局域参数 |
| 磁超精细与穆斯堡尔相关输入 | 接触项、芯层贡献、所用自旋归一化 | 超精细 / Fe例子 |
| 结构或振动散射 | 结构因子、模态与几何约定 | 理论XRD / 声子散射 |
| 连续介质 | 弹性张量、密度和材料坐标 | 方向模量与SI导出 |
这些步骤复用各自真实的计算结果;相同的“能量”或“谱”字样,不代表文件可以互换。
2. 用终端先导出一份有参照的数据#
本包的 inputs 提供实际Cu–Au六构型、Si弹性、Fe局域量和水屏蔽结果。运行:
bash run.sh
cat results/CuAu-reference-export.csv
cat results/Fe-local-parameter-export.csv
输出保留构型、成分、单胞原子数和能量定义。替换自己的输入时,先确认所有参照计算来自兼容的方法、PAW与数值精度。
3. Cu–Au:每摩尔原子与每摩尔化学式不同#
先用单胞能量除以原子数,再按Au原子分数x减去 (1−x)eCu+x eAu。本例纯Cu、Au参照分别为−3.7302235975与−3.2221795000 eV/atom;换成J/mol原子使用96485.33212331。
等原子比的有序CuAu构型得到约−5447.55 J/mol原子,另一种同成分排布约−3951.73 J/mol原子。它们成分相同、排布不同,不能删掉构型身份再混成一个数。
1 mol CuAu化学式含2 mol原子,单位相差2倍。当前数据只是静态电子能,尚未拟合成TDB;真正热力学建模还需要相模型、温度项和统一参照,见CALPHAD数据教程。

4. NMR:屏蔽与化学位移要分开#
水分子实际绝对屏蔽中,两H均为30.6804 ppm、O为318.7559 ppm,采用不含G=0宏观项的约定。包含该项时,本例各原子的值均增加0.0531 ppm;不能把不同约定的数据直接混在一张标定表。
这些不是已经参照标定的实验化学位移。要给出位移,应选择并用相同方法计算合适参照,保留张量与宏观项约定。不能仅因为原始输出单位是ppm就直接标成实验位移。
5. Fe:保留EFG、接触项与核态定义#
本例Fe的三种同体积畸变,固定z方向EFG分别为−4.250、0、4.552 V/Ų。原生 A_tot不包含芯层接触贡献 A_1c,两者相加后分别为−18.062、−17.215、−17.514 MHz。
cat inputs/Fe-source-local-inputs.json
cat results/Fe-local-parameter-export.csv
MHz接触参数还涉及VASP的电子自旋归一化,不能直接改成T。该数据提供的是局域电子结构输入,没有计算同质异能移、激发核态四极谱或完整穆斯堡尔谱。后续谱模型要明确基态/激发态参数、参照、磁/四极耦合、取向和线宽。
6. 自己计算时保留哪些原始文件#
从对应教程进入自己的计算目录,先检查INCAR、POSCAR、KPOINTS,提交后保存方法所需的OUTCAR、结构、DOS或声子结果,再运行那一篇匹配的提取脚本。不要用一个通用脚本替代所有响应方法。
grep -E 'ISPIN|LEFG|LHYPERFINE|LCHIMAG|NGYROMAG' INCAR
grep -n 'Electric field gradients' OUTCAR
grep -n 'Fermi contact' OUTCAR
这些检查只对包含相关输出的计算目录有意义。本下载包提供可直接处理的数据与脚本;需要完整VASP输入和计算步骤时,回到表中对应的原生实算教程。