| 查看: 2692 | 回复: 6 | |||
gaven0728新虫 (小有名气)
|
[求助]
【求助】castep计算弹性常数和热力学性质问题 已有1人参与
|
 我参考网上教程用castep计算了BN的弹性常数和Ge的热力学性质，发现原胞和晶胞（单胞）的结果有差距，不明白怎么回事。（1）计算BN的弹性常数时用了两种方法，第一种是先把数据库导入的BN晶胞转为原胞，然后结构优化，利用结构优化好之后的原胞直接进行弹性常数计算。第二种方法是利用原胞结构优化，优化好之后转为单胞，然后进行弹性常数计算。计算结果如下，请问各位原胞和单胞的计算结果差距大的原因？以及到底应该用原胞算还是用单胞算晶格常数？？ 单胞： id i j Cij (GPa) 1 1 1 713.55526 +/- 42.164 4 4 4 407.07452 +/- 5.991 7 1 2 128.73305 +/- 19.964 原胞： id i j Cij (GPa) 1 1 1 690.59274 +/- 26.342 4 4 4 381.91248 +/- 12.475 7 1 2 108.47395 +/- 10.319 （2）同样用Ge计算了热力学性质，参数一样的情况下，原胞和 单胞的德拜温度一样，但熵、含、比热容等单胞是原胞的4倍（因为Ge的单胞原子数为8，原胞的原子数为2，相差四倍）。请问我到底是用原胞还是单胞的结果？？ 总之，到底用原胞还是单胞来算晶格常数和热力学性质 ？？？？？ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
2楼2020-04-04 22:07:18
3楼2020-04-05 00:06:06
gaven0728
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5310.3
- 红花: 2
- 帖子: 144
- 在线: 106.1小时
- 虫号: 972400
- 注册: 2010-03-16
- 专业: 金属非晶态、准晶和纳米晶
4楼2020-04-05 09:43:01
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3561.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
5楼2020-04-06 06:31:14
gaven0728
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5310.3
- 红花: 2
- 帖子: 144
- 在线: 106.1小时
- 虫号: 972400
- 注册: 2010-03-16
- 专业: 金属非晶态、准晶和纳米晶
|
卡卡老师，感谢您的回答。我刚接触第一性原理计算，还有两个问题想请教下。 第一个是利用castep算phonons时，有些晶体结构的某个位置可能被多个原子占据，这个时候计算phonons时会提示disordered crystal，对于这种情况如何计算phonons？ 另外一个问题是，在计算phonos提示以下错误，请问如何改参数？ Warning in parameters_validate: current value of ELEC_EIGENVALUE_TOL = 0.100000E-08eV is too large to achieve desired level of convergence of response properties. This may cause convergence failures and/or inaccuracy of results of PHONON calculations - recommend you use a smaller value, e.g. ELEC_EIGENVALUE_TOL = 0.416667E-10eV |
6楼2020-04-10 16:34:52
 |
7楼2021-11-08 15:52:27
50