本脚本能够从单个 POSCAR 出发自动完成常压优化、压力点序列生成、结构迭代、能带与声子谱后处理的全流程。实现了从结构对称性分析到计算结果后处理的完整闭环适用于需要批量处理高压计算任务的体系。需提起配置好VASPKIT与PHONOPY环境。本文将详细介绍该脚本的设计思路、使用方法与注意事项准备文件仅需准备POSCAR操作流程$ sh db-4.sh[INFO] 2026-03-18 01:03:25 基础文件检查 [INFO] 2026-03-18 01:03:25 校验K点设置...[INFO] 2026-03-18 01:03:25 ✅ K点设置合理444 0 0[INFO] 2026-03-18 01:03:25 INCAR参数配置 是否设置磁性(y/n默认n)n是否取消对称性(y/n默认n)n是否增加截断能(y/n默认n)n自定义优化精度(y/n默认n格式1E-08;-2E-02)y输入EDIFF;EDIFFG1E-05;-2E-01添加NCORE(y/n默认n)yNCORE值4添加范德华修正(y/n默认n)nVASP运行命令示例mpirun -np 32 vasp_std默认mpirun -np 16 vasp_stdmpirun -np 32 vasp_std-w[INFO] 2026-03-18 01:04:05 ✅ VASP命令mpirun -np 32 vasp_std-w[INFO] 2026-03-18 01:04:05 初始结构设置 是否跳过常压0GPa优化(y/n默认n)n[INFO] 2026-03-18 01:04:18 常压0GPa结构优化 [INFO] 2026-03-18 01:04:18 ✅ 初始空间群International: Fm-3m[INFO] 2026-03-18 01:04:18 第1次常压优化...完成初次常压结构优化或可指定初始结构后可设置压力梯度与计算要求。[INFO] 2026-03-18 00:51:07 ✅ 已选择已有CONTCAR作为初始输入CONTCAR[INFO] 2026-03-18 00:51:07 ✅ 初始空间群基于已有CONTCARInternational: Fm-3m[INFO] 2026-03-18 00:51:07 压力梯度设置 压力输入类型1) 特殊压力点如1 5 10 GPa2) 常压至最大值如0-80GPa自动步长3) 分段压力如0-30,70-200GPa选择类型1/2/3默认11输入压力点空格分隔1 55[INFO] 2026-03-18 00:51:32 最终压力列表从小到大1 55 GPa[INFO] 2026-03-18 00:51:32 物性计算设置 是否计算能带(y/n默认n)y指定能带INCAR(y/n默认n)是否使用原胞进行能带计算(y/n默认n)y[INFO] 2026-03-18 00:51:40 ✅ 启用原胞进行能带计算[INFO] 2026-03-18 00:51:40 生成KPATH.invaspkit -task 303---------------------------------------------------------------| VASPKIT Standard Edition 1.5.1 (27 Jan. 2024) || Running VASPKIT Under Command-Line Mode |----------------------------------------------------------------------------------------- Summary ----------------------------Prototype: ABTotal Atoms in Input Cell: 8Lattice Constants in Input Cell: 6.411 6.411 6.411Lattice Angles in Input Cell: 90.000 90.000 90.000Total Atoms in Primitive Cell: 2Lattice Constants in Primitive Cell: 4.534 4.534 4.534Lattice Angles in Primitive Cell: 60.000 60.000 60.000Crystal System: CubicCrystal Class: m-3mBravais Lattice: cFExtended Bravais Lattice: cF2Space Group Number: 225Point Group: 32 [ Oh ]International: Fm-3mSymmetry Operations: 192Suggested K-Path: (shown in the next line)[ GAMMA-X-U|K-GAMMA-L-W-X ]----------------------------------------------------------------- (01) Written HIGH_SYMMETRY_POINTS File for Reference.-- (02) Written PRIMCELL.vasp file.-- (03) Written KPATH.in File for Band-Structure Calculation.o----------------------------WARNING----------------------------o| Do NOT forget to copy PRIMCELL.vasp to POSCAR unless you know || what you are doing. Otherwise you might get wrong results! |o---------------------------------------------------------------o---------------------------- Tip ------------------------------| Make sure you have installed python matplotlib successfully.|| The plot utility is under development in the preview release. || You CAN custom plot parameters in PLOT.in file if it exists, || or in ~/.vaspkit file if the PLOT.in file does NOT found. |----------------------------------------------------------------- (04) Reading Plot Setting from the PLOT.in File.findfont: Font family arial not found.fin-- (05) Written Brillouin_Zone_3D.jpg File.|---------------------------------------------------------------|| * CITATIONS * || When using VASPKIT in your research PLEASE cite the paper: || [1] V. WANG, N. XU, J.-C. LIU, G. TANG, W.-T. GENG, VASPKIT: A|| User-Friendly Interface Facilitating High-Throughput Computing|| and Analysis Using VASP Code, Computer Physics Communications || 267, 108033, (2021), DOI: 10.1016/j.cpc.2021.108033 |o---------------------------------------------------------------o[INFO] 2026-03-18 00:51:43 当前KPATH.in内容K-Path Generated by VASPKIT.20Line-ModeReciprocal0.0000000000 0.0000000000 0.0000000000 GAMMA0.5000000000 0.0000000000 0.5000000000 X0.5000000000 0.0000000000 0.5000000000 X0.6250000000 0.2500000000 0.6250000000 U0.3750000000 0.3750000000 0.7500000000 K0.0000000000 0.0000000000 0.0000000000 GAMMA0.0000000000 0.0000000000 0.0000000000 GAMMA0.5000000000 0.5000000000 0.5000000000 L0.5000000000 0.5000000000 0.5000000000 L0.5000000000 0.2500000000 0.7500000000 W0.5000000000 0.2500000000 0.7500000000 W0.5000000000 0.0000000000 0.5000000000 X确认KPATH.in是否合适(y/n默认y)y是否计算声子谱(y/n默认n)y指定声子INCAR(y/n默认n)[INFO] 2026-03-18 00:51:47 开始压力循环计算 [INFO] 2026-03-18 00:51:47 处理压力点1 GPa文件夹P1GPa [INFO] 2026-03-18 00:51:47 ✅ 设置PSTRESS10 kB1 GPa[INFO] 2026-03-18 00:51:47 第1次优化1 GPa结构...完整脚本附#!/bin/bash set -uo pipefail # 颜色定义 RED\033[0;31m GREEN\033[0;32m YELLOW\033[1;33m BLUE\033[0;34m NC\033[0m # 日志函数 log() { echo -e ${BLUE}[INFO] $(date %Y-%m-%d\ %H:%M:%S) $1${NC}; } error() { echo -e ${RED}[ERROR] $(date %Y-%m-%d\ %H:%M:%S) $1${NC}; exit 1; } warning() { echo -e ${YELLOW}[WARNING] $(date %Y-%m-%d\ %H:%M:%S) $1${NC}; } # 检查命令是否存在 check_command() { if ! command -v $1 /dev/null; then error 命令$1未找到请先安装并配置环境变量 fi } # 检查文件是否为非空白 check_non_empty_file() { local file$1 if [ ! -f $file ] || [ ! -s $file ]; then return 1 fi return 0 } # 核心函数提取空间群 get_space_group() { local poscar_file$1 if ! check_non_empty_file $poscar_file; then echo UNKNOWN return fi # 临时复制POSCAR到当前目录执行601 cp $poscar_file TMP_POSCAR mv POSCAR POSCAR_BAK 2/dev/null mv TMP_POSCAR POSCAR # 执行601并提取International空间群 space_group$(vaspkit -task 601 21 | grep International: | awk {print $2}) # 恢复原POSCAR mv POSCAR_BAK POSCAR 2/dev/null rm -f TMP_POSCAR [ -z $space_group ] space_groupUNKNOWN echo $space_group } # 第一步基础文件检查与修复 log 基础文件检查 # 2. 自动生成KPOINTS缺失/过小 if ! check_non_empty_file KPOINTS; then warning 未找到有效KPOINTS自动生成... (echo 102; echo 2; echo 0.04) | vaspkit fi # 3. 修正K点校验核心提取第4行 log 校验K点设置... kmesh_line$(sed -n 4p KPOINTS | awk {print $1, $2, $3} | tr -d | awk {for(i1;i3;i) printf %.0f , $i}) k1$(echo $kmesh_line | awk {print $1}) k2$(echo $kmesh_line | awk {print $2}) k3$(echo $kmesh_line | awk {print $3}) if [ $k1 -eq 1 ] [ $k2 -eq 1 ] [ $k3 -eq 1 ]; then warning K点过小1x1x1自动更新... (echo 102; echo 2; echo 0.04) | vaspkit kmesh_line$(sed -n 4p KPOINTS | awk {print $1, $2, $3}) fi log ✅ K点设置合理$kmesh_line # 1. 检查核心文件仅检查INCAR/POTCARPOSCAR可选因为可能用已有CONTCAR if ! check_non_empty_file INCAR; then error INCAR文件不存在或为空; fi if ! check_non_empty_file POTCAR; then error POTCAR文件不存在或为空需手动从VASP官网获取; fi # 第二步INCAR参数配置 log INCAR参数配置 cp INCAR INCAR.bak # 1. 磁性设置 read -p 是否设置磁性(y/n默认n) is_mag; is_mag${is_mag:-n} sed -i /ISPIN/d INCAR echo ISPIN $( [ $is_mag y ] echo 2 || echo 1 ) INCAR [ $is_mag y ] warning 已设ISPIN2请手动添加MAGMOM # 2. 对称性设置 read -p 是否取消对称性(y/n默认n) is_sym_off; is_sym_off${is_sym_off:-n} sed -i /ISYM/d INCAR echo ISYM $( [ $is_sym_off y ] echo 0 || echo 2 ) INCAR # 3. 截断能 read -p 是否增加截断能(y/n默认n) is_encut; is_encut${is_encut:-n} sed -i /ENCUT/d INCAR if [ $is_encut y ]; then while true; do read -p 输入截断能400整数 encut; if [[ $encut ~ ^[0-9]$ ]] [ $encut -gt 400 ]; then break; fi warning 输入无效 done else encut400 fi echo ENCUT $encut INCAR # 4. 优化精度 read -p 自定义优化精度(y/n默认n格式1E-08;-2E-02) is_ediff; is_ediff${is_ediff:-n} sed -i /EDIFF/d INCAR; sed -i /EDIFFG/d INCAR if [ $is_ediff y ]; then read -p 输入EDIFF;EDIFFG ediff_vals ediff$(echo $ediff_vals | cut -d; -f1) ediffg$(echo $ediff_vals | cut -d; -f2) else ediff1E-05; ediffg-2E-01 fi echo EDIFF $ediff INCAR; echo EDIFFG $ediffg INCAR # 5. NCORE/范德华修正可选跳过 read -p 添加NCORE(y/n默认n) is_ncore; is_ncore${is_ncore:-n} [ $is_ncore y ] read -p NCORE值 ncore sed -i /NCORE/d INCAR echo NCORE $ncore INCAR read -p 添加范德华修正(y/n默认n) is_vdw; is_vdw${is_vdw:-n} [ $is_vdw y ] read -p IVDW值1/2/10等 ivdw sed -i /IVDW/d INCAR echo IVDW $ivdw INCAR # 6. 强制参数 sed -i /ISTART/s/^/#/ INCAR; sed -i /ICHARG/s/^/#/ INCAR sed -i /LWAVE/d INCAR; sed -i /LCHARG/d INCAR echo LWAVE .FALSE. INCAR; echo LCHARG .FALSE. INCAR sed -i /NSW/d INCAR; sed -i /IBRION/d INCAR; sed -i /ISIF/d INCAR echo NSW 100 INCAR; echo IBRION 2 INCAR; echo ISIF 3 INCAR # 第三步VASP运行命令 read -p VASP运行命令示例mpirun -np 32 vasp_std默认mpirun -np 16 vasp_std vasp_cmd vasp_cmd${vasp_cmd:-mpirun -np 16 vasp_std} log ✅ VASP命令$vasp_cmd # 第四步初始CONTCAR设置核心支持跳过常压优化 log 初始结构设置 current_contcar read -p 是否跳过常压0GPa优化(y/n默认n) skip_0gpa; skip_0gpa${skip_0gpa:-n} if [ $skip_0gpa y ]; then # 跳过常压优化询问已有CONTCAR路径 while true; do read -p 请输入已有优化结果的CONTCAR文件路径如./CONTCAR.0GPa 或 ./P5GPa/CONTCAR.5GPa contcar_path if check_non_empty_file $contcar_path; then current_contcar$contcar_path # 备份初始CONTCAR cp $current_contcar CONTCAR.initial current_contcarCONTCAR.initial log ✅ 已选择已有CONTCAR作为初始输入$contcar_path # 提取初始空间群 initial_space_group$(get_space_group $current_contcar) log ✅ 初始空间群基于已有CONTCARInternational: $initial_space_group break else warning 指定的CONTCAR文件不存在或为空请重新输入 fi done else # 执行常压0GPa优化 log 常压0GPa结构优化 # 检查POSCAR常压优化必须有 if ! check_non_empty_file POSCAR; then error 常压优化需要有效POSCAR文件 fi # 提取初始空间群 initial_space_group$(get_space_group POSCAR) log ✅ 初始空间群International: $initial_space_group # 压力设置0GPa sed -i /PSTRESS/d INCAR; echo PSTRESS 0 INCAR optimization_donefalse max_attempts5; attempt0 while [ $attempt -lt $max_attempts ] [ $optimization_done false ]; do attempt$((attempt 1)) log 第$attempt次常压优化... log_filevasp_0GPa_$attempt.log if ! $vasp_cmd $log_file 21; then warning VASP运行失败日志$log_file read -p 重试(y/n默认n) retry; retry${retry:-n} [ $retry n ] error 用户终止 continue fi # 检查CONTCAR是否有效核心不管收敛与否优先用CONTCAR if check_non_empty_file CONTCAR; then cp CONTCAR CONTCAR.0GPa # 备份0GPa的CONTCAR current_contcarCONTCAR.0GPa log ✅ 生成有效CONTCAR0GPa强制使用该文件作为后续输入 optimization_donetrue else warning 未生成有效CONTCAR重试... continue fi # 可选检查收敛条件仅日志提示不影响CONTCAR使用 ionic_steps$(grep -E ^ [0-9] F $log_file | wc -l) reached$(grep -c reached required accuracy $log_file) if [ $ionic_steps -lt 3 ] [ $reached -ge 1 ]; then log ✅ 常压优化收敛离子步数$ionic_steps else warning 常压优化未收敛但仍使用CONTCAR作为后续输入 fi done [ $optimization_done false ] error 常压优化未生成有效CONTCAR fi # 第五步压力梯度设置排序校验 log 压力梯度设置 echo 压力输入类型 echo 1) 特殊压力点如1 5 10 GPa echo 2) 常压至最大值如0-80GPa自动步长 echo 3) 分段压力如0-30,70-200GPa read -p 选择类型1/2/3默认1 p_type; p_type${p_type:-1} pressure_list() case $p_type in 1) read -p 输入压力点空格分隔 p_points pressure_list($(echo $p_points | tr \n | sort -n | uniq)) ;; 2) read -p 压力最大值GPa max_p # 自动步长0-101, 10-202, 20-505, 5010 for p in $(seq 0 1 10); do [ $p -le $max_p ] pressure_list($p); done for p in $(seq 12 2 20); do [ $p -le $max_p ] pressure_list($p); done for p in $(seq 25 5 50); do [ $p -le $max_p ] pressure_list($p); done p60; while [ $p -le $max_p ]; do pressure_list($p); p$((p10)); done pressure_list($(echo ${pressure_list[]} | tr \n | sort -n | uniq)) ;; 3) read -p 分段压力逗号分隔如0-30,70-200 p_seg read -p 步长auto/数字默认auto step; step${step:-auto} IFS, read -ra segs $p_seg for seg in ${segs[]}; do start$(echo $seg | cut -d- -f1 | awk {print int($1)}) end$(echo $seg | cut -d- -f2 | awk {print int($1)}) [ $step auto ] step$( [ $end -le 10 ] echo 1 || [ $end -le 20 ] echo 2 || [ $end -le 50 ] echo 5 || echo 10 ) for p in $(seq $start $step $end); do pressure_list($p); done done pressure_list($(echo ${pressure_list[]} | tr \n | sort -n | uniq)) ;; esac # 移除初始压力如果跳过常压优化初始压力可能不是0需手动确认 if [ $skip_0gpa n ]; then # 跳过常压优化时移除0GPa已计算 pressure_list($(echo ${pressure_list[]} | tr \n | grep -v ^0$ | sort -n)) fi log 最终压力列表从小到大${pressure_list[*]} GPa # 第六步能带/声子谱设置原胞选项 log 物性计算设置 # 能带计算原胞选项 read -p 是否计算能带(y/n默认n) do_band; do_band${do_band:-n} band_incar use_prim_cellfalse kpath_updatedfalse if [ $do_band y ]; then read -p 指定能带INCAR(y/n默认n) has_band_incar; has_band_incar${has_band_incar:-n} [ $has_band_incar y ] read -p 能带INCAR路径 band_incar [ ! -f $band_incar ] error 能带INCAR不存在 # 询问是否用原胞计算能带 read -p 是否使用原胞进行能带计算(y/n默认n) use_prim; use_prim${use_prim:-n} [ $use_prim y ] use_prim_celltrue log ✅ 启用原胞进行能带计算 # 生成并输出KPATH.in到终端 log 生成KPATH.invaspkit -task 303 # 临时使用初始CONTCAR生成KPATH.in cp $current_contcar TMP_POSCAR mv POSCAR POSCAR_BAK 2/dev/null mv TMP_POSCAR POSCAR vaspkit -task 303 mv POSCAR_BAK POSCAR 2/dev/null rm -f TMP_POSCAR log 当前KPATH.in内容 cat KPATH.in read -p 确认KPATH.in是否合适(y/n默认y) confirm_kpath; confirm_kpath${confirm_kpath:-y} [ $confirm_kpath n ] read -p 请输入自定义KPATH.in路径 kpath_file cp $kpath_file KPATH.in fi # 声子谱计算 read -p 是否计算声子谱(y/n默认n) do_phono; do_phono${do_phono:-n} phono_incar if [ $do_phono y ]; then check_command phonopy read -p 指定声子INCAR(y/n默认n) has_phono_incar; has_phono_incar${has_phono_incar:-n} [ $has_phono_incar y ] read -p 声子INCAR路径 phono_incar [ ! -f $phono_incar ] error 声子INCAR不存在 fi # 第七步压力循环计算核心强制用CONTCAR log 开始压力循环计算 for p in ${pressure_list[]}; do dirP${p}GPa mkdir -p $dir log 处理压力点$p GPa文件夹$dir # 1. 复制初始/前一个压力点的CONTCAR核心强制用CONTCAR if ! check_non_empty_file $current_contcar; then warning 初始/前一个CONTCAR($current_contcar)无效请检查 read -p 是否终止计算(y/n默认y) stop; stop${stop:-y} [ $stop y ] error 用户终止计算 || continue fi cp $current_contcar $dir/POSCAR # 作为当前压力点的输入POSCAR cp KPOINTS POTCAR INCAR $dir/ # 2. 修改当前压力的PSTRESS p_kb$(echo $p * 10 | bc) sed -i s/PSTRESS .*/PSTRESS $p_kb/ $dir/INCAR log ✅ 设置PSTRESS$p_kb kB$p GPa # 3. 进入文件夹优化结构 cd $dir || { warning 无法进入$dir跳过; continue; } # 结构优化循环强制用CONTCAR opt_donefalse; attempt0; max_attempts3 while [ $attempt -lt $max_attempts ] [ $opt_done false ]; do attempt$((attempt 1)) log 第$attempt次优化$p GPa结构... log_filevasp_${p}GPa_$attempt.log if ! $vasp_cmd $log_file 21; then warning $p GPa第$attempt次VASP运行失败日志$log_file continue fi # 核心检查CONTCAR是否有效不管收敛与否 if check_non_empty_file CONTCAR; then cp CONTCAR CONTCAR.${p}GPa # 备份当前压力点的CONTCAR current_contcar../$dir/CONTCAR.${p}GPa # 更新下一个压力点的输入 log ✅ $p GPa生成有效CONTCAR强制使用该文件作为后续输入 opt_donetrue # 可选检查收敛条件仅日志提示 ionic_steps$(grep -E ^ [0-9] F $log_file | wc -l) reached$(grep -c reached required accuracy $log_file) if [ $ionic_steps -lt 3 ] [ $reached -ge 1 ]; then log ✅ $p GPa优化收敛离子步数$ionic_steps else warning $p GPa优化未收敛但仍使用CONTCAR作为后续输入 fi # 检查空间群是否变化仅能带计算时 if [ $do_band y ] [ $kpath_updated false ]; then current_sg$(get_space_group CONTCAR) if [ $current_sg ! $initial_space_group ] [ $current_sg ! UNKNOWN ]; then warning 空间群变化初始$initial_space_group → 当前$current_sg read -p 是否更新KPATH.in(y/n默认n) update_kpath; update_kpath${update_kpath:-n} if [ $update_kpath y ]; then # 基于当前CONTCAR生成新的KPATH.in vaspkit -task 303 cp KPATH.in ../../KPATH.in kpath_updatedtrue log ✅ 已基于当前结构更新KPATH.in fi fi fi else warning $p GPa第$attempt次优化未生成有效CONTCAR重试... continue fi done [ $opt_done false ] { warning $p GPa优化失败跳过; cd ..; continue; } # 能带计算原胞选项 if [ $do_band y ]; then log 开始$p GPa能带计算... mkdir -p band cd band || { warning 无法进入band跳过; cd ..; continue; } # 复制当前压力点的CONTCAR到band目录 cp ../CONTCAR POSCAR # 如果使用原胞先生成PRIMCELL.vasp并替换POSCAR再生成K点 if [ $use_prim_cell true ]; then log 生成原胞PRIMCELL.vaspvaspkit -task 602 vaspkit -task 602 if check_non_empty_file PRIMCELL.vasp; then mv POSCAR POSCAR_ori # 备份原结构 cp PRIMCELL.vasp POSCAR log ✅ 已替换为原胞PRIMCELL.vasp作为能带计算输入 # 重新生成K点基于原胞 log 基于原胞生成新的K点 (echo 102; echo 2; echo 0.04) | vaspkit else warning 生成原胞失败使用原CONTCAR继续 fi fi # 复制其他文件 cp ../../POTCAR ../../INCAR . # 处理能带INCAR if [ -n $band_incar ]; then cp $band_incar INCAR else sed -i /NSW/s/^/#/ INCAR; sed -i /IBRION/s/^/#/ INCAR sed -i /LWAVE/d INCAR; sed -i /LCHARG/d INCAR echo LWAVE .TRUE. INCAR; echo LCHARG .TRUE. INCAR echo LORBIT 10 INCAR; echo NEDOS 5001 INCAR fi # 第一步SCF计算 log 执行第一步SCF自洽计算 $vasp_cmd band_calc1.log 21 # 第二步能带计算替换KPATH.in cp ../../KPATH.in KPOINTS log 执行第二步能带计算 $vasp_cmd band_calc2.log 21 # 数据处理 log 处理能带数据 (echo 211; echo 1)|vaspkit (echo 215; echo 1)|vaspkit (echo 111; echo 2)|vaspkit (echo 211; echo 2)|vaspkit (echo 113; echo all;echo 2)|vaspkit (echo 213; echo 2;echo all)|vaspkit cd .. fi # 声子谱计算强制用CONTCAR if [ $do_phono y ]; then log 开始$p GPa声子谱计算... mkdir -p phono cd phono || { warning 无法进入phono跳过; cd ..; continue; } # 复制当前压力点的CONTCAR作为输入 cp ../CONTCAR POSCAR cp ../../POTCAR ../../KPOINTS ../../INCAR . # 处理声子INCAR if [ -n $phono_incar ]; then cp $phono_incar INCAR else sed -i /NCORE/s/^/#/ INCAR sed -i /NSW/d INCAR; sed -i /IBRION/d INCAR echo NSW 1 INCAR grep -q IVDW INCAR echo IBRION 6 INCAR || echo IBRION 8 INCAR fi # 运行VASPphonopy log 执行声子谱VASP计算 $vasp_cmd phono_vasp.log 21 if check_non_empty_file vasprun.xml; then # 生成KPATH.phonopy cat KPATH.phonopy EOF NPOINTS 501 DIM 1 1 1 BAND AUTO BAND_LABELS AUTO TETRAHEDRON .TRUE. PDOS AUTO BAND_CONNECTION .TRUE. FORCE_CONSTANTS READ EOF # 运行phonopy log 执行phonopy计算 phonopy --fc vasprun.xml phonopy -s -p KPATH.phonopy --paauto phonopy-bandplot --gnuplot | tee phono.txt rm -f *.yaml log ✅ 声子谱计算完成 else warning 未生成有效vasprun.xml跳过phonopy计算 fi cd .. fi cd .. # 返回主目录 done # 计算完成 log ${GREEN} 所有计算完成 ${NC} log ✅ 初始CONTCAR$current_contcar跳过常压优化/ CONTCAR.0GPa常压优化 log ✅ 压力点CONTCARP*GPa/CONTCAR.*GPa每个压力点的核心输出 log ✅ 关键文件 log - 能带计算P*GPa/band/使用CONTCAR/原胞PRIMCELL.vasp作为输入 log - 声子谱计算P*GPa/phono/使用CONTCAR作为输入 log ✅ 注意所有后续计算均强制使用CONTCAR无论是否收敛
VASP 自动压力梯度能带声子谱计算自动化脚本
发布时间：2026/6/6 20:15:53
本脚本能够从单个 POSCAR 出发自动完成常压优化、压力点序列生成、结构迭代、能带与声子谱后处理的全流程。实现了从结构对称性分析到计算结果后处理的完整闭环适用于需要批量处理高压计算任务的体系。需提起配置好VASPKIT与PHONOPY环境。本文将详细介绍该脚本的设计思路、使用方法与注意事项准备文件仅需准备POSCAR操作流程$ sh db-4.sh[INFO] 2026-03-18 01:03:25 基础文件检查 [INFO] 2026-03-18 01:03:25 校验K点设置...[INFO] 2026-03-18 01:03:25 ✅ K点设置合理444 0 0[INFO] 2026-03-18 01:03:25 INCAR参数配置 是否设置磁性(y/n默认n)n是否取消对称性(y/n默认n)n是否增加截断能(y/n默认n)n自定义优化精度(y/n默认n格式1E-08;-2E-02)y输入EDIFF;EDIFFG1E-05;-2E-01添加NCORE(y/n默认n)yNCORE值4添加范德华修正(y/n默认n)nVASP运行命令示例mpirun -np 32 vasp_std默认mpirun -np 16 vasp_stdmpirun -np 32 vasp_std-w[INFO] 2026-03-18 01:04:05 ✅ VASP命令mpirun -np 32 vasp_std-w[INFO] 2026-03-18 01:04:05 初始结构设置 是否跳过常压0GPa优化(y/n默认n)n[INFO] 2026-03-18 01:04:18 常压0GPa结构优化 [INFO] 2026-03-18 01:04:18 ✅ 初始空间群International: Fm-3m[INFO] 2026-03-18 01:04:18 第1次常压优化...完成初次常压结构优化或可指定初始结构后可设置压力梯度与计算要求。[INFO] 2026-03-18 00:51:07 ✅ 已选择已有CONTCAR作为初始输入CONTCAR[INFO] 2026-03-18 00:51:07 ✅ 初始空间群基于已有CONTCARInternational: Fm-3m[INFO] 2026-03-18 00:51:07 压力梯度设置 压力输入类型1) 特殊压力点如1 5 10 GPa2) 常压至最大值如0-80GPa自动步长3) 分段压力如0-30,70-200GPa选择类型1/2/3默认11输入压力点空格分隔1 55[INFO] 2026-03-18 00:51:32 最终压力列表从小到大1 55 GPa[INFO] 2026-03-18 00:51:32 物性计算设置 是否计算能带(y/n默认n)y指定能带INCAR(y/n默认n)是否使用原胞进行能带计算(y/n默认n)y[INFO] 2026-03-18 00:51:40 ✅ 启用原胞进行能带计算[INFO] 2026-03-18 00:51:40 生成KPATH.invaspkit -task 303---------------------------------------------------------------| VASPKIT Standard Edition 1.5.1 (27 Jan. 2024) || Running VASPKIT Under Command-Line Mode |----------------------------------------------------------------------------------------- Summary ----------------------------Prototype: ABTotal Atoms in Input Cell: 8Lattice Constants in Input Cell: 6.411 6.411 6.411Lattice Angles in Input Cell: 90.000 90.000 90.000Total Atoms in Primitive Cell: 2Lattice Constants in Primitive Cell: 4.534 4.534 4.534Lattice Angles in Primitive Cell: 60.000 60.000 60.000Crystal System: CubicCrystal Class: m-3mBravais Lattice: cFExtended Bravais Lattice: cF2Space Group Number: 225Point Group: 32 [ Oh ]International: Fm-3mSymmetry Operations: 192Suggested K-Path: (shown in the next line)[ GAMMA-X-U|K-GAMMA-L-W-X ]----------------------------------------------------------------- (01) Written HIGH_SYMMETRY_POINTS File for Reference.-- (02) Written PRIMCELL.vasp file.-- (03) Written KPATH.in File for Band-Structure Calculation.o----------------------------WARNING----------------------------o| Do NOT forget to copy PRIMCELL.vasp to POSCAR unless you know || what you are doing. Otherwise you might get wrong results! |o---------------------------------------------------------------o---------------------------- Tip ------------------------------| Make sure you have installed python matplotlib successfully.|| The plot utility is under development in the preview release. || You CAN custom plot parameters in PLOT.in file if it exists, || or in ~/.vaspkit file if the PLOT.in file does NOT found. |----------------------------------------------------------------- (04) Reading Plot Setting from the PLOT.in File.findfont: Font family arial not found.fin-- (05) Written Brillouin_Zone_3D.jpg File.|---------------------------------------------------------------|| * CITATIONS * || When using VASPKIT in your research PLEASE cite the paper: || [1] V. WANG, N. XU, J.-C. LIU, G. TANG, W.-T. GENG, VASPKIT: A|| User-Friendly Interface Facilitating High-Throughput Computing|| and Analysis Using VASP Code, Computer Physics Communications || 267, 108033, (2021), DOI: 10.1016/j.cpc.2021.108033 |o---------------------------------------------------------------o[INFO] 2026-03-18 00:51:43 当前KPATH.in内容K-Path Generated by VASPKIT.20Line-ModeReciprocal0.0000000000 0.0000000000 0.0000000000 GAMMA0.5000000000 0.0000000000 0.5000000000 X0.5000000000 0.0000000000 0.5000000000 X0.6250000000 0.2500000000 0.6250000000 U0.3750000000 0.3750000000 0.7500000000 K0.0000000000 0.0000000000 0.0000000000 GAMMA0.0000000000 0.0000000000 0.0000000000 GAMMA0.5000000000 0.5000000000 0.5000000000 L0.5000000000 0.5000000000 0.5000000000 L0.5000000000 0.2500000000 0.7500000000 W0.5000000000 0.2500000000 0.7500000000 W0.5000000000 0.0000000000 0.5000000000 X确认KPATH.in是否合适(y/n默认y)y是否计算声子谱(y/n默认n)y指定声子INCAR(y/n默认n)[INFO] 2026-03-18 00:51:47 开始压力循环计算 [INFO] 2026-03-18 00:51:47 处理压力点1 GPa文件夹P1GPa [INFO] 2026-03-18 00:51:47 ✅ 设置PSTRESS10 kB1 GPa[INFO] 2026-03-18 00:51:47 第1次优化1 GPa结构...完整脚本附#!/bin/bash set -uo pipefail # 颜色定义 RED\033[0;31m GREEN\033[0;32m YELLOW\033[1;33m BLUE\033[0;34m NC\033[0m # 日志函数 log() { echo -e ${BLUE}[INFO] $(date %Y-%m-%d\ %H:%M:%S) $1${NC}; } error() { echo -e ${RED}[ERROR] $(date %Y-%m-%d\ %H:%M:%S) $1${NC}; exit 1; } warning() { echo -e ${YELLOW}[WARNING] $(date %Y-%m-%d\ %H:%M:%S) $1${NC}; } # 检查命令是否存在 check_command() { if ! command -v $1 /dev/null; then error 命令$1未找到请先安装并配置环境变量 fi } # 检查文件是否为非空白 check_non_empty_file() { local file$1 if [ ! -f $file ] || [ ! -s $file ]; then return 1 fi return 0 } # 核心函数提取空间群 get_space_group() { local poscar_file$1 if ! check_non_empty_file $poscar_file; then echo UNKNOWN return fi # 临时复制POSCAR到当前目录执行601 cp $poscar_file TMP_POSCAR mv POSCAR POSCAR_BAK 2/dev/null mv TMP_POSCAR POSCAR # 执行601并提取International空间群 space_group$(vaspkit -task 601 21 | grep International: | awk {print $2}) # 恢复原POSCAR mv POSCAR_BAK POSCAR 2/dev/null rm -f TMP_POSCAR [ -z $space_group ] space_groupUNKNOWN echo $space_group } # 第一步基础文件检查与修复 log 基础文件检查 # 2. 自动生成KPOINTS缺失/过小 if ! check_non_empty_file KPOINTS; then warning 未找到有效KPOINTS自动生成... (echo 102; echo 2; echo 0.04) | vaspkit fi # 3. 修正K点校验核心提取第4行 log 校验K点设置... kmesh_line$(sed -n 4p KPOINTS | awk {print $1, $2, $3} | tr -d | awk {for(i1;i3;i) printf %.0f , $i}) k1$(echo $kmesh_line | awk {print $1}) k2$(echo $kmesh_line | awk {print $2}) k3$(echo $kmesh_line | awk {print $3}) if [ $k1 -eq 1 ] [ $k2 -eq 1 ] [ $k3 -eq 1 ]; then warning K点过小1x1x1自动更新... (echo 102; echo 2; echo 0.04) | vaspkit kmesh_line$(sed -n 4p KPOINTS | awk {print $1, $2, $3}) fi log ✅ K点设置合理$kmesh_line # 1. 检查核心文件仅检查INCAR/POTCARPOSCAR可选因为可能用已有CONTCAR if ! check_non_empty_file INCAR; then error INCAR文件不存在或为空; fi if ! check_non_empty_file POTCAR; then error POTCAR文件不存在或为空需手动从VASP官网获取; fi # 第二步INCAR参数配置 log INCAR参数配置 cp INCAR INCAR.bak # 1. 磁性设置 read -p 是否设置磁性(y/n默认n) is_mag; is_mag${is_mag:-n} sed -i /ISPIN/d INCAR echo ISPIN $( [ $is_mag y ] echo 2 || echo 1 ) INCAR [ $is_mag y ] warning 已设ISPIN2请手动添加MAGMOM # 2. 对称性设置 read -p 是否取消对称性(y/n默认n) is_sym_off; is_sym_off${is_sym_off:-n} sed -i /ISYM/d INCAR echo ISYM $( [ $is_sym_off y ] echo 0 || echo 2 ) INCAR # 3. 截断能 read -p 是否增加截断能(y/n默认n) is_encut; is_encut${is_encut:-n} sed -i /ENCUT/d INCAR if [ $is_encut y ]; then while true; do read -p 输入截断能400整数 encut; if [[ $encut ~ ^[0-9]$ ]] [ $encut -gt 400 ]; then break; fi warning 输入无效 done else encut400 fi echo ENCUT $encut INCAR # 4. 优化精度 read -p 自定义优化精度(y/n默认n格式1E-08;-2E-02) is_ediff; is_ediff${is_ediff:-n} sed -i /EDIFF/d INCAR; sed -i /EDIFFG/d INCAR if [ $is_ediff y ]; then read -p 输入EDIFF;EDIFFG ediff_vals ediff$(echo $ediff_vals | cut -d; -f1) ediffg$(echo $ediff_vals | cut -d; -f2) else ediff1E-05; ediffg-2E-01 fi echo EDIFF $ediff INCAR; echo EDIFFG $ediffg INCAR # 5. NCORE/范德华修正可选跳过 read -p 添加NCORE(y/n默认n) is_ncore; is_ncore${is_ncore:-n} [ $is_ncore y ] read -p NCORE值 ncore sed -i /NCORE/d INCAR echo NCORE $ncore INCAR read -p 添加范德华修正(y/n默认n) is_vdw; is_vdw${is_vdw:-n} [ $is_vdw y ] read -p IVDW值1/2/10等 ivdw sed -i /IVDW/d INCAR echo IVDW $ivdw INCAR # 6. 强制参数 sed -i /ISTART/s/^/#/ INCAR; sed -i /ICHARG/s/^/#/ INCAR sed -i /LWAVE/d INCAR; sed -i /LCHARG/d INCAR echo LWAVE .FALSE. INCAR; echo LCHARG .FALSE. INCAR sed -i /NSW/d INCAR; sed -i /IBRION/d INCAR; sed -i /ISIF/d INCAR echo NSW 100 INCAR; echo IBRION 2 INCAR; echo ISIF 3 INCAR # 第三步VASP运行命令 read -p VASP运行命令示例mpirun -np 32 vasp_std默认mpirun -np 16 vasp_std vasp_cmd vasp_cmd${vasp_cmd:-mpirun -np 16 vasp_std} log ✅ VASP命令$vasp_cmd # 第四步初始CONTCAR设置核心支持跳过常压优化 log 初始结构设置 current_contcar read -p 是否跳过常压0GPa优化(y/n默认n) skip_0gpa; skip_0gpa${skip_0gpa:-n} if [ $skip_0gpa y ]; then # 跳过常压优化询问已有CONTCAR路径 while true; do read -p 请输入已有优化结果的CONTCAR文件路径如./CONTCAR.0GPa 或 ./P5GPa/CONTCAR.5GPa contcar_path if check_non_empty_file $contcar_path; then current_contcar$contcar_path # 备份初始CONTCAR cp $current_contcar CONTCAR.initial current_contcarCONTCAR.initial log ✅ 已选择已有CONTCAR作为初始输入$contcar_path # 提取初始空间群 initial_space_group$(get_space_group $current_contcar) log ✅ 初始空间群基于已有CONTCARInternational: $initial_space_group break else warning 指定的CONTCAR文件不存在或为空请重新输入 fi done else # 执行常压0GPa优化 log 常压0GPa结构优化 # 检查POSCAR常压优化必须有 if ! check_non_empty_file POSCAR; then error 常压优化需要有效POSCAR文件 fi # 提取初始空间群 initial_space_group$(get_space_group POSCAR) log ✅ 初始空间群International: $initial_space_group # 压力设置0GPa sed -i /PSTRESS/d INCAR; echo PSTRESS 0 INCAR optimization_donefalse max_attempts5; attempt0 while [ $attempt -lt $max_attempts ] [ $optimization_done false ]; do attempt$((attempt 1)) log 第$attempt次常压优化... log_filevasp_0GPa_$attempt.log if ! $vasp_cmd $log_file 21; then warning VASP运行失败日志$log_file read -p 重试(y/n默认n) retry; retry${retry:-n} [ $retry n ] error 用户终止 continue fi # 检查CONTCAR是否有效核心不管收敛与否优先用CONTCAR if check_non_empty_file CONTCAR; then cp CONTCAR CONTCAR.0GPa # 备份0GPa的CONTCAR current_contcarCONTCAR.0GPa log ✅ 生成有效CONTCAR0GPa强制使用该文件作为后续输入 optimization_donetrue else warning 未生成有效CONTCAR重试... continue fi # 可选检查收敛条件仅日志提示不影响CONTCAR使用 ionic_steps$(grep -E ^ [0-9] F $log_file | wc -l) reached$(grep -c reached required accuracy $log_file) if [ $ionic_steps -lt 3 ] [ $reached -ge 1 ]; then log ✅ 常压优化收敛离子步数$ionic_steps else warning 常压优化未收敛但仍使用CONTCAR作为后续输入 fi done [ $optimization_done false ] error 常压优化未生成有效CONTCAR fi # 第五步压力梯度设置排序校验 log 压力梯度设置 echo 压力输入类型 echo 1) 特殊压力点如1 5 10 GPa echo 2) 常压至最大值如0-80GPa自动步长 echo 3) 分段压力如0-30,70-200GPa read -p 选择类型1/2/3默认1 p_type; p_type${p_type:-1} pressure_list() case $p_type in 1) read -p 输入压力点空格分隔 p_points pressure_list($(echo $p_points | tr \n | sort -n | uniq)) ;; 2) read -p 压力最大值GPa max_p # 自动步长0-101, 10-202, 20-505, 5010 for p in $(seq 0 1 10); do [ $p -le $max_p ] pressure_list($p); done for p in $(seq 12 2 20); do [ $p -le $max_p ] pressure_list($p); done for p in $(seq 25 5 50); do [ $p -le $max_p ] pressure_list($p); done p60; while [ $p -le $max_p ]; do pressure_list($p); p$((p10)); done pressure_list($(echo ${pressure_list[]} | tr \n | sort -n | uniq)) ;; 3) read -p 分段压力逗号分隔如0-30,70-200 p_seg read -p 步长auto/数字默认auto step; step${step:-auto} IFS, read -ra segs $p_seg for seg in ${segs[]}; do start$(echo $seg | cut -d- -f1 | awk {print int($1)}) end$(echo $seg | cut -d- -f2 | awk {print int($1)}) [ $step auto ] step$( [ $end -le 10 ] echo 1 || [ $end -le 20 ] echo 2 || [ $end -le 50 ] echo 5 || echo 10 ) for p in $(seq $start $step $end); do pressure_list($p); done done pressure_list($(echo ${pressure_list[]} | tr \n | sort -n | uniq)) ;; esac # 移除初始压力如果跳过常压优化初始压力可能不是0需手动确认 if [ $skip_0gpa n ]; then # 跳过常压优化时移除0GPa已计算 pressure_list($(echo ${pressure_list[]} | tr \n | grep -v ^0$ | sort -n)) fi log 最终压力列表从小到大${pressure_list[*]} GPa # 第六步能带/声子谱设置原胞选项 log 物性计算设置 # 能带计算原胞选项 read -p 是否计算能带(y/n默认n) do_band; do_band${do_band:-n} band_incar use_prim_cellfalse kpath_updatedfalse if [ $do_band y ]; then read -p 指定能带INCAR(y/n默认n) has_band_incar; has_band_incar${has_band_incar:-n} [ $has_band_incar y ] read -p 能带INCAR路径 band_incar [ ! -f $band_incar ] error 能带INCAR不存在 # 询问是否用原胞计算能带 read -p 是否使用原胞进行能带计算(y/n默认n) use_prim; use_prim${use_prim:-n} [ $use_prim y ] use_prim_celltrue log ✅ 启用原胞进行能带计算 # 生成并输出KPATH.in到终端 log 生成KPATH.invaspkit -task 303 # 临时使用初始CONTCAR生成KPATH.in cp $current_contcar TMP_POSCAR mv POSCAR POSCAR_BAK 2/dev/null mv TMP_POSCAR POSCAR vaspkit -task 303 mv POSCAR_BAK POSCAR 2/dev/null rm -f TMP_POSCAR log 当前KPATH.in内容 cat KPATH.in read -p 确认KPATH.in是否合适(y/n默认y) confirm_kpath; confirm_kpath${confirm_kpath:-y} [ $confirm_kpath n ] read -p 请输入自定义KPATH.in路径 kpath_file cp $kpath_file KPATH.in fi # 声子谱计算 read -p 是否计算声子谱(y/n默认n) do_phono; do_phono${do_phono:-n} phono_incar if [ $do_phono y ]; then check_command phonopy read -p 指定声子INCAR(y/n默认n) has_phono_incar; has_phono_incar${has_phono_incar:-n} [ $has_phono_incar y ] read -p 声子INCAR路径 phono_incar [ ! -f $phono_incar ] error 声子INCAR不存在 fi # 第七步压力循环计算核心强制用CONTCAR log 开始压力循环计算 for p in ${pressure_list[]}; do dirP${p}GPa mkdir -p $dir log 处理压力点$p GPa文件夹$dir # 1. 复制初始/前一个压力点的CONTCAR核心强制用CONTCAR if ! check_non_empty_file $current_contcar; then warning 初始/前一个CONTCAR($current_contcar)无效请检查 read -p 是否终止计算(y/n默认y) stop; stop${stop:-y} [ $stop y ] error 用户终止计算 || continue fi cp $current_contcar $dir/POSCAR # 作为当前压力点的输入POSCAR cp KPOINTS POTCAR INCAR $dir/ # 2. 修改当前压力的PSTRESS p_kb$(echo $p * 10 | bc) sed -i s/PSTRESS .*/PSTRESS $p_kb/ $dir/INCAR log ✅ 设置PSTRESS$p_kb kB$p GPa # 3. 进入文件夹优化结构 cd $dir || { warning 无法进入$dir跳过; continue; } # 结构优化循环强制用CONTCAR opt_donefalse; attempt0; max_attempts3 while [ $attempt -lt $max_attempts ] [ $opt_done false ]; do attempt$((attempt 1)) log 第$attempt次优化$p GPa结构... log_filevasp_${p}GPa_$attempt.log if ! $vasp_cmd $log_file 21; then warning $p GPa第$attempt次VASP运行失败日志$log_file continue fi # 核心检查CONTCAR是否有效不管收敛与否 if check_non_empty_file CONTCAR; then cp CONTCAR CONTCAR.${p}GPa # 备份当前压力点的CONTCAR current_contcar../$dir/CONTCAR.${p}GPa # 更新下一个压力点的输入 log ✅ $p GPa生成有效CONTCAR强制使用该文件作为后续输入 opt_donetrue # 可选检查收敛条件仅日志提示 ionic_steps$(grep -E ^ [0-9] F $log_file | wc -l) reached$(grep -c reached required accuracy $log_file) if [ $ionic_steps -lt 3 ] [ $reached -ge 1 ]; then log ✅ $p GPa优化收敛离子步数$ionic_steps else warning $p GPa优化未收敛但仍使用CONTCAR作为后续输入 fi # 检查空间群是否变化仅能带计算时 if [ $do_band y ] [ $kpath_updated false ]; then current_sg$(get_space_group CONTCAR) if [ $current_sg ! $initial_space_group ] [ $current_sg ! UNKNOWN ]; then warning 空间群变化初始$initial_space_group → 当前$current_sg read -p 是否更新KPATH.in(y/n默认n) update_kpath; update_kpath${update_kpath:-n} if [ $update_kpath y ]; then # 基于当前CONTCAR生成新的KPATH.in vaspkit -task 303 cp KPATH.in ../../KPATH.in kpath_updatedtrue log ✅ 已基于当前结构更新KPATH.in fi fi fi else warning $p GPa第$attempt次优化未生成有效CONTCAR重试... continue fi done [ $opt_done false ] { warning $p GPa优化失败跳过; cd ..; continue; } # 能带计算原胞选项 if [ $do_band y ]; then log 开始$p GPa能带计算... mkdir -p band cd band || { warning 无法进入band跳过; cd ..; continue; } # 复制当前压力点的CONTCAR到band目录 cp ../CONTCAR POSCAR # 如果使用原胞先生成PRIMCELL.vasp并替换POSCAR再生成K点 if [ $use_prim_cell true ]; then log 生成原胞PRIMCELL.vaspvaspkit -task 602 vaspkit -task 602 if check_non_empty_file PRIMCELL.vasp; then mv POSCAR POSCAR_ori # 备份原结构 cp PRIMCELL.vasp POSCAR log ✅ 已替换为原胞PRIMCELL.vasp作为能带计算输入 # 重新生成K点基于原胞 log 基于原胞生成新的K点 (echo 102; echo 2; echo 0.04) | vaspkit else warning 生成原胞失败使用原CONTCAR继续 fi fi # 复制其他文件 cp ../../POTCAR ../../INCAR . # 处理能带INCAR if [ -n $band_incar ]; then cp $band_incar INCAR else sed -i /NSW/s/^/#/ INCAR; sed -i /IBRION/s/^/#/ INCAR sed -i /LWAVE/d INCAR; sed -i /LCHARG/d INCAR echo LWAVE .TRUE. INCAR; echo LCHARG .TRUE. INCAR echo LORBIT 10 INCAR; echo NEDOS 5001 INCAR fi # 第一步SCF计算 log 执行第一步SCF自洽计算 $vasp_cmd band_calc1.log 21 # 第二步能带计算替换KPATH.in cp ../../KPATH.in KPOINTS log 执行第二步能带计算 $vasp_cmd band_calc2.log 21 # 数据处理 log 处理能带数据 (echo 211; echo 1)|vaspkit (echo 215; echo 1)|vaspkit (echo 111; echo 2)|vaspkit (echo 211; echo 2)|vaspkit (echo 113; echo all;echo 2)|vaspkit (echo 213; echo 2;echo all)|vaspkit cd .. fi # 声子谱计算强制用CONTCAR if [ $do_phono y ]; then log 开始$p GPa声子谱计算... mkdir -p phono cd phono || { warning 无法进入phono跳过; cd ..; continue; } # 复制当前压力点的CONTCAR作为输入 cp ../CONTCAR POSCAR cp ../../POTCAR ../../KPOINTS ../../INCAR . # 处理声子INCAR if [ -n $phono_incar ]; then cp $phono_incar INCAR else sed -i /NCORE/s/^/#/ INCAR sed -i /NSW/d INCAR; sed -i /IBRION/d INCAR echo NSW 1 INCAR grep -q IVDW INCAR echo IBRION 6 INCAR || echo IBRION 8 INCAR fi # 运行VASPphonopy log 执行声子谱VASP计算 $vasp_cmd phono_vasp.log 21 if check_non_empty_file vasprun.xml; then # 生成KPATH.phonopy cat KPATH.phonopy EOF NPOINTS 501 DIM 1 1 1 BAND AUTO BAND_LABELS AUTO TETRAHEDRON .TRUE. PDOS AUTO BAND_CONNECTION .TRUE. FORCE_CONSTANTS READ EOF # 运行phonopy log 执行phonopy计算 phonopy --fc vasprun.xml phonopy -s -p KPATH.phonopy --paauto phonopy-bandplot --gnuplot | tee phono.txt rm -f *.yaml log ✅ 声子谱计算完成 else warning 未生成有效vasprun.xml跳过phonopy计算 fi cd .. fi cd .. # 返回主目录 done # 计算完成 log ${GREEN} 所有计算完成 ${NC} log ✅ 初始CONTCAR$current_contcar跳过常压优化/ CONTCAR.0GPa常压优化 log ✅ 压力点CONTCARP*GPa/CONTCAR.*GPa每个压力点的核心输出 log ✅ 关键文件 log - 能带计算P*GPa/band/使用CONTCAR/原胞PRIMCELL.vasp作为输入 log - 声子谱计算P*GPa/phono/使用CONTCAR作为输入 log ✅ 注意所有后续计算均强制使用CONTCAR无论是否收敛
