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新功能

相关推荐
学习路上_write12 小时前
新版Pycharm添加导入anaconda的python解释器
开发语言·python·pycharm
光影少年13 小时前
AI大模型开发语言排行
开发语言·人工智能
爱和冰阔落13 小时前
【C++list】底层结构、迭代器核心原理与常用接口实现全解析
开发语言·数据结构·c++·list
编程岁月13 小时前
java面试0106-java什么时候会出现i>i+1和i<i-1?
java·开发语言·面试
TH888613 小时前
小麦赤霉病监测设备:通过多维度数据采集与智能分析,实现病害的早发现、早预警、早防控
python
练习时长一年14 小时前
Java开发者进阶之路
java·开发语言
去往火星14 小时前
文字转语音——sherpa-onnx语音识别离线部署C++实现
开发语言·c++
Q_Q196328847515 小时前
python+springboot+uniapp基于微信小程序的校园二手闲置二手交易公益系统 二手交易+公益捐赠
spring boot·python·django·flask·uni-app·node.js·php
semicolon_hello15 小时前
C++中 optional variant any 的使用
开发语言·c++
报错小能手16 小时前
C++笔记(基础)string基础
开发语言·c++·笔记