24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2769  |  回复: 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的回帖

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的回帖
查看全部 13 个回答

zzy870720z

荣誉版主 (文坛精英)

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

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

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿郑州大学085600求调剂 +19 吃的不少 2026-04-05 22/1100 2026-04-06 22:39 by tianyyysss
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +18 披星河 2026-04-03 18/900 2026-04-06 13:55 by BruceLiu320
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 男生,一志愿沪9生物学071000,初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 专硕310求调剂 +5 捞捞我…. 2026-04-04 6/300 2026-04-04 23:33 by barlinike
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 求调剂:085600材料与化工,考材科基,总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 0856调剂 +8 曲听筠 2026-03-30 8/400 2026-04-04 08:46 by tianyyysss
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 一志愿重庆大学085404,总分314分,求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 348求调剂 +6 吴彦祖24k 2026-04-02 6/300 2026-04-02 14:07 by 给你你注意休息
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 372求调剂 +3 jj涌77 2026-04-02 3/150 2026-04-02 09:57 by olim
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
信息提示
请填处理意见