应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】一定压强下的弹性常数计算方法
262/3返回列表上一页123下一页
查看: 4758  |  回复: 25
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fengshiquan

铁杆木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by babaleo at 2011-04-01 07:29:37:
搜一下侯博士的指南,(100000) 是来说明应变(e1,e2,e3,e4,e5,e6)里的e1不为零, 其他五个都是零, e1取几个值算一下,就可以stresss=cij*strain 拟合了. 文献里看到别人有的就算e=-0.01 和0.01下的stress, 然后直线 ...

谢谢babaleo的回答，我还是不是很明白，100000) 用来说明应变(e1,e2,e3,e4,e5,e6)是在INCAR里设置的吗，怎么设置，
还有你说用stresss=cij*strain 拟合不明白是什么意思
侯博士的指南我也搜到两个，不过没有弹性常数拟合这块，只有采用Murnaghan状态方程拟合得到晶格常数和体积模量，你能不能给我发一份你搜到那份侯博士的指南，我邮箱49273185@163.com,谢谢虫友了
赞一下(1人) 回复此楼
11楼2011-04-01 10:00:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gzqdyouxia(金币+2): 谢谢提示 2011-04-01 12:01:57
引用回帖:
Originally posted by fengshiquan at 2011-04-01 10:00:22:
谢谢babaleo的回答，我还是不是很明白，100000) 用来说明应变(e1,e2,e3,e4,e5,e6)是在INCAR里设置的吗，怎么设置，
还有你说用stresss=cij*strain 拟合不明白是什么意思
侯博士的指南我也搜到两个，不过没有弹 ...

加应变(e1,e2,e3,e4,e5,e6)通过 POSCAR来表示, 就是侯博士指南里的”defvector.f”那个程序. 不加应变优化后的CONTCAR拷成OLDPOS, 加注原子数, 执行那个程序可以生成不同应变后的POSCAR, 然后ISIF=2优化后再静态算, 就得到能量和应力了. 这一步不管是用应变-能量法, 或者应力-应变法都是一样的.

比如Orthorombic, 有9个非零 Cij, 加应变(00e00e)的话,
        s1                C11        C12        C13        0        0        0                0
        s2                C12        C22        C23        0        0        0                0
        s3        =        C13        C23        C33        0        0        0        X        e
        s4                0        0        0        C44        0        0                0
        s5                0        0        0        0        C55        0                0
        s6                0        0        0        0        0        C66                e

对应力s1, 从上面矩阵s1=c13xe, 加不同应变, 算出不同应力, 比如,
strain magnitute1=-0.006,        Stress s1== -30.775258
strain magnitute2=-0.0036,         stress s1=-30.376274
strain magnitute3=-0.0012,         Stress s1=-29.925712
strain magnitute4=0.0012,        Stress s1=-29.504233
strain magnitute5=0.0036, Stress s1=-29.063
strain magnitute6=0.006, Stress s1=-28.6

对上面应力-应变值 拟和, 拟和直线斜率就是c13班(c13=181.4GPa).同理,用s2, s3, s6可得到 c23,c33, c66.
赞一下(2人) 回复此楼
12楼2011-04-01 11:35:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)

您好！请问：

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by babaleo at 2011-04-01 11:35:40:
加应变(e1,e2,e3,e4,e5,e6)通过 POSCAR来表示, 就是侯博士指南里的”defvector.f”那个程序. 不加应变优化后的CONTCAR拷成OLDPOS, 加注原子数, 执行那个程序可以生成不同应变后的POSCAR, 然后ISIF=2优化后再静态 ...

请问在上面您提到的“...然后ISIF=2优化后再静态算, 就得到能量和应力了. 这一步不管是用应变-能量法, 或者应力-应变法都是一样的”中，如果用的是应力-应变法，优化后的静态计算的INCAR中，ISIF=?,要不要IBRION?还要加多哪些参数才能得到应力？（我用过应变-能量法）。

      盼回复，谢谢！
赞一下(1人) 回复此楼
好好学习
13楼2011-04-15 20:39:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)

fengshiquan，我这有一个你想要的文件

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1): 谢谢分享 2011-04-16 00:17:57
引用回帖:
Originally posted by fengshiquan at 2011-04-01 10:00:22:
谢谢babaleo的回答，我还是不是很明白，100000) 用来说明应变(e1,e2,e3,e4,e5,e6)是在INCAR里设置的吗，怎么设置，
还有你说用stresss=cij*strain 拟合不明白是什么意思
侯博士的指南我也搜到两个，不过没有弹 ...

我这有一个你想要的文件：就是侯老师那个
赞一下(2人) 回复此楼
好好学习
14楼2011-04-15 20:43:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by donggaomin at 2011-04-15 20:43:31:
我这有一个你想要的文件：就是侯老师那个

呵呵呵，谢谢，我已经找到了
赞一下(1人) 回复此楼
15楼2011-04-15 21:12:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)

请问babaleo

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by babaleo at 2011-04-01 11:35:40:
加应变(e1,e2,e3,e4,e5,e6)通过 POSCAR来表示, 就是侯博士指南里的”defvector.f”那个程序. 不加应变优化后的CONTCAR拷成OLDPOS, 加注原子数, 执行那个程序可以生成不同应变后的POSCAR, 然后ISIF=2优化后再静态 ...

请问babaleo，上面所说的 stress s1    是不是您前面给出的F1？（grep "in kB' OUTCAR,得F1，F2，F3，F4，F5，F6），stress s6就是F6吗？先谢谢啦！
赞一下(1人) 回复此楼
好好学习
16楼2011-04-15 23:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢指教 2011-04-16 08:32:55
引用回帖:
Originally posted by donggaomin at 2011-04-15 23:06:00:
请问babaleo，上面所说的 stress s1    是不是您前面给出的F1？（grep "in kB' OUTCAR,得F1，F2，F3，F4，F5，F6），stress s6就是F6吗？先谢谢啦！

是的。 只不过在vasp4.6版里，输出应力顺序是"

" Direction    X        Y        Z        XY       YZ       ZX"
这样一来，就变成了 “F1, F2, F3, F6, F4, F5".

[ Last edited by babaleo on 2011-4-16 at 06:44 ]
赞一下(2人) 回复此楼
17楼2011-04-16 03:37:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by babaleo at 2011-04-01 07:29:37:
搜一下侯博士的指南,(100000) 是来说明应变(e1,e2,e3,e4,e5,e6)里的e1不为零, 其他五个都是零, e1取几个值算一下,就可以stresss=cij*strain 拟合了. 文献里看到别人有的就算e=-0.01 和0.01下的stress, 然后直线 ...

请问虫友，
（1）计算弹性常数的程序defvector.f的输入文件是不是未加应力前优化晶格参数得到的CONCAR文件
（2）编译defvector.f后得到模块defvector.x,如何再由他得到应变后的初始POSCAR文件了，
谢谢虫友了，还望虫友不吝赐教
赞一下(1人) 回复此楼
18楼2011-04-26 20:33:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by fengshiquan at 2011-04-26 20:33:58:
请问虫友，
（1）计算弹性常数的程序defvector.f的输入文件是不是未加应力前优化晶格参数得到的CONCAR文件
（2）编译defvector.f后得到模块defvector.x,如何再由他得到应变后的初始POSCAR文件了，
谢谢虫友 ...

（1）输入文件是在你想做的压力下所优化晶格参数得到的CONCAR文件烤成名为OLDPOS并在第一行上加注总原子数（参考侯博士指南）
（2）执行模块defvector.x后便生成一个fort.3文件，然后铐成POSCAR就行了。侯博士指南里已经把这步脚本化进去了。
另外，坛内虫友dongguomin刚刚发了一个贴， http://muchong.com/bbs/viewthread.php?tid=3123071&fpage=1 步骤写得很详细。
赞一下(1人) 回复此楼
19楼2011-04-26 21:21:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by babaleo at 2011-04-26 21:21:15:
（1）输入文件是在你想做的压力下所优化晶格参数得到的CONCAR文件烤成名为OLDPOS并在第一行上加注总原子数（参考侯博士指南）
（2）执行模块defvector.x后便生成一个fort.3文件，然后铐成POSCAR就行了。侯博士 ...

谢谢赐教，去看一下
赞一下(1人) 回复此楼
20楼2011-04-27 11:13:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xh512 的主题更新
262/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭