使用vasp计算α-SiO2弹性模量(个人记录)
使用vasp计算α-SiO2弹性模量
文件下载
由于vasp wiki的example里有alpha-SiO2的样例,故INCAR、KPOINTS、POSCAR、POTCAR文件可以从vasp wiki直接下载
样例链接:https://www.vasp.at/wiki/index.php/Alpha-SiO2
输入文件下载链接:https://www.vasp.at/wiki/images/4/4e/SiO2_NMR.tgz
输入文件选用
因为需要计算弹性模量,所以修改(添加)如下INCAR文件参数
IBRION = 6
ISIF = 3
NFREE = 4
NSW = 1
ENCUT收敛性测试
选用官网样例中默认的KPOINTS.4x4x3作为KPOINTS文件,测试ENCUT收敛性
4x4x30
G4 4 30 0 0
修改INCAR中的ENCUT参数为AAA并另存为INCAR0,编写bash如下,测试ENCUT收敛性, 将不同ENCUT值时的能量写入energy文件
#!/bin/bash
for i in 200 250 300 350 400 450 500 550 600 # 候选ENCUT值dosed "s/AAA/$i/g" INCAR0 > INCAR # 文本替换mkdir $icp vasp_std INCAR POTCAR POSCAR KPOINTS $i/cd $iulimit -s unlimited # 关闭堆栈大小限制mpirun -np 8 vasp_std 2>&1 | tee log # vasp主程序运行echo $icd ../doneecho "#encut Energy" > energy
for i in 200 250 300 350 400 450 500 550 600doE=`grep -A 2 "free\ energy" $i/OUTCAR | tail -1 | cut -c33-44` # 字符串处理,读取能量值echo $i " " $E >> energydone
计算得到energy文件如下
#encut Energy
200 -69.62312758
250 -71.76037373
300 -70.73853595
350 -69.96559936
400 -69.75761344
450 -69.65430805
500 -69.61260083
550 -69.62444570
600 -69.64168679
规律不明显,故选取ENCUT=400作为截断能(和官网样例相同)
KPOINTS收敛性测试
如上一步所示,选取ENCUT=400作为截断能
修改KPOINTS中的晶格参数为AAA并另存为KPOINTS0,编写bash如下,测试KPOINTS收敛性, 将不同KPOINTS值时的能量写入energy文件
#!/bin/bash
for i in "3 3 2" "4 4 3" "5 5 4" # 几个候选KPOINTS值dosed "s/AAA/$i/g" KPOINTS0 > KPOINTS # 文本替换mkdir "$i"cp vasp_std INCAR POTCAR POSCAR KPOINTS "$i"/cd "$i"ulimit -s unlimited # 关闭堆栈大小限制mpirun -np 8 ./vasp_std 2>&1 | tee log # 主程序运行echo $icd ../doneecho "#k-points Energy" >energy
for i in "3 3 2" "4 4 3" "5 5 4"doE=`grep -A 2 "free\ energy" "$i"/OUTCAR | tail -1 | cut -c33-44` # 字符串处理,读取能量值echo $i " " $E >> energydone
计算得到energy文件如下
#k-points Energy
3 3 2 -69.75780682
4 4 3 -69.75761344
5 5 4 -69.76029387
规律不明显,故选取4*4*3作为KPOINTS(也和官网样例相同)
结构弛豫计算 优化晶格
修改INCAR文件中的以下参数如下,准备做结构弛豫计算
IBRION = 6
NFREE = 4
ISIF = 3
NSF = 100
ENCUT = 520 # 400*1.3
计算后使用vimdiff对比POSCAR、CONTCAR两个文件的差异如下
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-hKOnaCPG-1665304896101)(/home/ian/.config/Typora/typora-user-images/image-20221009160546639.png)]
但由于结构优化最终并未收敛,故并未采取CONTCAR优化过的晶格替代初始POSCAR进行最终计算
最终计算 计算弹性模量
修改INCAR文件中的以下参数如下,进行力学性能计算,计算弹性模量
IBRION = 6
NFREE = 4
ISIF = 3
NSW = 1
ENCUT = 400
计算完成后可在OUTCAR文件中找到弹性矩阵如下:
TOTAL ELASTIC MODULI (kBar)Direction XX YY ZZ XY YZ ZX--------------------------------------------------------------------------------XX 1866.1602 729.1319 756.9221 -0.0000 44.6970 -0.0000YY 729.1319 1866.1602 756.9221 0.0000 -44.6970 0.0000ZZ 756.9221 756.9221 2800.3034 -0.0000 0.0000 0.0000XY -0.0000 0.0000 0.0000 568.5142 -0.0000 44.6970YZ 44.6970 -44.6970 0.0000 -0.0000 835.5214 -0.0000ZX -0.0000 0.0000 0.0000 44.6970 -0.0000 835.5214--------------------------------------------------------------------------------
之后便可以使用力学性质在线作图工具(链接:http://progs.coudert.name/elate)读取数据并作图
使用vasp计算α-SiO2弹性模量(个人记录)相关推荐
- VASP计算新手礼包:基于CentOS 7.6 系统安装VASP 5.4.1
2020-9-24更新:更新VASP5.4.1网盘链接: 链接:https://pan.baidu.com/s/1n9C9CmF0j7iMilo1cPe1cA 提取码:69dj 2019-9-9 更新 ...
- 使用VASP计算材料能带和态密度
前言 一.前期准备 二.计算流程 1.静态自洽(scf) 2.能带计算 3.态密度计算 总结 前言 能带结构是目前采用第一性原理计算所得到的常用信息,可用来结合解释金属.半导体和绝缘体的区别.能带用来 ...
- VASP计算光学性质
PS:请见文末的打赏选项. VASP计算光学性质的INCAR是根据HSE06计算的自洽INCAR来改的 SCF INCAR: SYSTEM = System ISTART = 1 ICHARG = 1 ...
- 分页总数计算,根据总记录数获取分页总数_ZZGoals
分页总数计算,根据总记录数获取分页总数,本文章提供5种写法,大家感觉哪个好用就拿去用吧. int pageSize = 20; //每页记录数量 int totalCount = 100005; // ...
- VASP计算笔记_声子谱计算
声子谱计算 主要使用phonopy软件计算声子谱.版本:phonopy-2.1.3 密度泛函微软理论/线性响应方法 (DFPT) 必要的输入文件: INCAR KPOINTS POSCAR-unitc ...
- VASP计算AgGaS2能带及态密度及光学性质
第一步:构型优化 1.准备四个输入文件 POSCAR INCAR POTCAR KPOINTS POSCAR: 从ms中导入AgGaS2结构,选择CASTEP,file,save,并保存成原包.这样, ...
- VASP计算弹性常数
记录一下最近计算的不同压力下Cr3C2的弹性常数 1. vasp直接计算 vasp5.1版本后支持设置IBRION=6来直接计算材料的弹性常数,最直接的做法就是在INCAR中设置以下三个参数 IBRI ...
- vasp计算之输入输出文件
一.vasp文件: INCAR in ** STOPCAR in stout out POTCAR in ** KPOINTS in ** IBZKPT out POSCAR in ** CONTCA ...
- 自定义起始时间的时间戳计算(个人代码记录)
个人代码记录,这个时间戳计算,由于项目数据特殊需求,自己设置的时间戳起始时间为2019年1月1日开始, 然后根据时间戳数据计算出现在时间,具体使用如下,这里setFullYear的月份是0到11,所以 ...
最新文章
- 安全也要“易”,谈NAC的硬件化
- 【数字逻辑入门】计算机如何存储1位二进制数
- C# 利用反射动态创建对象[摘录]
- Java时间类之间的转换案例——计算自己活了多少天
- 【Java并发编程】—–“J.U.C”:ArrayBlockingQueue
- python4.2_python4.2参数传入
- 登陆服务器老出现“达到最大连接数解决方法
- photoshop 新建文档尺寸预设如何导出保存
- ***X档案2007配套光盘
- 智慧城市综合管控平台
- winrm java客户端_Ansible for Windows:WinRM HTTPS设置
- R语言中的rgl包的安装问题
- 【MLDL】logistics regression理解
- dataBinding 配置问题 Could not find method buildFeatures()for arguments xxx on extension ‘android
- python pie图
- 时间触发以太网TTE
- 论ArcGIS10.2的Band Collection Statistics工具计算相关系数的正确性
- 高效学英语 - 统计英文书词频
- PCIE高性能RAID固态存储卡
- 十分钟带你了解阿里、美团、滴滴、头条等互联网头部大厂面经