24小时热门版块排行榜    

查看: 3509  |  回复: 13

skydark

木虫 (正式写手)

[求助] WAVECAR波函数的疑问 已有2人参与

我最近读取了 WAVECAR里的 波函数数据出来
但却觉得数据跟物理直觉差太远

假设 同一条band 的两点K点  (K1+G=K2) 且 (K1 = -K2) ,
那波函数 of K1点 应该与波函数 of K2点 相等或只差个负号,
但WAVECAR给出的波函数数据 却 既不相等 也非只是差了个负号

如下
   0.0000000      0.50000000       0.0000000   <-- K点座标
           1 <--第几条band           5468
(   -13.1314     ,    0.00000     )    1.00000
     0     0     0  (  -0.161791     ,   0.579942E-01 )  <-- 平面波 与其 CG 系数 (实数与虚数)
     1     0     0  (   0.134249     ,  -0.303834E-01 )
     2     0     0  (  -0.457166E-02 ,   0.163872E-02 )
     3     0     0  (   0.835663E-03 ,  -0.189125E-03 )

   0.0000000     -0.50000000       0.0000000
           1        5468
(   -13.1314     ,    0.00000     )    1.00000
     0     0     0  (   0.754764E-01 ,  -0.482369     )
     1     0     0  (   0.529809E-01 ,   0.283576E-01 )
     2     0     0  (   0.213271E-02 ,  -0.136301E-01 )
     3     0     0  (   0.329789E-03 ,   0.176500E-03 )

请问一下各位高手  为何会出这样的错误呢?
该如何修正?

注: 我确定我的band structure 是对的 ( 有跟文献比较 )   我的原子位置也是对的 (完全符合该空间群的对称性)
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+4, 谢谢指导 2015-04-17 21:32:41
skydark: 金币+5, 有帮助 2015-04-18 18:18:16
??? PAW POTENTIAL ???P?S
???~?????????????????? NORM ?????
??????? 1 (?????????? 1)
???????? 0.9 (?????? 0.8)
???????~?@??????? PAW????
PAW ?o??????Ч??~????Z?????Щ???_??
4楼2015-04-17 15:35:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
franch: 金币+3, 谢谢回帖交流。。。。 2015-03-31 09:12:52
skydark: 金币+10, ★★★很有帮助 2015-04-17 11:31:55
不是错误,差的不一定是正负号,而是一个相因子。因为平面波基组的基组在每个k点是不同的,形式为exp(i (G+k) r),最后的相因子同时和G/k/晶格都有关系,VASP中基组的排列顺序是按G固定的,换了k点以后,对应的基组在排列中会变掉,所以不能直接对比前几个。
你可以把所有波函数的模全部列出来(大概有几千个,需要计算一下),就能看出来对应关系了。
2楼2015-03-27 21:56:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

skydark

木虫 (正式写手)

???????:
2?: Originally posted by KalaShayminS at 2015-03-27 21:56:20
?????????????????????????????????????????沨?????????????k?????????????exp(i (G+k) r)????????????????G/k/?????й????VASP?л??????????????G??????????k???????????? ...

????  ?????????沨?????????
?????real space ?? wavefunction????????
????е????

???????gamma point ???? ???????????? ???????

??????
?????unit cell ???200?? X??????
????????????X??????????
??????????????????????  ???????????????????????
?·?????????????????

???????????????????????????

????????????????????
???????L??????????
???????????? +0.04  ,??2????????? -0.05

??????????????????????

?: ????WAVECAR???
????wavetrans ?? WaveTransSpinor.f90
WAVECAR????????????
?D?13.png

3楼2015-04-17 12:15:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+2, 谢谢回复 2015-04-17 21:32:59
skydark: 金币+5, 有帮助 2015-04-18 18:18:21
引用回帖:
3楼: Originally posted by skydark at 2015-04-17 12:15:04
大神  我把所有平面波数据加总了
并且把real space 的 wavefunction画出来了
但图有点怪怪的

我画的是gamma point 那点的 最高价带那条能带 的波函数

如图所示
在一个unit cell 内取200点 X点座标
画出 ...

这个我只能怀疑实空间求和有问题
因为所有基组都是符合周期性的 不管求和用的什么数据 都不可能破坏周期性
5楼2015-04-17 16:38:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skydark

木虫 (正式写手)

???????:
5?: Originally posted by KalaShayminS at 2015-04-17 16:38:44
????????????????????????
??????л??鶼???????????? ???????????????? ?????????????????...

????

?????????wave function?????????????????
?????????????band ?Щ???????wave function???????????

??????????????????????? (unitcell?????????? ???4???????)
band structure ?? ???λ????????

????highest valence band ?? Y???????????????---??μ????????????????ò????

??????????????????????????(????????)

?????highest valence band ??gamma point ??wavefunction (??X?????)

?????highest valence band ??x point ??wavefunction (??X?????)

?4??highest valence band ??M point ??wavefunction (??X?????)

?5??highest valence band ??Y point ??wavefunction (??X?????)

?6??highest valence band ??Y point ??wavefunction (??X?????)(??2?μ????????)

?????lowest conduction band ??Y point ??wavefunction (??X?????)

?????highest valence band ??"near"Y point ??wavefunction (??X?????)

?????? X??G??M?? ??????near Y?????????????????
"exactly"??Y point??? lowest conduction band ??wavefunction??????

???"exactly"??Y?????highest valence band?????????????
???????????CHGCAR ????????VASP plot band????????? ?ó???????????????


????????????????????????
WAVECAR????????????-1
δ????3.png


WAVECAR????????????-2
δ????G.png


WAVECAR????????????-3
δ????x.png


WAVECAR????????????-4
δ????M.png


WAVECAR????????????-5
δ????y.png


WAVECAR????????????-6
δ????y2.png


WAVECAR????????????-7
δ????y_conduction_band.png


WAVECAR????????????-8
δ????neary.png

6楼2015-05-06 01:16:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

引用回帖:
6楼: Originally posted by skydark at 2015-05-06 01:16:27
大神

我想这回我wave function的傅利叶转换有做对了
但不知为何在某条band 某些高对称点上wave function却会全然是错误的

我同样还是模拟单层黑磷这种简单结构 (unitcell是正长方形 只有4颗磷原子)
band s ...

不好意思,这个过于细节了,我也看不出这是什么问题
7楼2015-05-06 10:47:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skydark

木虫 (正式写手)

引用回帖:
7楼: Originally posted by KalaShayminS at 2015-05-06 10:47:22
不好意思,这个过于细节了,我也看不出这是什么问题...

沒關係  

我還是非常謝謝你 多次幫我解答問題
8楼2015-05-06 11:50:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goodluck1982

铁虫 (初入文坛)

【答案】应助回帖

引用回帖:
6楼: Originally posted by skydark at 2015-05-06 01:16:27
大神

我想这回我wave function的傅利叶转换有做对了
但不知为何在某条band 某些高对称点上wave function却会全然是错误的

我同样还是模拟单层黑磷这种简单结构 (unitcell是正长方形 只有4颗磷原子)
band s ...

你画的波函数是分实部和虚部分开画的还是画的模方?
如果是分实部虚部画的,不同应该是正常的,因为不同次计算得出的波函数可以整体上差一个相因子,这就会导致实部虚部的不同,但是模方应该相同。你有没有比较过不同计算的模方是否一样?
9楼2015-08-18 10:37:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zuocuiping

木虫 (职业作家)

你的wavecar用什么脚本读的呢
谢谢
10楼2015-08-19 22:23:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 skydark 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 5/250 2026-08-31 02:00 by hZiFeudZyoGR
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 6/300 2026-08-31 01:48 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:15 by hZiFeudZyoGR
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 00:56 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-31 00:45 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 ASdOkHsho7FD 2026-08-28 12/600 2026-08-31 00:34 by hZiFeudZyoGR
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-31 00:32 by hZiFeudZyoGR
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 12/600 2026-08-31 00:31 by hZiFeudZyoGR
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 19:01 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:40 by hZiFeudZyoGR
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见