vasp计算弹性常数

  1. 准备优化好的POSCAR文件

  2. 在POSCAR目录下运行下面代码,得到KPOINTS和INCAR:

bash 复制代码
vaspkit → 102 → 2 → 0.03

3.拖入VPKIT.in,注意第一行写1,用于前处理

bash 复制代码
1                     ! 1 for pre-processing; 2 for post-processing
3D                    ! 2D for slab, 3D for bulk                                       
11                    ! number of strain case
-0.010 -0.007 -0.005 -0.003 -0.001  0.000 0.001 0.003 0.005 0.007 0.010   ! Strain range

在当前目录下运行:

bash 复制代码
vaspkit → 200

可以看到类似下面的输出:

  1. 生成的文件夹下批量提交vasp作业即可,我的脚本如下:
python 复制代码
import os

def batch_submit_strain_jobs():
    for folder in sorted(os.listdir()):
        if folder.startswith("strain_") and os.path.isdir(folder):
            print(f"Submitting VASP job in: {folder}")
            os.chdir(folder)
            os.system("runvasp -q short")  # 根据集群需求可改为 vasp_std > out 或 sbatch job.sh
            os.chdir("..")
    print("✅ All strain jobs submitted.")

batch_submit_strain_jobs()
  1. 等VASP全部计算完成之后,再次修改VPKIT.in文件中第一行为2 (后处理),然后再次运行VASPKIT并选择200,得到以下结果

参考链接:基于应力-应变关系计算弹性常数---VASPKIT v1.2.0新功能

相关推荐
极客小张13 小时前
【项目思路】基于STM32+ZigBee的智能家居--浴室场景设计
c语言·python·stm32·智能家居·课程设计·项目设计·企业项目
墩墩分墩15 小时前
【Go语言入门教程】 Go语言的起源与技术特点:从诞生到现代编程利器(一)
开发语言·后端·golang·go
CHANG_THE_WORLD15 小时前
并发编程指南 同步操作与强制排序
开发语言·c++·算法
仰泳之鹅16 小时前
【C语言】深入理解指针(5)
c语言·开发语言
无为之士17 小时前
君正交叉编译链工具mips-gcc540-glibc222-64bit-r3.3.0.smaller.bz2编译st-device-sdk-c
c语言·开发语言
源力祁老师18 小时前
深入分析 json2(新)与标准的 jsonrpc的区别
开发语言
小wanga18 小时前
C++知识
java·开发语言·c++
学渣6765618 小时前
文件传输工具rsync|rust开发环境安装|Ascend实验相关命令
开发语言·后端·rust
木心爱编程18 小时前
C++容器内存布局与性能优化指南
开发语言·c++·性能优化