简介CASTEP-19.11.tar.gz是Materials Studio等软件中CASTEP模块的19.11版源码与文档压缩包面向从事固体材料电子结构计算的研究人员、材料科学/物理专业高年级学生。该版本基于密度泛函理论支持几何优化、能带结构、频率分析、电荷密度、反应路径等典型计算适合需要在Linux集群上自行编译部署CASTEP的用户。压缩包共4542个文件大小20.12MB内含233个C源码、222个Fortran源码、1317个param输入参数、1067个cell结构文件、99个cif晶体结构示例、95个usp赝势及大量makefile、configure编译脚本、tex说明文档和py后处理工具目录结构完整便于按模块查阅。目前已有1595人学习下载。解压后即可获得完整源代码、编译指南、文档与示例用户可依据指南安装依赖并编译或直接参考示例输入文件快速上手是学习CASTEP编译与使用、了解新版算法改进的实用资料。1. 你拿到的 CASTEP-19.11.tar.gz到底比旧版多了什么做计算材料的人对 CASTEP 应该不陌生但真正动手把 CASTEP-19.11.tar.gz 从压缩包变成能算任务的二进制是另一回事。这个包本身是 SourceForge 上 CASTEP 官方发布的 19.11 版本源码压缩包Linux 下最常见的分发形式里面是完整的 Fortran 源码、文档和构建脚本。19.11 这个版本号对应的是 2019 年 11 月的周期性维护版本跟更早的 18.1、17.2 相比它对 Intel 编译器、AMD 平台的兼容性更好还默认启用了 OpenMP MPI 混合并行也支持了 ENCUT 自动判断、TS 搜索的改进等一批实用功能。对大多数课题组来说这意味着同样的硬件条件下旧版本跑不动或收敛不了的任务在这个版本里有更高的概率直接跑通。这篇文章就是顺着这个 tar.gz 包展开的实战笔记。读者是刚拿到源码包的学生或者是组里负责维护计算服务器的工程师——你们关心的不是 CASTEP 是什么而是这个包怎么编译、参数怎么设、任务怎么往下算。后面章节会从解压检查开始依次走完依赖库准备、编译配置、性能调优、常见报错排查最后落到日常使用习惯上。如果你手里正好有个编译到一半就卡住的 19.11或者打算从旧版本迁移过来这篇文章应该能帮你省下几天的折腾时间。2. 解压之前先看包内结构不要急着敲 make2.1 解压命令与包结构检查我见过不少人拿到 CASTEP-19.11.tar.gz 就迫不及待地 tar -zxvf解压完直接 ./configure 或者 make然后被一堆缺失的头文件劝退。其实这个包的结构是固定的先花两分钟检查一下能避免后面一半的编译错误。最常见的解压方式是tar -zxvf CASTEP-19.11.tar.gz cd CASTEP-19.11 ls -la解压后你应该看到一个典型布局根目录下是 charset、Docs、Examples、Src还有 install 相关的脚本。打开 Docs 目录里的 ReleaseNotes.txt第一件事是确认你下载的版本号和日期避免拿到的文件被传输过程损坏。ReleaseNotes 里还会写明这个版本需要的编译器版本下限——19.11 对 Intel 编译器要求 2017 版本以上对 GCC 要求 6.0 以上这是个硬门槛低于这些版本编译时会有大量莫名的语法报错。包内校验也是值得养成的习惯。源码包在传输过程中偶尔会出问题很多人忽略这一步等编译到一半报错才开始怀疑是不是包坏了。我的做法是先在下载目录里检查压缩包的完整性gzip -t CASTEP-19.11.tar.gz md5sum CASTEP-19.11.tar.gzgzip -t 只校验压缩流完整性md5sum 用于记录一个常规基线。如果你的文件来源渠道本身提供了哈希值直接比对即可如果没提供至少确认 gzip -t 输出正常。这一步不需要花很多时间但能直接过滤掉一类最让人抓狂的“随机编译错误”。2.2 源码包内目录如何决定你的编译策略CASTEP-19.11 源码包的 Src 目录内部结构基本决定了你接下来要走哪条编译路径。打开 Src 目录会有 common、modules、io、parser 这些子目录共同点是全部代码在编译前由一个定制脚本预处理生成机器相关的 Fortran 源文件。这意味着 CPATH 或 Fortran 编译器版本不对你在预处理阶段就会被卡住而不是等到链接阶段。我一般会在解压后先执行一次./install.sh --check这是了解当前环境最直接的手段。如果包内没有这个脚本就去 Docs 目录看安装说明不同途径拿到的包可能在打包粒度上有差异但构建系统本质上是同一套。这里刻意强调“先检查”因为常见做法是直接跑./install.sh靠脚本的交互式配置去设置编译器和数学库路径而交互式脚本一旦碰到默认路径不对就会自动跳过去编译时才发现某些库没链接上。./install.sh --check这个命令会输出一个环境检查快照包括哪些编译器可见、哪些数学库在当前环境能找到。注意--check参数不是所有 CASTEP 分发版本都自带拿不准就去 Docs 里翻一下说明。它检查到缺失项时不会中断安装流程只会留下几行警告你要做的是记下这些警告在下一步编译命令里手动指定库位置。这一步是后面所有编译配置的基础别偷懒跳过。3. 编译 CASTEP-19.11 之前的三板斧编译器、数学库、环境变量3.1 先定编译器Intel 还是 GCCCASTEP-19.11 的源码对 Fortran 标准特性依赖较重编译器选错了后面全是泪。业界在这个版本上默认推荐 Intel 编译器如果你有正版或者学术许可选 ifort/ifx 会让你后面的日子好过很多。Intel 编译器对 CASTEP 的优化程度最好生成的二进制效率也更高。用 GCC 倒是也能编译但有两个明显的坑一是 gfortran 版本低了编译报错版本太高了链接阶段被某些新选项卡住二是同样数量的核GCC 编译出来的二进制跑起来会比 Intel 慢 10%~20% 不等。我的建议是先看看你的集群上装了什么。如果节点上已经有 Intel oneAPI 或者 Parallel Studio XE直接走 ifort MKL 路线。如果只有开源的 GCC 环境关注gfortran --version是否 6.0。需要说明的是系统的 GCC 是 4.8.x 这种老版本的话编译 CASTEP-19.11 基本等于在十字路口迷路该考虑的是怎么装新一点的编译器而不是在旧版本上浪费时间调参数。# 检查当前环境里可用的编译器 which ifort ifx gfortran ifort --version 2/dev/null || ifx --version 2/dev/null gfortran --version | head -n 1输出结果里如果有 ifort 或者 ifx那就优先用 Intel 路线。如果只有 gfortran继续往前走也要心里有数后面配置数学库的时候BLAS/LAPACK 要么用纯参考实现要么用 OpenBLAS折腾程度比 MKL 稍高。3.2 数学库是编译成败的分水岭MKL、FFTW 和 ScaLAPACKCASTEP 是平面波 DFT 代码它的大部分时间花在 FFT 和线性代数上所以数学库的选择直接决定你编译出来的东西到底能不能用。CASTEP-19.11 支持两种数学库组合一种是 Intel MKL 全家桶包含 BLAS、LAPACK、FFTW 接口和 ScaLAPACK另一种是 OpenBLAS/参考 BLAS FFTW 独立的 ScaLAPACK。对大多数课题组来说MKL 是最省心的方案因为接口一致性最好不用手工拼装。这里有个容易误解的地方CASTEP 源码包本身不捆绑数学库也不捆绑 MPI 库。你以为是“解压即可编译”实际是“解压只是开始”。编译前你必须确认以下东西都已经在系统里并且版本兼容# 假定用户用 Intel oneAPI source /opt/intel/oneapi/setvars.sh --force # 检查 MKL 是否可用 echo $MKLROOT ls ${MKLROOT}/lib/intel64/libmkl_scalapack_lp64.a如果$MKLROOT是空的说明 Intel 环境没有正确 source 进来或者根本没有装 oneAPI。这时候就得考虑 MPI 是不是也缺失了which mpifort mpirun同样需要检查。MPI 方面19.11 版本并不会强制使用某一种 MPI但 Intel MPI 配 Intel 编译器、OpenMPI 配 gfortran 是更不出错的组合交叉使用偶尔能跑通遇到问题排查起来却要double工作量。3.3 把环境变量一次性配齐后面少受罪编译 CASTEP-19.11 的时候环境变量没有配好最常见的现象是cs_make能开始跑但过一会儿就报cannot find -lfftw3或者libmpi_mpifh.so: cannot open shared object file。理论上数学库路径都写进cs_make的配置里就不会出问题但你不可能每次都手动改 Makefile 里的几十行变量所以你需要的是把这套环境固化到 shell 配置里。# 以 bash 为例把以下内容放到 ~/.bashrc 或 ~/.bash_profile export PATH/opt/intel/oneapi/compiler/latest/linux/bin/intel64:/opt/intel/oneapi/mpi/latest/bin:$PATH source /opt/intel/oneapi/setvars.sh --force /dev/null 21 # FFTW 接口由 MKL 提供无需再独立编译 FFTW export CASTEP_FFTW1 export CASTEP_ARCHlinux_x86_64_ifort注意CASTEP_ARCH这个变量不是官方制定的必选设置但它是很多编译脚本里用来识别平台和编译器组合的识别符如果版本匹配会自动选择对应的预编译配置。实际使用中还是以你手里包内提供的 ARCH 选择为准不要照抄这里。设置好后新开一个 shell 确认echo $MKLROOT有输出再进行下一步。我再强调一遍这一步的目的是把“已知好用的环境”固定下来不是为了装样子。如果你跟别人共用一台服务器这些变量应该写进你自己的工作目录下的环境脚本里每次编译前source env.sh即可别去动系统级的 profile否则会给其他用户带去困扰。4. 从 tar.gz 到 castep.mpi编译配置与两个关键命令4.1 使用命令行配置的两种方式交互脚本与非交互参数CASTEP-19.11 的构建方式有种“半自动”的感觉它不像 cmake 那样有个统一的 CMakeLists.txt 入口而是提供了一套基于 make 的构建流程。实际项目里我通常用./install.sh交互式配置但交互脚本有时候在后台服务器上跑会卡在等待输入的地方或者某些终端环境下会把选项显示得乱七八糟所以更推荐的方式是把参数一次性写在命令行里即使某个参数不生效也比卡住要好排查。./install.sh --compiler ifort --mpi intelmpi --math mkl --arch linux_x86_64_ifort --prefix /your/local/install--compiler指定 Fortran 编译器--mpi指定 MPI 实现--math指定数学库--prefix指定最终安装路径。这些参数的命名和可选值会因分发版本略有差异但思路一致。如果你手里的包没有install.sh或者参数对不上就去 Docs 里的构建文档里查对应的 make 目标名通常在Src目录下直接执行make target也能完成编译只是没有自动安装步骤。配置过程本身并不保证万事大吉。参数识别是一回事实际的编译错误又是另一回事。我见过最多的情况是--mpi intelmpi已经指定了但 make 阶段依然用了 mpif90 去连接 MPI 程序导致一堆 MPI 符号找不到。这时候你需要进入Src目录找到生成的 Makefile检查里面MPI_F90或者F90变量到底指向哪个路径必要时手动修正后再编译。4.2 编译命令与并行选项cs_make 是入口编译命令的核心是cs_make。这个名字在 CASTEP 语境里等同于“Castep make”你用make直接编译也不是完全不行但cs_make会额外处理一些依赖顺序和平台相关的配置所以推荐入口是cs_make。cd Src ./cs_make -j8-j8表示同时用 8 个进程编译具体数值按机器 CPU 核数来定一般取物理核心数或者稍低。注意这里不要用make -j的惯用思维直接把数字拉满因为 CASTEP 编译过程中预处理和编译之间的文件依赖很长并行进程太多会把 IO 拖爆反而更慢。你可以边编译边看输出里有没有 with warnings 字样这个阶段最常见的错误是把数学库路径设错表现为cannot find -lfftw3这类信息。如果cs_make中途失败不要立刻重新运行一遍先找到最后一屏的报错行。根据我维护计算服务器的经验90% 的情况下是配置阶段埋下的问题——某个库路径不对或者某个编译选项跟编译器版本不匹配。清掉之前编译产生的所有中间文件再重试是标准操作./cs_make distclean ./cs_make -j4distclean会删除所有之前生成的 .o 和 .mod 文件给你一个干净的编译现场。重新编译时建议降一下并行数便于观察错误输出来自哪一步。4.3 二进制产物验证拿到 castep.mpi 后先做静态确认编译成功后的目录下会出现几个可执行文件最核心的是castep.mpi另外还有castep.serial在某些配置下也会生成。先用file命令确认二进制格式正确再执行一个最简单的--help或者无参数运行看它能不能打印版本信息file bin/castep.mpi bin/castep.mpi --help 21 | head -n 20如果输出的是 ELF 格式信息并且能看到 CASTEP 版本信息或者命令行参数列表说明编译阶段已经通过了。这时候先不要高兴得太早建议编一个包含三个原子的简单体系跑一次单点能计算验证结果是否合理。CASTEP 的 Docs 目录里自带了一些 Examples你可以选Examples/xi这类最基础的工作目录直接复制过来运行。cd Examples/xi cp ../InputFiles/xi.param . cp ../InputFiles/xi.cell . mpirun -np 2 ../../bin/castep.mpi xi如果能顺利跑完并且.castep文件里出现收敛信息才可以认为这一步编译成功。第一次跑通后把这个命令和运行目录记录下来后面调优和部署都用得上。5. 性能调优与边界从跑通到省时的关键四步5.1 先看 castep.mpi 的并行效率再说调参编译只是第一步跑得快才是目的。19.11 版本的 CASTEP 支持混合并行但默认情况下计算效率未必最优。当你在单机或集群上编译好之后第一次实际运行任务时重点观察castep输出文件里每一步计算时间的分布。大部分时间应该花在 FFT 和矩阵对角化上如果某个环节的时间占比异常高那通常说明某类参数需要调整。常见做法是先用一个中小体系做并行测试。选取一个 10~30 原子的结构分别用 1、2、4、8、16 个核跑一遍单点能看加速比曲线。CASTEP 的并行是跨 k 点和 band 的体系越大并行度越好体系太小核数开多了通信开销反而拖垮整体时间。我刚上手时犯过一个错拿一个 8 原子的分子用 64 核跑跑完比 8 核还慢一倍纯粹是通信成本超过了计算收益。# 典型并行测试命令 for np in 1 2 4 8 16; do mpirun -np $np ../../bin/castep.mpi Si16 grep -E Total CPU time|Total wall time Si16.castep done对比不同核数下的 wall time如果从 8 核到 16 核加速比低于 1.5 倍基本可以判断这个体系不适合开到 16 核或者你在编译时没有启用 OpenMP 层的混合并行支持。5.2 设置合理的 ENCUT 与 k 点不要照抄默认值CASTEP-19.11 有几个参数经常造成翻车最常见的当属cut_off_energy。默认值有时偏小结构优化或者计算能量学量时精度不够。对于一般金属体系400 eV 起步是底线氧化物或半导体体系通常要 500 eV 以上。收敛测试是判断cut_off_energy是否合适的最直接手段但没有必要每个体系都从 200 eV 开始扫合理的做法是先在 300、400、500、600 eV 各算一次单点能看总能量变化是否已经进入你想要的误差范围内。k 点密度同理。金属体系和小带隙半导体需要的 k 点密度通常很高绝缘体则相对稀疏。CASTEP 里通过kpoint_mp_grid设置# xxx.cell 文件里设置 kpoint_mp_grid 8 8 8 kpoint_mp_offset 0.05 0.05 0.05kpoint_mp_grid是 3 个整数表示在倒空间三个方向上生成多少分点。kpoint_mp_offset是 Monkhorst-Pack 网格的偏移量大多数体系用默认的 0 0 0 就行金属体系有时会需要加一点偏移来改善 Fermi 面附近的分辨率。不要小看这两行配置很多论文里的结构优化结果不收敛问题就出在这里。5.3 合理控制空带数别让 Dipp 算法滑向边缘CASTEP-19.11 中空带数由num_dump_cycles和nextra_bands参数间接控制。默认情况下 CASTEP 会根据体系自动分配一部分空带但如果你计算的是宽带隙半导体、绝缘体或者含有过渡金属/稀土元素的体系默认空带数常常不够。空带数不足的后果是 SCF 循环极难收敛或总能量出现非物理的缓慢漂移。一个实用的尝试是逐步增加nextra_bands比如从默认值增加到 10% 或 20% 的能带数观察 SCF 收敛代数和最终能量变化。不要一次性给到 50% 以上会白白拉长计算时间。输出电压致密度、介电函数这类需要较高空带数的性质时nextra_bands的作用就很关键可以直接决定结果是否可信。5.4 固定好一批运行脚本避免每次手工拼命令集群环境中跑 CASTEP 通常要提交作业脚本不要每次手动敲 mpirun。一个比较通用的作业脚本模板如下#!/bin/bash #$ -N ca_opt #$ -pe mpich 16 #$ -cwd source /your/env/path/env.sh mpirun -np $NSLOTS ../../bin/castep.mpi Si16脚本的作用是固定住环境变量、并行配置和输入文件路径。许多课题组里的初学者在相同节点上跑了很多次都失败原因就是这次忘了 source 环境变量那次 mpirun 参数写错做法上是白费了许多时间。我把这套模板固化下来后组里的学生基本上开箱即用。6. 新手最容易掉的坑19.11 编译运行中的排查清单6.1 现象编译报错compilation aborted for file xxxx.f90原因编译器版本过低或预处理阶段宏定义冲突。CASTEP-19.11 在代码里使用了较新的 Fortran 特性GCC 6.0 以下基本过不了预处理。解决首先确认编译器版本如果是 gfortran建议换一台有新版编译器的节点或者用 Intel 编译器如果版本没问题把报错文件路径和行号贴出来搜一下大概率是某个宏需要关闭例如-DCOMMS_MPI之类的定义冲突。6.2 现象mpirun无法启动或一启动就崩原因MPI 库与编译器匹配不上。比如用 ifort 编译却用 OpenMPI 去跑或者集群上不同节点 MPI 版本不一致。解决严格按编译时使用的 MPI 启动方式运行。Intel MPI 用mpirun或mpiexec.hydraOpenMPI 用的命令是mpirun但是两者在进程管理上机制不同。最直接的检查办法是ldd castep.mpi | grep mpi看清楚它到底链接的是哪个 MPI 库再决定用什么命令拉起任务。6.3 现象SCF 循环一直收敛不了原因这个太多重可能。最常见的是cut_off_energy太低或者空带数不足特别是磁性体系没有设置spin_polarized。解决先提高cut_off_energy往 500 eV 附近靠再把nextra_bands加到 10% 试一遍如果体系带有强关联电子检查U设置是否生效。不要先急着去调mix_scheme这些电子步参数那是在浪费调试时间。6.4 现象几何优化前几步很顺利后几步能量震荡、跑飞原因原子坐标或者晶胞设置有问题常见的是多个原子位置重叠也有可能是最大步长设置太大结构往不合理方向闯。解决先检查原始结构用可视化软件重新观察原子间距然后看看.castep文件里的受力是否在合理范围。实在找不到原因就把geom_max_step调小比如设成 0.01 Å然后再慢慢恢复。6.5 现象运行中途NaN或Inf出现原因不同体系的数值稳定性问题通常是初始磁性状态或者波函数初猜不当也有时候是 k 点网格设置与体系对称性不匹配。解决修改输入文件加入初始磁性设置spin 3*2.0这类初始磁矩如果是纯非磁体系确保结构文件中没有多余的对称性重复原子必要时降低kpoint_mp_grid到较低的网格先跑通再逐渐加密。7. 一台服务器上的日常技巧CASTEP-19.11 的隐性设置与备用手段编译完成、算例能跑这只是迈过了门槛。日常维护中还有几个细节值得养成习惯。第一个习惯是不要往系统路径里装 CASTEP而是每个项目维护一个独立的运行目录里面放一个env.sh和一组作业脚本。这样的好处是别人在你这台机器上跑任务或者你升级了编译器环境不会因为路径冲突导致重建整个编译链。具体做法是把castep.mpi放在固定的安装目录然后在每个项目目录下用一个符号链接指过去ln -s /your/install/bin/castep.mpi /project/jobs/castep.mpi第二个是学会利用--dryrun或者并行检查类命令去审查输入文件而不是浪费节点资源去试错。某些情况下一个明显错写的cell文件会直接在 SCF 第一步就爆炸但如果你先用快速检查模式让它报出输入文件错误就能省下大量排队时间。这个检查模式在不同分发版本中的名字不完全一样有些版本里是castep.serial seed --check有些是castep.mpi seed --check具体以你编译出的二进制--help输出为准。第三个习惯跟效率调优有关CASTEP-19.11 带的O(N)算法在大型非周期性体系里很有价值但默认运行并不一定开启。如果你的体系是上千原子的非晶结构或者大表面模型且只关心基态能量和结构弛豫可以尝试在cell文件里加入O(N)相关关键字。它的收敛性比常规方法敏感但一旦收敛计算成本明显低于平面波对角化方案。这部分是加分项新手还是先把标准流程跑通。还有一个小技巧是对计算结果做重复利用——把算好的波函数文件.castep_bin保留下来后续需要接着算态密度、能带时在参数文件里加上# xxx.param 文件里 continuation : default只要顺序变化不大continuation可以沿续旧的波函数初猜跳过重新 SCF 的成本。这招在长流程的几何优化和过渡态搜索里非常省时但要注意体系也发生了剧烈变化时沿用旧波函数反而会拖慢收敛这时删掉.castep_bin重新从头算反而是最快的。回到开始时的那个关键词CASTEP-19.11.tar.gz不只是个文件名它是你进入一套材料模拟工作流的入场券。编译它、跑通它、理解它每个环节背后都对应着不同的计算思维习惯。希望这篇笔记能帮你在自己的服务器上少踩几个坑早日跑出像样的数据。我自己的习惯是每次新拿到一个版本先在一个独立目录里编译再把默认配置记录下来之后所有项目都从这份记录起步——这个习惯帮我省了不少时间希望能帮到你。本文还有配套的精品资源点击获取