24小时热门版块排行榜    

查看: 3498  |  回复: 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一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 4/200 2026-08-29 22:27 by ZPa0EcMwuECS
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-29 21:35 by 4qO8zhfjr2AX
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-29 21:32 by 4qO8zhfjr2AX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:11 by 4qO8zhfjr2AX
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 20:07 by winsaint
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 17:30 by 4FFAWE8HcgUD
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 13:51 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:29 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
信息提示
请填处理意见