ASE命令行建模、计算与绘图
18个实际CLI入口的操作例,包含Cu经验势、真实Si能带、NEB和终端查看。
阅读本篇目录
下载完整脚本、实际输入、参考输出和四张高清图。在集群新解压目录执行 bash run.sh,结果写入 work/;使用ASE3.29.0、单进程和单线程,不调用GPU或新DFT计算。
1 不写Python也可以先开始#
ASE命令行适合快速建模、转换和查看;需要复杂流程时再写Python。本篇核心命令对应附件中的实际调用。先确认环境:
module load ase/3.29.0
ase --version
ase --help
ase help build
ase info
本地CLI提供18个入口。ase info --formats列出格式,ase info --all-calculators列出接口名称;接口列出不代表外部程序已安装,也不代表它适合当前材料。
2 建立32原子铜结构并转换#
在一个新目录执行,避免覆盖已有文件:
mkdir cu-cli
cd cu-cli
ase build -x fcc -a 3.6 --cubic -r 2,2,2 Cu cu32.traj
ase convert cu32.traj cu32.extxyz
ase info --files cu32.traj cu32.extxyz
-a 3.6的单位是Å,--cubic选择FCC常规胞,-r 2,2,2重复三方向。结果为32个Cu、7.2Å立方晶胞。这里把选项放在前面,最后两个位置参数 Cu cu32.traj相邻,避免命令行解析歧义。
3 用EMT做一个快速单点#
ase run emt cu32.traj --properties efs -o cu32-emt.extxyz
e、f、s分别请求能量、力、应力。EMT是内置经验势,不是DFT;本例完整32原子胞能量为 −0.214040597945eV,不是每原子能量。理想FCC对称结构的最大力接近数值零,这不等于晶格常数已经优化。
本篇小例用于掌握命令,材料研究应选择对应方法和收敛条件。VASP、ABACUS、CP2K和NEP实际联用分别见ASE接口教程。
4 用命令行完成状态方程#
另建一个Cu原胞,再对晶格长度±3%取九个点:
ase build -x fcc -a 3.6 Cu cu1.traj
ase run emt cu1.traj -E 9,3 -t eos --properties ef -o cu-fit.extxyz
ase eos cu1-eos.traj
-E 9,3指9个晶格长度点、±3%,不是体积±3%。默认EOS类型是 sj;如需其他方程,先看 ase eos --help,并在生成与提取时保持一致。
本例默认拟合得到平衡原子体积 11.565543ų、体积模量 134.315411GPa。与13点Birch–Murnaghan示例相比,采样范围和拟合方程不同,末位数值无需相同。不能用拟合程序成功作为参数收敛的证据。

5 读取真实Si能带和已有NEB#
先回到解压根目录执行过 bash run.sh,再进入输出目录:
cd work/output
ase band-structure si-bands.json -o my-bands.svg -r -6 6
ase reciprocal cu1.traj my-reciprocal.svg
ase nebplot --nimages 7 ../inputs/neb.traj my-neb.pdf
Si数据来自集群既有VASP计算,Γ–X–L共40个k点、32条带;附件从原始EIGENVAL转换成ASE BandStructure JSON,逐值保留本征值。能量参考使用OUTCAR中的E-fermi=5.8961eV;CLI的 -r -6 6表示相对该参照的显示窗口,图中原生坐标仍为对应的能量数值。共享X端点按两段各保留一次,图上可能显示X,X。

Cu NEB轨迹保存了优化过程中的多组图像,每组7帧,所以明确指定 --nimages 7;最后7帧是最终路径。不能把整个轨迹的帧数误认为NEB图像数。这里复用原有EMT路径,不重复进行NEB计算。
6 无显示窗口也能查看和导出#
ase gui -t -g i,e ../inputs/md.traj
ase gui -t -r 2,1,1 cu32.traj -o cu64.traj
ase ulm cu32.traj -n 0 -v
gui -t使用终端模式;第一个命令逐帧打印索引和能量,第二个把32原子胞沿x重复成64原子并导出。ulm查看ASE轨迹的底层存储信息。交互窗口需要本地显示或远程图形环境,本页终端步骤不要求显示服务。
7 其余入口怎样使用#
数据库、结构搜索、差异、逐帧执行和维度分析见ASE数据管理。还可以查询:
ase completion
ase test --list-calculators
completion只打印Bash补全配置行,不自动修改个人文件。test --list-calculators只查询计算器名称;运行完整测试套件需要pytest,本环境未安装,因此不能把这条查询解释成测试套件通过。
18个入口为help、info、test、gui、db、run、band-structure、build、dimensionality、eos、ulm、find、nebplot、convert、reciprocal、completion、diff、exec。附件保存了26条实际调用及其输出;本页不把入口演示等同于全部参数、GUI操作或外部计算器覆盖。
接口依据:ASE官方命令行说明,本地行为以 ase 子命令 --help 为准。