应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】单胞与原包能量不一样？
191/2返回列表12下一页
查看: 3550  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

superest

金虫 (小有名气)

[交流] 【求助】单胞与原包能量不一样？ 已有7人参与

我在做一个体系的计算，单胞得到的能量是 -191.12291424 eV，而2X2X1supercell得到的能量是-761.878750eV（-191.12291424 X4=764.49165696）。两者计算条件相同，单胞K点442，supercell的K点222. 结构确定没有写错。supercell的结构是用Ms从单胞扩出来的。请问是什么原因导致两者之间4eV的能量差呢？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-16 21:51:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

stractor

金虫 (著名写手)
★ ★ ★ ★
mazuju028(金币+4):说得很详细。 2010-10-17 18:54:23
超胞能量怎么能理解成原胞的能量的简单叠加呢？！
使用超胞后，晶格尺寸改变，影响总能。例如一维情况，假设原胞的晶格常数为a,使用原胞的时候1.5a处的原子是周期性镜像原子，0.5a和1.5a的原子相同，而构成超胞后，1.5a处的原子和0.5a处的原子就不互为镜像原子。我不好表达我的意思，就是说构成超薄后，布洛赫波会改变，从而能量不是简单原胞能量的叠加。

[ Last edited by stractor on 2010-10-17 at 11:28 ]
赞一下(4人) 回复此楼
10楼2010-10-17 11:21:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
mazuju028(金币+2):谢谢提示 2010-10-17 18:53:26
建议楼主把OUTCAR里面的总能量的各项都给出来，就是有EDWARD的那些。看看是那一项的改变不到4倍
赞一下(3人) 回复此楼
8楼2010-10-17 10:30:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zxzj05

荣誉版主 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bingmou(金币+1):感谢参与 2010-10-17 10:27:39
mazuju028(金币+2):谢谢分享。 2010-10-17 18:52:05
mazuju028:谢谢交流 2010-10-17 18:52:25
请搞清楚单胞和原胞的定义。

你的2*2*1不是原胞，是超晶胞。
你所给的K点不同，足以影响计算结果了。

[ Last edited by zxzj05 on 2010-10-16 at 22:23 ]
赞一下(4人) 回复此楼
储氢家族欢迎储氢研究者！
2楼2010-10-16 22:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+1):谢谢回复 2010-10-18 14:22:48
呵呵，不知道你算得是不是LiFePO4，感觉能量很接近

[ Last edited by sunyang1988 on 2010-10-16 at 22:55 ]
赞一下(2人) 回复此楼
心之所向无惧无悔
3楼2010-10-16 22:54:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by zxzj05 at 2010-10-16 22:17:59:
请搞清楚单胞和原胞的定义。

你的2*2*1不是原胞，是超晶胞。
你所给的K点不同，足以影响计算结果了。

[ Last edited by zxzj05 on 2010-10-16 at 22:23 ]

就是说相同的K点密度并不能保证计算结果相近了？
赞一下(1人) 回复此楼
心之所向无惧无悔
4楼2010-10-16 22:58:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superest

金虫 (小有名气)
引用回帖:
Originally posted by zxzj05 at 2010-10-16 22:17:59:
请搞清楚单胞和原胞的定义。

你的2*2*1不是原胞，是超晶胞。
你所给的K点不同，足以影响计算结果了。

[ Last edited by zxzj05 on 2010-10-16 at 22:23 ]

2*2*1的supercell测试了K点2×2×2，3×3×2，4×4×4，能量均相等（误差在1meV）.30个原子的cell，测试了K点4×4×2，5×5×2和6×6×6，能量误差在0.1meV左右，应该不是K点密度问题
赞一下 回复此楼
5楼2010-10-17 09:33:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superest

金虫 (小有名气)
引用回帖:
Originally posted by sunyang1988 at 2010-10-16 22:58:30:

就是说相同的K点密度并不能保证计算结果相近了？

2*2*1的supercell测试了K点2×2×2，3×3×2，4×4×4，能量均相等（误差在1meV）.30个原子的cell，测试了K点4×4×2，5×5×2和6×6×6，能量误差在0.1meV左右，应该不是K点密度问题
赞一下 回复此楼
6楼2010-10-17 09:33:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
恩，这样看来超胞能量就是不一样。另外超胞的带隙与原胞是不同的
赞一下(1人) 回复此楼
心之所向无惧无悔
7楼2010-10-17 09:48:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
引用回帖:
Originally posted by sunyang1988 at 2010-10-17 09:48:15:
恩，这样看来超胞能量就是不一样。另外超胞的带隙与原胞是不同的

带隙会不同？举个例子吧
赞一下 回复此楼
9楼2010-10-17 10:31:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 superest 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 纯娱乐，不喜欢勿喷 +4 Tide man 2026-08-04 7/350 2026-08-04 11:23 by Tide man
[基金申请] 娱乐 +3 Tide man 2026-08-03 3/150 2026-08-04 09:58 by gltch
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 jRl3mE6ddZGq 2026-08-03 10/500 2026-08-04 09:26 by v1TJC4qJauAB
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化，科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭