| 查看: 2761 | 回复: 6 | ||||||
| 【奖励】 本帖被评价3次，作者watermall增加金币 2.25 个 | ||||||
| 当前主题已经存档。 | ||||||
[资源]
【原创】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 |
» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
国自然面上复盘~欢迎讨论
已经有15人回复
-
中青基了要发朋友圈吗？
已经有3人回复
2楼2009-08-28 09:47:04
3楼2009-08-28 09:58:13
4楼2009-08-28 12:57:15
5楼2009-08-31 10:14:09
6楼2009-10-04 11:25:29
7楼2009-10-04 15:21:38
投票:
30