| 查看: 3427 | 回复: 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超软赝势进行计算，结果依然有错误提示： 顺便看了一下结果图，的确这个突变有点太那啥了……  由于是工科背景，对于赝势之类的理论比较空白，虽然在恶补，但是暂时依然杯水车薪。哪位高人可以帮忙指出一下，上面的那些步骤，都有哪些问题么…… |
回复此楼
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