利用VASP进行弹性模量的计算笔记

本文是基于《基于能量-应变关系计算弹性常数—VASPKIT v1.00新功能》、《基于应力-应变关系计算弹性常数—VASPKIT v1.2.0新功能》，并在此基础上进行整理记录，在此对原分享者表示感谢。

1 准备优化彻底的POSCAR文件，注意通常采用标准的惯用原胞计算弹性常数，如果不确信POSCAR文件中是否是标准的惯用原胞，可以用vaspkit-603/604生成标准结构；

补充：在Materials Project中下载所需材料的CIF文件，利用VESTA导出为XX.VASP文件，重命名为POSCAR

2 结构优化（弛豫）

INCAR 获取：vaspkit → 101 → LR

KPOINTS 获取： vaspkit →102→ 2→ 0.03 （由于计算弹性常数对K-mesh要求很高，因此对于半导体（金属体）体系，生成K点的精度应不小于0.03(0.02) * 2π Å-1）

2.1 如果弛豫发现报错

最好重新走一遍结构优化，或者INCAR中加上【ISYM=0】忽略对称性（但是增加计算量）

3 弹性模量计算文件准备

POSCAR : 使用结构优化后的新坐标， cp CONTCAR POSCAR

KPOINTS : vaspkit →102 →2 →0.03

INCAR : vaspkit →101 →DC , 可以除了PREC参数，也可以手动设置ENCUT=1.5*ENMAX

Global Parameters
ISTART = 0
LREAL = F
PREC   = High （截断能设置默认值1.5-2倍）
LWAVE = F
LCHARG = F
ADDGRID= .TRUE.
Electronic Relaxation
ISMEAR = 0
SIGMA = 0.05
NELM   = 40
NELMIN = 4
EDIFF = 1E-08
Ionic Relaxation
NELMIN = 6
NSW   = 100
IBRION = 2
ISIF   = 2   （切记选择2,如果选择3会把施加应变后原胞重新优化成平衡原胞）
EDIFFG = -1E-02

4 准备VPKIT.in文件，内容如下

1                   ! 设置1将产生计算弹性常数的输入文件，2则计算弹性常数
3D                   ! 2D为二维体系，3D为三维体系
7                   ! 7个应变
-0.015 -0.010 -0.005 0.000 0.005 0.010 0.015 ! 应变变化范围

运行vaspkit 200（或201），分别对应【用应力-应变法计算弹性常数】、【能量-应变法计算弹性常数】

根据POSCAR产生对应张量元文件夹和计算文件

5 准备批量提交vasp作业的脚本（txml.sh）

#这条代码运行后，会在每个张量元的应变计算后，把超算脚本vasp.job复制到下一文件夹，继续张量元的计算
#!/bin/bash
root_path=`pwd`
for cij in `ls -F | grep /$`
do
cd ${root_path}/$cij
for s in strain_*
do
cd ${root_path}/$cij/$s
echo `pwd`cp ../../vasp.job .            ./vasp.job
# 上面两行vasp.job 替换为 超算运算的脚本文件名
# Add here your vasp_submit_job_script     把这一行替换为提交运算的命令
done
done

上述脚本要授予权限，否则如下图

利用chmod命令：chmod 777 txml.sh
6 计算结束后，再次修改VPKIT.in文件中第一行为2（后处理），然后再次运行vaspkit并选择200/201，即可得到弹性模量矩阵

利用VASP进行弹性模量的计算笔记相关推荐

板式橡胶支座弹性模量怎样计算_Midas Civil支座模拟的详细解决方法
MIDAS支座模拟误区误区1.在支座位置建立节点,并将所有的支座节点与主梁节点刚性连接,荷载会按照就近原则分配,导致离主梁越近的支座反力越大. 误区2.当用弹性连接的一般类型模拟单支座时,如果没有定 ...
python基础教程：Python中利用sqrt()方法进行平方根计算的教程
这篇文章主要介绍了Python中利用sqrt()方法进行平方根计算的教程,是Python学习的基础知识,需要的朋友可以参考下 sqrt()方法返回x的平方根(x>0). 语法以下是sqrt() ...
【深度学习】利用深度可分离卷积减小计算量及提升网络性能
[深度学习]利用深度可分离卷积减小计算量及提升网络性能文章目录 1 深度可分离卷积 2 一个深度可分离卷积层的代码示例(keras) 3 优势与创新3.1 Depthwise 过程3.2 Point ...
如何写python脚本抓取数据并计算_【小工具】利用Python脚本从Gaussian计算结果中提取信息...
1.前言高斯(Gaussian)是一个功能强大的量子化学综合软件包,所有从事计算化学相关领域的科研工作者应该都使用或者了解过这个软件.它的输出文件(.log文件)是一个文本文件,可以利用文本工具打开 ...
TCP/IP学习——利用主机IP与子网掩码计算网络号
阿里巴巴的一道题目,是关于利用主机IP与子网掩码计算网络号. 在此之前,我们需要先了解下五类IP地址. 在同一个物理网络上的主机的地址都有一个相同前缀,即IP地址分成两个部分:netid,hostid ...
[并发编程]并发编程第二篇：利用并发编程，实现计算大量数据的和
利用并发编程,实现计算大量数据的和实现代码: package tj.pojo.generate.main;import java.util.ArrayList; import java.util.L ...
利用blink CEP实现流计算中的超时统计问题
案例与解决方案汇总页: 阿里云实时计算产品案例&解决方案汇总一. 背景介绍如<利用blink+MQ实现流计算中的延时统计问题>一文中所描述的场景,我们将其简化为以下案例: 实时 ...
python利用tushare下载数据并计算当日收益率
python利用tushare下载数据并计算当日收益率计算股票收益率的程序主要有以下几部分构成: 1.获取股票接口数据函数:pro_daily_stock() 2.计算收益率函数:cal_stock ...
利用球谐系数计算函数值及利用EGM球谐系数计算重力异常
利用球谐系数计算函数值及利用EGM球谐系数计算重力异常参考文章: (1)利用球谐系数计算函数值及利用EGM球谐系数计算重力异常 (2)https://www.cnblogs.com/heaventi ...

利用VASP进行弹性模量的计算笔记

利用VASP进行弹性模量的计算笔记相关推荐

最新文章

热门文章