24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 584  |  回复: 0
当前主题已经存档。

xillo

新虫 (初入文坛)

[交流] 【求助】【请教】如何在ms和vasp中固定单一应变

为了做应变应力曲线,看文献中的做法是要求对晶体施加一个单一的应变,并且固定住这个应变,优化其他的应变和晶格参数(也就是驰豫这一步)。

我目前的做法是在castep中手动修改晶格参数(立方晶格,或者重新定义基矢成互相垂直):
1、对于单轴拉伸,只需要将某一个基矢比如a拉长,然后固定住a的长度和晶轴夹角。然后优化另外两个基矢b,c,即可得到拉伸应力。这个应力就只会出现在stress张量的11分量。

2、对于单剪切,比如进行13方向的剪切。我目前是对c基矢手动增加一个a方向的矢量,比如c'=c+0.02a,其中每一个都是矢量。所以在设置上就是a和c的夹角变小,c轴变长。
但是这样就出现了一个问题,无法像单轴拉伸一样固定住这个剪切应变(shear strain),只有固定住ac的夹角以及c轴的长度,但是这样一来计算出的stress tensor就不只是13分量有值,而且33分量也有。不符合文献的要求(需要保证所有其他应力分量都为零,出了应变的13分量)

就此想请教朋友们:在CASTEP中或者VASP&MedeA中,用什么办法可以fix 单一剪切应变,而做晶格优化的时候可以使别的应变分量变化,从而得到合理的stress tensor。
谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xillo 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见