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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » gk.f 的输入文件syml如何准备
71/1返回列表
查看: 2797  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wangjiang12

铁虫 (小有名气)

[求助] gk.f 的输入文件syml如何准备 已有3人参与

各位高手好，

我准备用vasp计算的体系是六方层状结构，但是有缺陷的，如何准备syml文件啊，求问syml文件的格式。非常感谢
回复此楼

» 收录本帖的淘帖专辑推荐

精华网帖收集 第一性原理相关文档

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-04-29 15:29:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
六方结构可以到MS找到相应的对称性，加入syml即可。这样的对称性是最全的，若是按照文献给出的，一般采用相同的即省去的原则，可以采用如下syml.

6   
30  20  30  20  30
G 0.0    0.0    0.0
M 0.0    0.5    0.0
K -0.333333333333   0.6666666666667   0.0
G 0.0    0.0    0.0
A 0.0    0.0    0.5
L 0.0    0.5    0.5
0.000000000 -4.052756964  0.000000000     0.142458646 -0.246745613  0.000000000
3.509790486  2.026378482  0.000000000     0.284917292  0.000000000  0.000000000
0.000000000  0.000000000  5.572964774     0.000000000  0.000000000  0.179437703
-20.0  20.0
0.0

上面每一行意思：
#############comment
1 line: nhighk   !how many high symmetry k-points
2     : ndiv(i)  !number of grid points between two symmetry k-points
3     : labhk(i), phighk(i,j)   !label of high sysmmetry k-points, coordinates
......
9     : a(i,j),b(i,j)   !direct & reciprocal lattice vectors
......
12    : emin, emax   ! energy width of band plot
13    : efermi       ! fermi level
##############
\\\ G--M--K--G--A---L   NKPT=131
\\\ 1--31--51--81--101---131
--------------------------------------------
Hex
M:0.0    0.5    0.0
G:0.0    0.0    0.0
H:-1/3   2/3    0.5
K:-1/3   2/3    0.0
L:0.0    0.5    0.5
A:0.0    0.0    0.5
赞一下(1人) 回复此楼
2楼2014-04-29 18:31:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

古泉月

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
楼主解决了吗？我也在学习ＶＡＳＰ用的侯博士的教程，也是遇到了这问题，还没解决
赞一下 回复此楼
3楼2014-04-30 10:38:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)
引用回帖:
2楼: Originally posted by mywai520 at 2014-04-29 18:31:45
六方结构可以到MS找到相应的对称性，加入syml即可。这样的对称性是最全的，若是按照文献给出的，一般采用相同的即省去的原则，可以采用如下syml.

6   
30  20  30  20  30
G 0.0    0.0    0.0
M 0.0    0.5  ...

请问费米能级那一栏怎么修改？
比如我从静态计算OUTCAR grep Fermi 为4.5326
0.0 是改为4.5236 还是 -4.5326 ？
赞一下 回复此楼
相互学习，共同进步。
4楼2014-05-13 15:13:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苦行客

木虫 (正式写手)

【答案】应助回帖

你grep Fermi 的到的是4.5326，那就把0.0改为4.5326。
赞一下 回复此楼
行有不得，反求诸己～
5楼2014-05-13 15:56:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)
引用回帖:
4楼: Originally posted by 对马拉二胡 at 2014-05-13 15:13:39
请问费米能级那一栏怎么修改？
比如我从静态计算OUTCAR grep Fermi 为4.5326
0.0 是改为4.5236 还是 -4.5326 ？...

得到多少就改多少。你的是4.5326
赞一下 回复此楼
6楼2014-05-13 18:24:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)
引用回帖:
6楼: Originally posted by mywai520 at 2014-05-13 18:24:52
得到多少就改多少。你的是4.5326...

大神，还想请教一个问题。为什么我算出来的能带数据BANDS文件都要手动减去费米能级的值才能得到正确的图？（步骤都是按手册来的，gk.x生成KPOINTS文件，pbnd.x生成BANDS文件）
赞一下 回复此楼
相互学习，共同进步。
7楼2014-07-17 17:50:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangjiang12 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭