24小时热门版块排行榜    

查看: 3282  |  回复: 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的回帖

franch

版主 (著名写手)

souledge(金币+5): 剩余的钱没地儿解决了,哈哈~ 2011-04-26 08:56:36
引用回帖:
Originally posted by souledge at 2010-04-28 17:17:55:
咋没人回帖了捏……

问题解决没?
遇弱则弱
10楼2011-04-03 19:38:08
已阅   回复此楼   关注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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:45 by KM5EcsNQRBPn
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +3 process2012 2026-08-23 4/200 2026-08-23 12:54 by huixian257
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
信息提示
请填处理意见