24小时热门版块排行榜    

查看: 2690  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

souledge

专家顾问 (著名写手)

[交流] 【求助】CASTEP做Dynamic时晶体易变…… 已有4人参与

本人课题与热力学等问题联系比较多,所以想用CASTEP中的Dynamic来求得不同晶体的热力学数据(关键是吉布斯函数)。
后来看到了小木虫计算版期刊第一期中wenhua luo和lei0736的《自由能的分子动力学计算》一文,开始按照里面的步骤学习使用Dynamic来模拟计算。
在内能模拟时,先从MS的数据库中直接导入了Ag,然后扩展到8x8x8的结构,再Make P1。然后按照文章中的步骤,选择NPT系综,每个时间步为2fs,温度为300K,然后开始计算。结果总是在进行几个时间步后就提示晶体易变,任务失败。
为了加快全流程速度,改为设置了10fs为一个时间步,共100时间步,Thermostat使用Nose,Barostat使用Andersen,并且为固定了质心(Fix center of mass)。并且为了增强收敛性,加大了截断能,并且使用了8x8x8的K点设置。采用LDA超软赝势进行计算,结果依然有错误提示:
引用回帖:
Error in md_vv1_npt_nose - cell inconsistency?
Current trace stack:
md_vv1_npt_nose
md_run
castep

顺便看了一下结果图,的确这个突变有点太那啥了……

由于是工科背景,对于赝势之类的理论比较空白,虽然在恶补,但是暂时依然杯水车薪。哪位高人可以帮忙指出一下,上面的那些步骤,都有哪些问题么……
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

思想重于技巧,内涵重于表象
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主

★ ★
souledge(金币+1):补加一个,呵呵 2010-04-28 13:02
hometownlove(金币+2):谢谢交流 2010-04-28 18:38
不是什么高人,呵呵
提高k点能提高精度,不过进行计算前最好测一下k点和cutoff,选取合适的k点和cutoff就能保证精度,没有必要直接就设那么高的k点
用1fs试试吧
博学、审问、慎思、明辨、笃学
4楼2010-04-28 12:37:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主

souledge(金币+1):感谢回复,不过还是不明白,能详细说说么? 2010-04-28 12:33
时间步长缩小试试,为增加速度,减小点k点
博学、审问、慎思、明辨、笃学
2楼2010-04-28 12:23:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

引用回帖:
Originally posted by zzy870720z at 2010-04-28 12:23:31:
时间步长缩小试试,为增加速度,减小点k点

可是之前的2fs也一样出这个问题啊,用1fs试试?
另外,记得增加K点会提高计算精度和增加收敛,但是会造成过大的计算负担,可是为什么要减少K点呢?
现在的最重要问题就是总出现晶格易变,无法继续计算。速度只是为了能够快点看到是否出现相同问题……算是try error了……
刚刚接触这个的新人,还望高人详细指点迷津

[ Last edited by souledge on 2010-4-28 at 12:34 ]
思想重于技巧,内涵重于表象
3楼2010-04-28 12:33:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

引用回帖:
Originally posted by zzy870720z at 2010-04-28 12:37:26:
不是什么高人,呵呵
提高k点能提高精度,不过进行计算前最好测一下k点和cutoff,选取合适的k点和cutoff就能保证精度,没有必要直接就设那么高的k点
用1fs试试吧

把截断能降低到了300,把K点也降低到了3x3x3,时间步降低到了1fs,甚至整个的Ag的晶体阵也降低到了4x4x4,速度到是明显上去了,可是,晶体易变的错误依旧……
引用回帖:
Error in md_vv1_npt_nose - cell inconsistency?
Current trace stack:
md_vv1_npt_nose
md_run
castep
Error in md_vv1_npt_nose - cell inconsistency?
Current trace stack:
md_vv1_npt_nose
md_run
castep
MPI Application rank 3 exited before MPI_Finalize() with status 1

为了防止这个问题,还在此前特地几何优化了一下4x4x4的结构。不过,由于内存限制,无法把这个结构做成Supercell再计算,只能make P1。
再上两张计算出来的图,这个变化依然是那么强烈……

思想重于技巧,内涵重于表象
5楼2010-04-28 12:54:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 9/450 2026-03-16 17:35 by 沐霖12138
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
信息提示
请填处理意见