| 查看: 2766 | 回复: 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超软赝势进行计算，结果依然有错误提示： 顺便看了一下结果图，的确这个突变有点太那啥了……  由于是工科背景，对于赝势之类的理论比较空白，虽然在恶补，但是暂时依然杯水车薪。哪位高人可以帮忙指出一下，上面的那些步骤，都有哪些问题么…… |
回复此楼» 猜你喜欢
304求调剂
已经有8人回复
-
308求调剂
已经有8人回复
-
362求调剂
已经有6人回复
-
26自然地理学303分求调剂
已经有9人回复
-
312求调剂
已经有15人回复
-
化工求调剂！
已经有15人回复
-
284求调剂
已经有8人回复
-
一志愿哈工大，初试329，求环境科学与工程调剂！
已经有5人回复
-
327求调剂
已经有3人回复
-
一志愿吉大化学327求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- CASTEP对优化后的晶体进行分态密度分析，应该怎么操作
已经有14人回复
- 用特定的面去切割CASTEP建立的晶体
已经有5人回复
- castep计算的晶体结构，bands无带隙而dos里面有，如何认定材料？castep遇到导体?
已经有9人回复
- CASTEP建模方法：虚拟晶体近似(VCA)
已经有498人回复
- MS 中的CASTEP可以计算强关联体系的晶体结构（如VO2)吗？
已经有7人回复
- 【求助】CASTEP 增加压强是晶体结构变化？如何设置？
已经有11人回复
- 【求助】CASTEP中怎样构建晶体聚合物呢？
已经有3人回复
- 【求助】用MS中的castep模块计算晶体
已经有12人回复
- 【求助】CASTEP 分子晶体建结构
已经有10人回复
- 【求助】用CASTEP计算晶体时结构问题
已经有10人回复
- 【求助】castep建立晶体结构时出现孤立的原子是怎么回事
已经有3人回复
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
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
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
5楼2010-04-28 12:54:04
zzy870720z
荣誉版主 (文坛精英)
- 1ST强帖: 3
- 应助: 47 (小学生)
- 贵宾: 9.05
- 金币: 30914.3
- 散金: 5613
- 红花: 68
- 沙发: 99
- 帖子: 12592
- 在线: 23567.6小时
- 虫号: 745488
- 注册: 2009-04-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
6楼2010-04-28 13:03:58
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
 有变化/波动到是不在乎，可是问题在于，这个错误提示：Error in md_vv1_npt_nose - cell inconsistency? Current trace stack: md_vv1_npt_nose md_run castep Error in md_vv1_npt_nose - cell inconsistency? Current trace stack: md_vv1_npt_nose md_run castep MPI Application rank 3 exited before MPI_Finalize() with status 1 这是失败后的castep文件的末尾显示的错误信息……偶就是不明白为什么总会有这么个信息…… 至于那两个图，是因为通过观察，有上述错误信息时，在Motion和Temp的两个图中，至少有一个会发生剧烈变化……不过，这个没关系，可以不care的说……重点还是在于为什么MD的计算会以cell inconsistency的问题失败…… |
7楼2010-04-28 13:20:44
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
咋没人回帖了捏…… 8楼2010-04-28 17:17:55
zzy870720z
荣誉版主 (文坛精英)
- 1ST强帖: 3
- 应助: 47 (小学生)
- 贵宾: 9.05
- 金币: 30914.3
- 散金: 5613
- 红花: 68
- 沙发: 99
- 帖子: 12592
- 在线: 23567.6小时
- 虫号: 745488
- 注册: 2009-04-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
9楼2010-04-28 18:11:51
franch
版主 (著名写手)
- 1ST强帖: 2
- 应助: 62 (初中生)
- 贵宾: 0.185
- 金币: 7766.8
- 散金: 2986
- 红花: 23
- 沙发: 2
- 帖子: 1964
- 在线: 731.1小时
- 虫号: 587928
- 注册: 2008-08-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
10楼2011-04-03 19:38:08
