📍

选择城市

立即咨询
精选文章 / 深度解读

bash+vasp+vaspkit能量应变计算弹性常数

脚本分享计算理论弹性常数常用的方法是应力应变,而另一种方法是根据能量变化,对平衡晶格结构施加小应变。能量-应变法对应于弹性刚度张量由总能量对应变的二阶导数导出。一般来说,应力应变法需要更高的计算精度才 ......

脚本分享


计算理论弹性常数常用的方法是应力应变,而另一种方法是根据能量变化,对平衡晶格结构施加小应变。能量-应变法对应于弹性刚度张量由总能量对应变的二阶导数导出。一般来说,应力应变法需要更高的计算精度才能达到与能量应变法相同的精度。然而,能量应变法需要的变形比后者小得多。能量应变法比应力应变法具有更小的应力敏感性,而且该方法已经在VASPKIT程序中实现。

独立弹性常数的数目取决于晶体的对称性。对称性越低意味着独立的弹性常数越多。例如,立方晶体有三个弹性常数,而三斜晶体有21个独立的弹性常数。



https://vaspkit.com/tutorials.html#mechanical-properties

而通过VASPKIT使用能量应变计算理论弹性常数时,生成的计算文件目录众多,需要辅助以bash脚本进行计算。


计算所需要的INCAR(参照vaspkit网站例子

Global Parameters
ISTART = 0
LREAL = F
PREC = High
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 (Must be 2)
EDIFFG = -1E-02

预处理和后处理所需要准备的http://VPKIT.in文件,如同后面注释所说明的,预处理时第一行开头为1,后处理时为2.

  • 1 ! 1 for pre-processing; 2 for post-processing
    3D ! 2D for two-dimentional, 3D for bulk
    7 ! number of strain case
    -0.015 -0.010 -0.005 0.000 0.005 0.010 0.015 ! Strain range

笔者准备了一个空间群为Pmm2(序号为25,属于斜方晶系)的POSCAR。



通过vaspkit的201功能生成了能量应变计算的目录




总计有9个大文件夹,分别以该空间群所需要的独立弹性常数来命名



在每个大文件夹内则是以应变强度命名的实际计算执行的文件夹,文件夹的数量可根据预处理的http://VPKIT.in文件的最后一行的应变设置来更改。




在每一个应变文件夹内则是施加了应变的POSCAR和事先准备好的INCAR、KPOINTS和POTCAR文件的链接,可直接进行vasp计算。

这里通过一个简单的两重for循环命令依次序完成所需要的计算任务。

for i in C11 C11_C12_C22 C11_C13_C33 C22 C22_C23_C13 C33 C44 C55 C66
#i为由独立弹性常数所命名的上级文件夹名称
do
cd $i
for k in strain_0.000 strain_-0.005 strain_+0.005 strain_-0.010 strain_+0.010 strain_-0.015 strain_+0.015
#k 为单独应变命名的文件的名称,可根据需要灵活改写
do
cd $k
mpirun -np 24 vasp|tee runlog
#vasp的执行命令,请根据实际计算环境和需要进行改写
cd ..
echo $i-$k
echo $i-$k
done
cd ..
done

待脚本执行完毕,可更改http://VPKIT.in文件进行后处理并导出。

为了可以将vaspkit所导出的关于声速、德拜温度等结果完整记录,建议执行如下命令

vaspkit -task 201|tee vpkitlog

运行过程记录则会保存在vpkitlog文件(可根据需要灵活命名)中。


愿有所成


更多请关注微信公众号:IEchoQ

引喻失义 妄自菲薄

阅读重点

这篇内容能帮你快速理解什么

帮助快速理解主题

通过更完整的主题说明和结构表达,帮助用户更快抓住重点,也让搜索系统更容易识别页面主题。

帮助判断下一步动作

让访问者快速理解当前问题、可行方法以及下一步应该继续看案例、看服务还是直接沟通。

帮助继续浏览与沟通

文章页不只是获取流量,也承担继续阅读、查看服务和发起咨询的承接作用。

继续了解

继续了解这个主题前,你可能还关心这些问题

为什么这类主题适合写成文章?

因为很多用户会通过问题词、对比词和方案词进入网站,文章页越清楚,越容易覆盖更具体的需求。

为什么文章页不能只有正文?

仅有正文不利于继续浏览和转化,文章页还需要总结、问答、相关推荐与咨询入口来承接用户。

看完之后下一步可以做什么?

可以继续看同类文章、服务页与案例页,也可以直接沟通官网升级与搜索优化需求。

这篇文章能帮助我解决什么具体问题?

这篇文章围绕当前主题提供了详细的解决方案、操作步骤和注意事项,帮助你快速理解核心要点并应用到实际场景中。

如何判断这篇文章的内容是否权威可靠?

内容基于实际项目经验和技术实践编写,结合行业标准和最佳实践,同时提供案例数据和方法论支撑,确保专业性和可操作性。

这类内容对SEO和网站排名有什么帮助?

优质的长文内容和FAQ结构能够提升页面主题相关性、增加用户停留时间、降低跳出率,这些都有助于搜索引擎评估页面质量并提升排名表现。

AI搜索引擎会如何理解和引用这类内容?

AI搜索系统会提取文章的实体信息、观点结论和结构化问答,当用户提出相关问题时,可能会引用本文作为答案来源或参考依据。

如果我有更多相关问题可以咨询谁?

可以通过页面底部的联系方式直接咨询我们的专业团队,包括电话、QQ或在线表单,我们会根据你的具体情况提供针对性的建议和方案。

这篇文章和同类内容有什么不同之处?

本文不仅提供理论知识,还包含实战经验、避坑指南和可执行的行动建议,同时兼顾传统SEO和新兴的GEO生成式搜索优化视角。

多久需要更新一次这类内容以保持时效性?

建议每季度审查并更新一次关键数据和案例,如果涉及技术工具或算法变化则需要更频繁地维护,确保内容持续为用户提供准确价值。