24小时热门版块排行榜

返回列表

【奖励】本帖被评价3次，作者watermall增加金币 2.25 个

当前主题已经存档。

watermall

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1981
帖子: 286
在线: 295.6小时
虫号: 118944

[资源] 【原创】pwscf 电子和声子能带计算中的K点的生成

这是我的方法来设定pwscf中计算能带所需要的K点
1. 对结构进行优化，然后用xcrysden看结构，看特殊k点
2. 结合Materials studio 中特殊K点的设置
3. 写syml 文件（下面详述syml文件的写法）
4. 运行附件中的bk.exe，得到所需要的K点。注意将最后一个K点去掉，（这是一个小bug，这是一个多余的K点，然后把K点数减1就行了）
5. 用pw.x计算能带的时候，注意，kpoint的单位为 “crystal”，不能用默认的tpiba单位
6. 在进行声子谱的计算的时候，matdyn.x 只认识tpiba单位，怎么办呢？这里我提出一个办法，非常简单，就是用pw.x 进行能带结构的计算，没必要让能带计算进行完，大概5秒钟就能强制停止，然后看out文件，会发现pw.x会将crystal 单位的K点转换为 tpiba单位的k点。手动删掉些东西，就能得到matdyn.x所需要的以tpiba为单位的K点。

syml文件的写法：
它的输入文件为syml，输出文件为inp.kpt。其中syml输入文件的格式如下：
8
15 15 15 15 15 15 15
G 0.0 0.0 0.0
K -0.33333333333 0.6666666666667 0.000000000
H -0.33333333333 0.6666666666667 0.500000000
A 0.0 0.0 0.5
G 0.0 0.0 0.0
M 0.0 0.5 0.0
L 0.0 0.5 0.5
A 0.0 0.0 0.5
第一行用来标记有多少个特殊k点，下面是这些特殊k点之间每个要分多少个k点，接着就是这些特殊k点的坐标。
产生的inp.kpt可以之间拷贝到pw.x在计算能带时的输入文件中。

[ Last edited by watermall on 2009-8-28 at 09:56 ]

回复此楼

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

第一性原理计算辅助工具

vasp

QE-pwscf

» 猜你喜欢

1楼 2009-08-28 08:59:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★★★★★ 五星级,优秀推荐

非常好
顶一个；）

但有个问题，要请教一下；）
1：您是用两个特殊k点的坐标相减，然后将其等分的方法得到k点么？
2：xcrysden和MS给出的特殊k点的名称一致么？
谢谢；P

[ Last edited by xirainbow on 2009-8-28 at 09:53 ]

赞一下

回复此楼

2楼2009-08-28 09:47:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

watermall

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1981
帖子: 286
在线: 295.6小时
虫号: 118944

★
freshgirl(金币+1,VIP+0):谢谢热心分享经验 9-2 11:59

引用回帖:

Originally posted by xirainbow at 2009-8-28 09:47:
非常好
顶一个；）

但有个问题，要请教一下；）
1：您是用两个特殊k点的坐标相减，然后将其等分的方法得到k点么？
2：xcrysden和MS给出的特殊k点的名称一致么？
谢谢；P

[ Last edited by xirainbow o ...

1，是的
2，建议你用用xcrysden和MS试试看。一般xcrysden是不给出特殊K点的名称的，所以要通过MS来确定名称

赞一下(10人)

回复此楼

3楼2009-08-28 09:58:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖