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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】CASTEP做Dynamic时晶体易变……
51/1返回列表
查看: 2701  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

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

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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
思想重于技巧，内涵重于表象
1楼 2010-04-28 12:11:40
已阅   回复此楼   关注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试试吧
赞一下(2人) 回复此楼
博学、审问、慎思、明辨、笃学
4楼2010-04-28 12:37:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工304求B区调剂 +6 邱gl 2026-03-11 7/350 2026-03-17 01:27 by 学员FnSWZj
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 0703化学调剂，求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 304求调剂 +6 小熊joy 2026-03-14 6/300 2026-03-16 12:59 by Iveryant
[考研] 311求调剂 +6 冬十三 2026-03-15 6/300 2026-03-16 08:00 by wang_dand
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
关闭
请填处理意见
关闭