| 查看: 3288 | 回复: 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超软赝势进行计算，结果依然有错误提示： 顺便看了一下结果图，的确这个突变有点太那啥了……  由于是工科背景，对于赝势之类的理论比较空白，虽然在恶补，但是暂时依然杯水车薪。哪位高人可以帮忙指出一下，上面的那些步骤，都有哪些问题么…… |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
didi5158
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 482.4
- 帖子: 250
- 在线: 89.9小时
- 虫号: 273449
- 注册: 2006-08-20
- 性别: GG
- 专业: 无机材料化学
12楼2011-06-10 14:18:28
zzy870720z
荣誉版主 (文坛精英)
- 1ST强帖: 3
- 应助: 47 (小学生)
- 贵宾: 9.05
- 金币: 30914.3
- 散金: 5613
- 红花: 68
- 沙发: 99
- 帖子: 12592
- 在线: 23567.6小时
- 虫号: 745488
- 注册: 2009-04-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
2楼2010-04-28 12:23:31
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
3楼2010-04-28 12:33:17
zzy870720z
荣誉版主 (文坛精英)
- 1ST强帖: 3
- 应助: 47 (小学生)
- 贵宾: 9.05
- 金币: 30914.3
- 散金: 5613
- 红花: 68
- 沙发: 99
- 帖子: 12592
- 在线: 23567.6小时
- 虫号: 745488
- 注册: 2009-04-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
4楼2010-04-28 12:37:26
40