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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【原创】pwscf 电子和声子能带计算中的K点的生成
71/1返回列表
查看: 2753  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价3次，作者watermall增加金币 2.25 个
当前主题已经存档。

watermall

金虫 (小有名气)

[资源] 【原创】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

金虫 (小有名气)
★
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的回帖

gavinliu7390

木虫 (著名写手)

★ 一星级,一般

在xcrysden 生成K点后。   在写个脚本， 把k点提出来，再晶体坐标转直接坐标就ok了。 不用那么麻烦。
赞一下 回复此楼
4楼2009-08-28 12:57:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuzhaonan

铁杆木虫 (著名写手)

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

楼主提供的方法还是比较正规的，但楼上的方法简单些
赞一下 回复此楼
5楼2009-08-31 10:14:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guizi~


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

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

[ Last edited by xirainbow o ...

您说的"用两个特殊k点的坐标相减，然后将其等分的方法得到k点"是什么意思呀?有点不懂了
我是新来的,请多指教
   谢谢!
赞一下 回复此楼
6楼2009-10-04 11:25:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
引用回帖:
Originally posted by guizi~ at 2009-10-4 11:25:


  您说的"用两个特殊k点的坐标相减，然后将其等分的方法得到k点"是什么意思呀?有点不懂了
我是新来的,请多指教
   谢谢!

就是平面几何中：已知两个特殊点，求他们连线的方法
赞一下 回复此楼
7楼2009-10-04 15:21:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 watermall 的主题更新
71/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 05:13 by winsaint
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中，共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 基金未中，这种答复是模板吗？ +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
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
关闭
请填处理意见
关闭