| 查看: 2869 | 回复: 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 是对的 ( 有跟文献比较 ) 我的原子位置也是对的 (完全符合该空间群的对称性) |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有21人回复
-
垃圾破二本职称评审标准
已经有15人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有20人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- pwscf 自洽计算，电子步难收敛的问题
已经有12人回复
- 请做过Partical charge density的前辈进来帮忙看看输入文件有何问题？
已经有8人回复
- WARNING:stress and force are not correct
已经有6人回复
- vasp计算弹性时的错误（使用defvector方法）
已经有25人回复
- 计算好像出错了，CONTCAR和WAVECAR空空如也！
已经有10人回复
- vasp 计算 partial charge 时k点选取的问题
已经有19人回复
- 求助VASP能带计算的中断原因
已经有6人回复
- 各位大侠，vasp下计算完能带如何作图
已经有12人回复
- 关于赝势与自洽计算
已经有15人回复
- hse计算中ICHARG=0，不收敛
已经有10人回复
- VASP结构优化时出现的POTIM问题
已经有4人回复
- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- vasp5.2计算半导体Si光学性质的总结
已经有106人回复
- 流化床曳力系数UDF
已经有9人回复
- vasp计算在进行静态计算时，出错。
已经有3人回复
- 【求助】band structrue 出错
已经有13人回复
- 【讨论】自洽计算与非自洽计算的含义、所做的工作
已经有10人回复
- 【求助】vasp5.2光学性质问题
已经有18人回复
- 【求助】Hse06计算时KPOINTS出错
已经有30人回复
- 【分享】用VASP进行Partial Charge分析实例
已经有22人回复
skydark
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 44 (小学生)
- 金币: 3652.2
- 散金: 30
- 红花: 6
- 帖子: 302
- 在线: 249.8小时
- 虫号: 1927082
- 注册: 2012-08-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
大神 我想这回我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的那一步后 得出的波函数都不一样 请问这会是什麽原因所造成的呢?  未命名3.png  未命名G.png  未命名x.png  未命名M.png  未命名y.png  未命名y2.png  未命名y_conduction_band.png  未命名neary.png |
6楼2015-05-06 01:16:27
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
franch: 金币+3, 谢谢回帖交流。。。。 2015-03-31 09:12:52
skydark: 金币+10, ★★★很有帮助 2015-04-17 11:31:55
感谢参与，应助指数 +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
skydark
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 44 (小学生)
- 金币: 3652.2
- 散金: 30
- 红花: 6
- 帖子: 302
- 在线: 249.8小时
- 虫号: 1927082
- 注册: 2012-08-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
大神 我把所有平面波数据加总了 并且把real space 的 wavefunction画出来了 但图有点怪怪的 我画的是gamma point 那点的 最高价带那条能带 的波函数 如图所示 在一个unit cell 内取200点 X点座标 画出波函数在X轴上的空间分?� 左上角的是波函数的实数部份 右上角的是波函数的虚数部份 下方是波函数的绝对值平方 虚数部分的波函数看起来挺正常的 但实数部分的波函数很奇怪 不符合週期性边界条件 且第一个峰值去到 +0.04 ,第2个峰值却去到 -0.05 不知道这是什麽原因所造成的? 注: 我读取WAVECAR档时 是用wavetrans 的 WaveTransSpinor.f90  圖片13.png |
3楼2015-04-17 12:15:04
4楼2015-04-17 15:35:37