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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】一定压强下的弹性常数计算方法
262/3返回列表上一页123下一页
查看: 4887  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 07化学280分求调剂 +6 722865 2026-03-23 6/300 2026-03-24 20:58 by allen-yin
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-24 4/200 2026-03-24 16:41 by peike
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 289材料与化工（085600）B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 一志愿南昌大学，327分，材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 0856调剂，是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
信息提示
关闭
请填处理意见
关闭