VASP和Linux这两个词凑在一起，基本就是国内做第一性原理计算新人的日常开局。你大概率是材料、物理或化学方向的研究生，导师丢过来一个超算账号，说了句"跑个结构优化看看"，然后你就对着黑色终端发呆——ls、cd、vim、INCAR、KPOINTS，每个单词都认识，拼在一起就不知道从哪下手。这种情况太常见了，我当年也是从"连文件都传不上去"开始熬过来的。
这篇内容想做的事情很具体：把VASP和Linux这两块最劝退的门槛，用一条完整的实操链路串起来。前半段讲Linux到底要会多少就够用，哪些命令是必须肌肉记忆的，哪些可以现查现用；中段把VASP四大输入文件逐个拆开，讲清楚每个标签背后在控制什么物理量；后半段给一个能直接抄的MgO结构优化实例，从建POSCAR到判断收敛，一步不落。适合完全没碰过Linux命令行的新手，也适合跑过几个算例但一直不太清楚参数含义的半熟手。我不会给你堆一堆手册原文，而是把"为什么这么设""设错了会怎样""报错怎么救"这三件事讲透。
1. 为什么VASP新手的第一课是Linux而不是VASP本身
1.1 计算软件的运行逻辑决定了你的学习顺序
VASP不是那种双击图标就能弹窗的软件。它是一套跑在Linux环境下的并行计算程序，输入是几个纯文本文件，输出也是一堆纯文本文件，中间没有任何图形界面帮你兜底。你所有的操作——建目录、改参数、提交任务、看日志、判断收敛——都要靠命令行完成。这就意味着Linux不是"顺带学一下"的附属技能，而是你操作VASP的唯一通道。
很多人卡住不是因为不懂密度泛函理论，而是因为不会用grep在几千行的OUTCAR里找那一行关键信息，或者不会用cp备份上一版INCAR，改崩了只能重来。这种卡顿极其消耗耐心，一天下来算例没跑成，人先烦了。我自己带过几个低年级同学，前两周专门花时间练Linux基本操作，后面的效率比直接上手VASP的人高出不止一倍。所以顺序很清楚：先让命令行不吓人，再让VASP不神秘。
1.2 三种常见上手环境，选哪个更省事
新手面对的第一个现实问题是：到底在哪里跑。常见有三条路，各有取舍。
本地虚拟机。在自己笔记本上装个VMware或VirtualBox，再装Ubuntu或CentOS。优点是随便折腾，装崩了重装就行，不用怕搞坏生产环境；缺点是笔记本性能有限，只能跑两三个原子的体系练手，稍微大一点就卡死。
课题组或学院的服务器。通过SSH登录，多核CPU和内存都够用，是绝大多数人真正干活的平台。你需要熟悉ssh登录、scp传文件、nohup后台跑任务这几件事。缺点是多用户共享，资源要抢，管理上有约定俗成的规矩。
超算中心。资源最强，但一般用SLURM或PBS作业调度系统，提交任务要写作业脚本，入门门槛最高。等你把服务器用熟了再迁过去，节奏比较顺。
我的建议是虚拟机练手一周，熟悉命令和文件操作，然后立刻转到服务器上做真实算例。虚拟机跑不出有意义的科学结果，它的价值只是让你练习手速。
顺带一提，国产Linux发行版这两年在科研服务器上的占比在上升，openEuler、统信UOS这类系统底层还是标准的Linux内核和shell，你学的那套命令完全通用，不用担心换个发行版就白学。
1.3 新手头一周最容易卡住的三个动作
按我的观察，新人报废时间最多的不是算例本身，而是这三个动作：
第一，文件路径搞混。在~/calc/MgO目录里建了文件，人却在~/calc目录下执行VASP，程序找不到POTCAR直接报错退出。养成用pwd确认当前位置的习惯，能省掉一半的困惑。
第二，文本编辑器逃不出去。用vim打开INCAR，改了内容，结果不知道怎么保存，按了一堆键，文件变成一堆乱码。这是最难堪也最普遍的一幕。解决办法很简单，先只记四个操作：i进入编辑、Esc退出编辑、:wq保存退出、:q!不保存强退。够用了。
第三，不知道任务在不在跑。提交之后终端没反应，以为死机了，其实是正常在算。这时候要用ps -ef | grep vasp确认进程还在，用tail -f OSZICAR实时看能量往下掉。
这三个坎迈过去，后面就是纯粹的VASP知识问题了。
2. 够用就好的Linux命令集，附速查表
2.1 导航与文件操作：每天都要用的十个命令
Linux命令有几百个，但跑VASP真正高频的也就十来个。别一上来就抱着"常用100个命令"背，记不住也用不上，用到再查才有效。下面这张表是我自己常用的清单，建议先熟到不用看表。
| 命令 | 作用 | 典型用法 |
|---|---|---|
pwd | 显示当前所在目录 | pwd |
ls -lh | 列出文件，带大小和权限 | ls -lh *.vasp |
cd | 切换目录 | cd ~/calc/MgO |
mkdir -p | 递归创建目录 | mkdir -p bulk/relax |
cp | 复制文件或目录 | cp INCAR INCAR.bak |
mv | 移动或重命名 | mv POSCAR POSCAR.v1 |
rm -rf | 强制递归删除 | rm -rf tmpdir |
cat | 输出文件内容或拼接 | cat POTCAR_Mg POTCAR_O > POTCAR |
du -sh | 查看目录占用空间 | du -sh ~/calc |
df -h | 查看磁盘剩余空间 | df -h |
这里有两个坑必须提前说。rm -rf删掉的东西不进回收站，rm -rf ~/calc手滑一下，一个月的计算数据就没了。我的习惯是在删之前先把命令敲进终端，停三秒，确认路径拼写没错再回车。另一个是cat拼接POTCAR时的重定向符号，>是覆盖，>>是追加，第一次拼POTCAR如果误用>>，会把旧内容全留下，VASP读进去直接元素顺序对不上。
提示：所有涉及删除的命令，先用
ls看一遍目标路径里到底是什么，再执行删除。这个习惯能救你无数次。
2.2 查看与编辑：less加vim的最小可用组合
看OUTCAR、OSZICAR这类文件，别用cat直接刷屏，几千行在你眼前飞过去，什么都没看清。用less OUTCAR分页看，方向键上下翻，/关键词搜索，按n跳到下一个匹配，q退出。想看文件末尾的收敛信息，tail -n 50 OSZICAR比翻页快得多。想在OUTCAR里定位那行最终的收敛提示，grep "reached required accuracy" OUTCAR一行就出来。
编辑文件用vim，新手先掌握这些：
vim INCAR # 打开文件 # 按 i 进入插入模式，开始打字 # 改完按 Esc 回到普通模式 # :wq 回车 保存并退出 # :q! 回车 不保存强退再进阶一点的操作：普通模式下dd删一整行，yy复制一整行，p粘贴，u撤销，/ENCUT搜索，:%s/400/520/g全局把400替换成520，:set nu显示行号，:80跳到第80行。这几招组合起来，改INCAR、修POSCAR完全够用。有人喜欢用nano，操作更直观，但服务器上vim几乎一定装了，nano不一定，还是练vim更保险。
2.3 权限、进程与后台运行：让任务自己跑
VASP算一个中等体系动辄几个小时，你不可能守着终端。这时候nohup就派上用场了：
nohup mpirun -np 24 vasp_std > vasp.log 2>&1 &末尾的&让命令转入后台，nohup保证你断开SSH后进程不被打断，> vasp.log 2>&1把标准输出和报错都写进日志文件。提交完可以用ps -ef | grep vasp确认进程在，或者top看CPU占用。要停掉任务，先查到进程号，再kill -9 进程号。
权限方面，新手最容易遇到的是脚本无法执行。./run.sh报"权限不够"，chmod +x run.sh加上执行权限就好。文件权限那串rwxr-xr-x的含义可以用一句话记住：三段分别对应文件所有者、同组用户、其他人，r读w写x执行。
2.4 环境变量与软件加载：解决"找不到vasp命令"
装好VASP后运行vasp_std提示"command not found"，说明程序路径没进环境变量。两种做法：
export PATH=$PATH:/opt/vasp/bin这条命令只在当前终端生效，关掉就没了。想永久生效，写进~/.bashrc，然后source ~/.bashrc重新加载。这是最通用的方式。
超算平台更常见的是用module管理：
module avail # 看有哪些软件 module load vasp/6.4.2 # 加载VASP module load openmpi/4.1 # 加载并行库 module list # 确认已加载注意VASP的并行版本必须和加载的MPI库匹配，用Intel MPI编译的VASP配OpenMPI的mpirun，十有八九崩在启动阶段。这个坑后面排查章节还会细说。
3. VASP四大输入文件到底在控制什么
3.1 POSCAR：结构怎么被程序读懂
POSCAR描述晶胞和原子位置，格式看着简单，错一个数字整个计算就废了。它分几块：第一行注释（随便写，但建议写清楚体系名称），第二行是全局缩放因子，接着三行是三个晶格矢量，然后一行元素符号，一行对应元素数量，再一行决定坐标是Direct（分数坐标）还是Cartesian（笛卡尔坐标），最后是各原子的坐标。
一个MgO岩盐结构的例子：
MgO bulk 1.0 4.2100000000 0.0000000000 0.0000000000 0.0000000000 4.2100000000 0.0000000000 0.0000000000 0.0000000000 4.2100000000 Mg O 1 1 Direct 0.0000000000 0.0000000000 0.0000000000 0.5000000000 0.5000000000 0.5000000000几个容易翻车的地方。元素数量那一行必须和你实际坐标行数对得上，写1 1就得有2行坐标。坐标单位要看清是Direct还是Cartesian，Direct是相对于晶格矢量的分数，Cartesian是绝对长度（Å），两者换算靠晶格矩阵，混用会导致结构完全错乱。还有，POSCAR里的元素顺序就是后面POTCAR必须遵循的顺序，这一点极其关键，下一节细说。
如果你习惯用Materials Studio建模型，导出成cif或xsd后可以用VESTA、ASE或pymatgen转成POSCAR。转完一定要用VESTA打开看一眼，确认原子没跑到晶胞外面、没有重叠。我见过有人转完忘了去掉周期性镜像原子，晶胞里多了好几个重复原子，算出来的能量完全没意义。
3.2 POTCAR：赝势拼装与元素顺序陷阱
POTCAR是赝势文件，描述每个元素的芯电子近似。它不是你自己写的，是从赝势库里复制拼接出来的。以MgO为例：
cat POTCAR_Mg POTCAR_O > POTCAR问题来了：这里Mg在前O在后，POSCAR里元素那行也必须是Mg O，顺序完全一致。如果POSCAR写O Mg而POTCAR拼成Mg O，VASP会在启动时直接报错：
ERROR: the type information in POTCAR is not consistent with POSCAR或者更隐晦一点，不报错但结果全错——因为程序按顺序把Mg的赝势套到了O的坐标上，算出来的东西毫无物理意义，而且你很难第一时间发现。
另一个关键是赝势必须来自同一套库。VASP常见的有potpaw_PBE、potpaw_PBE.54、potpaw_LDA等。你不能Mg用PBE的、O用LDA的，混着用会导致元素间描述精度不一致，能量不可比。目前主流是PBE泛函，配合.54版本的赝势。选哪个元素用哪个版本（比如Mg用Mg还是Mg_pv），取决于你要不要把那几个半芯态电子显式当价电子处理，_pv、_sv表示把更多半芯电子纳入价电子，精度更高但计算量更大。新手做常规结构优化，用默认版本就行。
拼完POTCAR，建议grep ENMAX POTCAR看一眼每个元素的截断能建议值，下一步设ENCUT要用到。也建议grep ZVAL POTCAR看价电子数，估算体系总电子数，对判断KPOINTS和内存有用。
3.3 INCAR：常用标签逐条拆开讲
INCAR是控制计算类型的核心文件，标签有上百个，但结构优化常用的就十几个。一个个过。
SYSTEM只是个注释，任意字符串都可以，方便你区分不同算例，我一般写成体系名加任务类型。
ENCUT是平面波截断能，单位eV。它决定了基组的大小，设得越高越精确也越慢。经验规则是取所有元素中最大ENMAX的1.3倍。Mg和O的PBE赝势ENMAX通常是400 eV左右，取520 eV比较稳妥。低于建议值结果不可靠，高太多纯属浪费机时。
PREC控制精度档次，Normal是默认，够用；Accurate更严但更慢，对能量差要求高的场合（比如做相稳定性比较）才用。
EDIFF是电子自洽收敛判据，能量变化小于这个值就认为电子步收敛。常规设1E-5，高精度设1E-6。设太松结果粗糙，设太紧收敛慢甚至难收敛。
EDIFFG是离子步收敛判据。取正值表示按能量判据，取负值表示按受力判据，单位eV/Å。结构优化普遍用负值，-0.01意味着所有原子受力都小于0.01 eV/Å才算收敛，这是比较严格的标准；-0.02或-0.03更宽松更快。
ISMEAR和SIGMA控制占据数的展宽方式，这个参数经常被新手忽略但影响很大。金属体系用ISMEAR=1（Methfessel-Paxton一阶）或2，SIGMA取0.1到0.2；半导体和绝缘体用ISMEAR=0（高斯展宽），SIGMA取0.05；做态密度计算时为了让DOS平滑，反而要用ISMEAR=-5（四面体方法）。MgO是绝缘体，用ISMEAR=0、SIGMA=0.05。
NSW是最大离子步数，给个上限防止死循环，结构优化设100就够。
IBRION决定离子怎么移动，2是共轭梯度法，最常用最稳；1是准牛顿法，靠近极小值时收敛快；-1是静态计算不动离子。
ISIF控制优化哪些自由度，这是结构优化的关键。2只优化原子位置，晶胞固定；3同时优化原子位置、晶胞形状和体积；4优化形状和体积但原子位置不动；7只优化体积。做晶格常数优化用ISIF=3。
LREAL决定投影算符是实空间还是倒空间处理，小体系用.FALSE.（倒空间），大体系（超过20个原子）用Auto或.TRUE.加快速度。
LWAVE和LCHARG控制是否输出WAVECAR和CHGCAR，这两个文件很大。结构优化阶段不需要，设.FALSE.省空间；做后续静态计算或态密度时要打开，用来续算。
NELM是最大电子步数，默认60，收敛困难的体系可以加到100甚至更多。
把这些拼起来，一个MgO结构优化的INCAR长这样：
SYSTEM = MgO bulk relax ISTART = 0 ICHARG = 2 ENCUT = 520 PREC = Normal ALGO = Normal EDIFF = 1E-6 EDIFFG = -0.01 ISMEAR = 0 SIGMA = 0.05 NSW = 100 IBRION = 2 ISIF = 3 LREAL = Auto LWAVE = .FALSE. LCHARG = .FALSE. NELM = 80ISTART=0表示从头算（有WAVECAR时会读取）；ICHARG=2表示从原子电荷密度叠加开始，是最干净的起点；ALGO=Normal是标准对角化算法，稳定优先。
3.4 KPOINTS：k点密度怎么定才不吃亏
KPOINTS控制布里渊区积分采样密度。密度太低结果不准，太高浪费机时。自动生成的格式最省心：
K-Points 0 Gamma 8 8 8 0 0 0第一行注释随便写；第二行写0表示用自动网格；第三行Gamma表示Gamma中心网格（也可以写Monkhorst-Pack，两者对六方和三角晶系略有差异，立方晶系差别不大）；第四行是三个方向的网格数；第五行是偏移量，一般全零。
网格怎么选？经验做法是按倒空间间距算。对MgO这种立方晶系，晶格常数4.21 Å，用8×8×8基本够收敛。粗略的判断标准是每个倒格矢方向至少要有30到40个k点采样，能量对k点加密到1.5倍时变化小于1 meV/atom就算收敛。结构优化时k点可以略稀一点（省时间），最后一步做静态能量时加密算一次。
VASP 6有个方便的标签KSPACING，写进INCAR就不用KPOINTS文件了，程序自动按间距生成网格。但新手期我还是建议老老实实写KPOINTS，直观看到自己在设什么。
4. 一次完整的MgO结构优化实操
4.1 目录规划与文件准备
先把工作目录搭起来，养成"一个算例一个目录"的习惯：
mkdir -p ~/calc/MgO/relax cd ~/calc/MgO/relax然后把POSCAR、INCAR、KPOINTS、POTCAR四个文件放进来。POSCAR就是3.1节那个立方MgO。POTCAR拼接时注意顺序和POSCAR一致：
cat /path/to/potpaw_PBE.54/Mg/POTCAR \ /path/to/potpaw_PBE.54/O/POTCAR > POTCAR grep ENMAX POTCAR grep ZVAL POTCARgrep ENMAX确认截断能建议值，grep ZVAL看价电子数——Mg的ZVAL一般是2，O是6，总共8个价电子，方便你后面判断计算规模。
准备就绪后，四个文件ls一下确认都在：
ls # INCAR KPOINTS POSCAR POTCAR注意：POTCAR路径里的元素目录名必须和你POSCAR里写的元素符号完全匹配，大小写敏感，
Mg不能写成MG。
4.2 分两步走：先粗后精的收敛策略
一上来就ISIF=3、EDIFFG=-0.01、ENCUT=520地猛算，风险是如果初始结构不合理，第一步离子步就可能跑飞，几十个核时打水漂。我自己的做法是分两步：
第一步，粗优化确认结构没大问题。用ISIF=2只动原子位置，EDIFFG=-0.05放宽受力判据，NSW=30，快速跑完看CONTCAR里的原子有没有跑到很怪的位置，能量是不是合理。这一步通常十几分钟到半小时。
第二步，正式优化。把INCAR换成ISIF=3、EDIFFG=-0.01、NSW=100，从粗优化的CONTCAR接着算。这一步优化晶格常数和原子位置一起，收敛更彻底。
cp CONTCAR POSCAR # 把粗优化的结果作为下一步起点这个cp CONTCAR POSCAR是极高频操作，记住它。CONTCAR就是优化后的结构，格式和POSCAR一样，直接覆盖就能续算。
4.3 提交任务与并行设置
单机或服务器上直接用mpirun：
nohup mpirun -np 24 vasp_std > vasp.log 2>&1 &-np 24是并行核数，一般设成机器可用核数。核数不是越多越好，超过物理核心数反而变慢，而且VASP的并行效率随核数增加而下降，24核跑MgO原胞已经绰绰有余。体系大、k点多的时候，可以考虑-np 32配合k点并行（在INCAR里加KPAR=4把k点分成4组并行）。
如果用的是SLURM调度的超算，写成作业脚本：
#!/bin/bash #SBATCH --job-name=MgO_relax #SBATCH --nodes=1 #SBATCH --ntasks=32 #SBATCH --time=24:00:00 #SBATCH --partition=normal module purge module load openmpi/4.1.4 module load vasp/6.4.2 cd $SLURM_SUBMIT_DIR srun -n 32 vasp_std > vasp.log 2>&1提交用sbatch run.sh，查状态squeue -u 你的用户名，取消scancel 任务号。
4.4 输出文件解读与收敛判定
任务跑起来后，最常看的是OSZICAR和OUTCAR。OSZICAR实时更新，每行开头是DAV或RMM表示电子步，出现F=的那行是离子步结束的能量记录：
DAV: 1 -0.12345678E+02 ... DAV: 2 -0.23456789E+02 ... ... 1 F= -.34567890E+02 E0= -.34567890E+02 d E =-.123456E-03F=是自由能，E0=是外推到零展宽的能量，d E=是这步离子移动带来的能量变化。看着能量一路往下掉，d E越来越小，说明在收敛。当OUTCAR里出现这一行：
reached required accuracy - stopping structural energy minimisation就说明受力判据满足了，结构优化收敛。用命令直接抓：
grep "reached required accuracy" OUTCAR最终结果在CONTCAR里，最终能量和受力在OUTCAR里：
grep "energy without entropy" OUTCAR | tail -1 grep "external pressure" OUTCAR | tail -1ISIF=3优化完成后，external pressure应该接近0（几个kBar以内都算合理，十几kBar说明晶格还没完全松弛）。把CONTCAR和原始POSCAR用VESTA对比一下，看晶格常数变化了多少，MgO的实验晶格常数是4.21 Å，PBE泛函通常会略微高估，算出来4.24 Å左右属于正常范围。
5. 报错速查与踩坑经验实录
5.1 高频报错对照表
VASP报错千奇百怪，但新手遇到的翻来覆去就那几个。整理成表，遇到直接查：
| 报错信息 | 常见原因 | 处理办法 |
|---|---|---|
POSCAR and POTCAR are incompatible | 元素顺序不一致 | 检查POSCAR元素行和POTCAR拼接顺序 |
command not found: vasp_std | 环境变量或module未加载 | module load或检查PATH |
ZBRENT: fatal error in bracketing | KPOINTS太少或结构严重不合理 | 加密k点，检查POSCAR有无重叠原子 |
LAPACK: Routine ZPOTRF failed | 并行或内存问题 | 减小核数，改ALGO=Normal |
Error EDDDAV: Call to ZHEGV failed | k点或算法问题 | 降低k点密度，或设ALGO=Fast试探 |
internal error in subroutine SGRCON | 对称性判断失败 | INCAR加ISYM=0关闭对称性 |
your PSP version is not compatible | 赝势版本与VASP不匹配 | 换用对应版本的赝势库 |
Segmentation fault | 内存不足或MPI不匹配 | 减少核数，检查MPI与VASP编译版本 |
DENTET: can't reach specified precision | SIGMA设置不当 | 调整ISMEAR/SIGMA，金属用1或2 |
ISYM=0这个技巧特别值得记住。有些低对称结构、缺陷结构或表面模型，VASP的对称性检测会出错，加一行ISYM=0关掉对称性就能绕过，代价是计算量增加一些，但总比跑不起来强。
5.2 几条用时间换来的心得
关于k点和ENCUT的收敛测试。很多人图省事直接套别人的参数，结果体系不一样，算出来的能量差了个几十meV，还以为是物理效应。正确做法是每个新体系都做一次收敛测试：固定ENCUT，把k点从4×4×4加密到6×6×6、8×8×8，看能量变化；固定k点，把ENCUT从400加到500、600，看能量变化。当再加密能量变化小于1 meV/atom，就认为收敛了。这个测试花不了多少时间，但能让你后面所有结果都站得住脚。
关于SCF不收敛。电子步跑到NELM上限还没收敛，OSZICAR里能量来回震荡，这是最让人头疼的情况。常见原因和对应处理：SIGMA设太小（尤其金属体系），适当加大到0.1到0.2；ALGO用Normal还不够就试All或Fast；NELM加到200；实在不行用AMIX=0.2和BMIX=0.0001降低混合参数，让电荷密度更新得慢一点。这些都是在反复试错中攒下来的。
关于WAVECAR续算。结构优化跑完，如果要做静态能量或态密度，不用从头再算一遍。把LWAVE打开跑完优化，得到WAVECAR，然后复制到新目录做静态计算，INCAR里设ISTART=1、ICHARG=11，直接读取已有波函数和电荷密度，能省掉大量收敛时间。这个操作在大体系里省下的机时非常可观。
关于文件备份。每次改INCAR之前，cp INCAR INCAR.bak。听起来很啰嗦，但当你改了三个参数结果算崩了，想回到上一版却发现原内容忘光了的时候，就知道这个习惯多值钱。我现在的做法是每个算例目录下保留INCAR.v1、INCAR.v2这种带版本号的副本，配合一个简单的note.txt记录每次改了什么、为什么改。
关于中文文件名乱码。在Linux上解压Windows打包的压缩文件，中文名经常变成乱码，原因是压缩时用的字符集和你的locale不一致。可以用unzip -O GBK xxx.zip指定原始编码，或者干脆用7z x xxx.zip，兼容性更好。但最省事的办法是——建算例目录时全用英文和数字，MgO_relax而不是氧化镁优化，从源头避免这类问题。
关于别在共享服务器上乱来。多用户服务器上，kill -9别人的进程、占满磁盘、把任务提交到登录节点，这三件事都是禁忌。跑大规模任务前先df -h看一眼磁盘剩余，free -h看内存余量，确认资源够再提交。这些不是技术问题，但踩了会让你在组里很难做人。
我个人在实际操作中最大的体会是，VASP这东西的难点从来不在理论有多深，而在于每个环节都有一堆默认值和小陷阱，你只有真的跑过十几个体系、被各种报错教育过，才会形成肌肉记忆。上面这个MgO的例子虽然简单，但把它完整走通一遍，从建目录到判断收敛，再到用CONTCAR续算、做静态能量、画个能带或态密度，这条链路捋顺了，后面换成自己的体系就是换个POSCAR和POTCAR的事。真正值得花时间的其实是收敛测试和参数含义的理解，这两块扎实了，报错就成了查表的小问题。
