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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】CASTEP做Dynamic时晶体易变……
51/1返回列表
查看: 2688  |  回复: 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的回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主
这个还没遇到过，把你的设置和xsd文件发我邮箱，我也测一下，呵呵
zzy870720z@foxmail.com
赞一下 回复此楼
博学、审问、慎思、明辨、笃学
9楼2010-04-28 18:11:51
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 ，六级已过，有科研经历 +7 曦熙兮 2026-03-15 7/350 2026-03-16 16:34 by houyaoxu
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
信息提示
关闭
请填处理意见
关闭