24小时热门版块排行榜    

查看: 3252  |  回复: 12

阿飞2012

铜虫 (初入文坛)

[求助] 用castep进行高压下铁的第一性原理计算,设置和结果对比的问题?

我初学castep,请问:1,自旋极化对能量和结构有什么影响。2,为什么我设置K-point为16. 16 16时老是失败?最多只能是15 15 15    .3,在一篇文献上看到别人计算出来的体积弹性模量为200左右,而我计算出来是350左右,这是为何?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shmilylxc

木虫 (正式写手)

【答案】应助回帖

★
感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2012-05-12 10:12:24
MS算出来的体积模量误差比较大,你要是算的话可以先算弹性常数
2楼2012-05-11 21:50:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿飞2012

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by shmilylxc at 2012-05-11 21:50:02:
MS算出来的体积模量误差比较大,你要是算的话可以先算弹性常数

非常感谢
3楼2012-05-11 23:41:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿飞2012

铜虫 (初入文坛)

还有没有人做这方面的研究,欢迎参加讨论!我做铁的相变,结果在12GPA就相变了,理论上是15GPa左右,这又是为何?请高手们发表下高见!!
4楼2012-05-14 08:07:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿飞2012

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by 阿飞2012 at 2012-05-14 08:07:10:
还有没有人做这方面的研究,欢迎参加讨论!我做铁的相变,结果在12GPA就相变了,理论上是15GPa左右,这又是为何?请高手们发表下高见!!

第二个问题我解决了一部分,就是16 16 16 的k-point 和400的cut-off energy 可以计算,550就不行,但不清楚这是为什么。用14 14 14 的k-point和550的cut-off 反而能较好的得到相变压强。
5楼2012-05-14 08:12:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qx0813

木虫 (正式写手)

【答案】应助回帖

★
感谢参与,应助指数 +1
WDD880227: 金币+1, 多谢交流提示 2012-05-14 13:55:48
引用回帖:
5楼: Originally posted by 阿飞2012 at 2012-05-14 08:12:46:
第二个问题我解决了一部分,就是16 16 16 的k-point 和400的cut-off energy 可以计算,550就不行,但不清楚这是为什么。用14 14 14 的k-point和550的cut-off 反而能较好的得到相变压强。

这可能跟你选择的精度关系很大,k点越多,对计算机的配置要求会更好,但是单独提高K点对计算精度影响的改善可能意义不大,关键是协调各个参数,共同达到最佳效果。另外,你所说的相变是两种不同结构之间的转化吗?
遇到最好的自己
6楼2012-05-14 10:38:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljyang320

木虫 (正式写手)

小书虫爱吃糖

★
WDD880227: 金币+1, 多谢提示 2012-05-14 13:56:07
引用回帖:
4楼: Originally posted by 阿飞2012 at 2012-05-14 08:07:10:
还有没有人做这方面的研究,欢迎参加讨论!我做铁的相变,结果在12GPA就相变了,理论上是15GPa左右,这又是为何?请高手们发表下高见!!

我觉得你在计算之前应该做收敛性测试
不然你的计算都是没有任何用处的。
但使岁月静好,现世安稳!
7楼2012-05-14 10:54:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿飞2012

铜虫 (初入文坛)

引用回帖:
7楼: Originally posted by ljyang320 at 2012-05-14 10:54:06:
我觉得你在计算之前应该做收敛性测试
不然你的计算都是没有任何用处的。

我看过一些文献,基本都是用1 6 16 16 的k-point 和 400的cut-off,这样也不行吗
8楼2012-05-14 13:45:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿飞2012

铜虫 (初入文坛)

引用回帖:
6楼: Originally posted by qx0813 at 2012-05-14 10:38:36:
这可能跟你选择的精度关系很大,k点越多,对计算机的配置要求会更好,但是单独提高K点对计算精度影响的改善可能意义不大,关键是协调各个参数,共同达到最佳效果。另外,你所说的相变是两种不同结构之间的转化吗?

恩,我已经看到了Fe在15GPa左右实现了bcc结构向hcp结构的相变,不过这都是别人做过的了,我只是验证了一下而已。我想做出新东西来,比如计算压强对弹性系数的影响,目前正在学习中。有问题的话希望不吝赐教,那顺便就问问你对这个研究有什么看法吗?或者意见,非常感谢
9楼2012-05-14 13:48:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qx0813

木虫 (正式写手)

【答案】应助回帖

★ ★
fzx2008: 金币+2, 谢谢指教 2012-05-14 20:12:47
如果别人研究过了,估计在研究的意义不大,不知道你是不是用的心方法研究的,如果是的话,还是有一定意义的。
遇到最好的自己
10楼2012-05-14 15:48:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿飞2012 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 SX5mJmPWDfOO 2026-11-13 4/200 2026-11-15 07:34 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:33 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 9/450 2026-11-15 02:04 by xSrQNSSjcgoX
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-15 00:35 by Equinoxhua
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1Nb0IVlmV6v0 2026-11-11 4/200 2026-11-13 11:18 by Equinoxhua
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:11 by Equinoxhua
[考博] 27博士申请(市政/环境) +3 ladf 2026-11-12 3/150 2026-11-13 08:19 by Equinoxhua
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 11:22 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:14 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 3/150 2026-11-12 08:26 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 06:03 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
信息提示
请填处理意见