24小时热门版块排行榜    

查看: 3293  |  回复: 13
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yanhy

木虫 (正式写手)

[交流] 【求助】castep中同一晶体为什么make p1 后计算结果差别很大 已有3人参与

我前几天练习使用castep,计算Ag的态密度,金属Ag的对称群为fm-3m(编号225),我先按照这个结构计算了一下,然后make p1, 再进行同样的计算,两次得到的态密度图很不一样,计算结果详见附件,哪位高手帮忙指点一下这是为什么,应该如何处理,此外,如果make p1 对计算结果影响这么大的话,那么掺杂计算计算的结果的可信度会不会受到影响,因为进行原子取代时往往要进行make p1 操作。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
freshgirl(金币+1,VIP+0):谢谢参与~ 9-2 12:35
两次的DOS其实是一样的,只是MAKE P1后,对称性低了,高能量的那部分没有计算(矩阵变大了,高能量的那些本证值没有计算出来)。
7楼2009-08-21 16:32:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

xirainbow

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★
freshgirl(金币+1,VIP+0):谢谢你的认真回答~ 8-19 22:29
yanhy(金币+6,VIP+0):十分感谢 8-19 22:32
首先,你必须把两次的K点设置告诉大家
其次,由于对称性的差异,在相同的K点设置前提下,不同的对称性会给出不同的K点(具体的K点数值MS的输入文件中有)。这同样能导致DOS的差别。
再次,如果你所用 的K点并不收敛,给出上面的图形就并不奇怪了;)

所以,建议你首先检查两次计算使用的真正的K点是否一致(是真正的K点,不是你设置的K点)
然后进行K点的收敛测试,应该就不会出现这样的问题了

[ Last edited by xirainbow on 2009-8-19 at 22:11 ]
2楼2009-08-19 22:09:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyb4142

金虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuli8(金币+3,VIP+0):3ks 8-20 09:58
对称性不同的两个体系,怎么比较?
1. 默认情况下,不同对称性K点取得不同,DOS和能带图自然不一样,特别是态密度图原子多了,密度就会大。
2. 虽然都是Ag,当你改变对称性时,一个事我们常见的Ag金属,另一个就不知何物了,显然晶胞不同其(物质不同),自然结果不可比较了
结论:若果你要比较掺杂体系,自然要保证体系不变(或者说对称性不变)一个可行的方法就是都makeP1,在相同规模下用超晶胞来计算,保持体系一致性
3楼2009-08-20 09:18:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guifan

木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuli8(金币+2,VIP+0):感谢讨论 8-20 17:20
make P1将使结构失去对称性,这样的话相应的参数就要进行相应的修改,只要参数选择得当,一般对计算速度影响较大,对计算结果影响不大
4楼2009-08-20 15:54:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:03 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 00:51 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[基金申请] 国社科系统bug了,是不是要放榜了? +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
信息提示
请填处理意见