24小时热门版块排行榜    

查看: 2685  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +9 大米饭! 2026-03-15 9/450 2026-03-16 16:41 by 我的船我的海
[文学芳草园] 伙伴们,祝我生日快乐吧 +16 myrtle 2026-03-10 25/1250 2026-03-16 16:21 by 火星超人xi
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 304求调剂(085602一志愿985) +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
信息提示
请填处理意见