24小时热门版块排行榜    

查看: 4320  |  回复: 7

乌洛托品

金虫 (小有名气)

[求助] 用cp2k怎么算固体材料的band gap 和 band structure?

有人使用知道怎么用cp2k算固体的能带和带隙么,并且在设置参数时一些需要留心的地方。我看了软件自带的test 里面的例子,没找到相关的例子
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunshoutian168

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
乌洛托品: 金币+5, ★★★很有帮助 2013-01-03 14:48:23
chaizhm: 金币+2, 谢谢~ 2013-04-02 08:18:46
可以算带隙, 参考manual里的  DFT 部分的 print, 可以输出 DOS。
但是 band structure 可能不能计算, 因为cp2k是单K点的,没有多K点功能
想算band structure的话可以用 CPMD, VASP 等等 很多软件都可以

» 本帖已获得的红花(最新10朵)

Ilovechemistry
2楼2013-01-03 13:38:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

乌洛托品

金虫 (小有名气)

送鲜花一朵
引用回帖:
2楼: Originally posted by sunshoutian168 at 2013-01-03 13:38:16
可以算带隙, 参考manual里的  DFT 部分的 print, 可以输出 DOS。
但是 band structure 可能不能计算, 因为cp2k是单K点的,没有多K点功能
想算band structure的话可以用 CPMD, VASP 等等 很多软件都可以

我印象中用vasp只有算分子才可以用单k点,cp2k用单k点算固体,对固体的结构有什么要求么
3楼2013-01-03 14:42:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunshoutian168

银虫 (小有名气)


chaizhm: 金币+1, 谢谢~ 2013-04-02 08:19:10
要求就是  用大的 supercell  减少 k点的影响
Ilovechemistry
4楼2013-01-07 16:51:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

乌洛托品

金虫 (小有名气)

引用回帖:
4楼: Originally posted by sunshoutian168 at 2013-01-07 16:51:33
要求就是  用大的 supercell  减少 k点的影响

我还想问一下,算晶体的时候potential用哪个啊。我用test 里面的example测试,我选了QS里面NE-hybrid-HSE06-lda.inp作为例子计算,就是把其中的SUBSYS选项的Periodic 从NONE改成了XYZ,以及POISSON的Periodic从NONE改成的XYZ,solver改成了periodic。但是输出文件有一句警告:
*** 11:46:40 WARNING in hfx_types:hfx_create_neighbor_cells Periodic    ***
*** Hartree Fock calculation requested without use of a truncated       ***
*** potential. This may lead to unphysical total energies. Use a        ***
*** truncated  potential to avoid possible problems. hfx_types.F line   ***
*** 1348                                                                ***

还有就是,我输出DOS的时候,他没有把空轨道的能级输出来,要怎么设置输出一些空轨道的能量
5楼2013-01-30 12:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)


★ ★ ★ ★ ★
chaizhm: 金币+5, 谢谢~ 2013-04-02 08:19:34
例子中的HSE06是对气相中的一个Ne原子做计算,采用了MT方法消除了周期性边界条件,此Ne原子与其镜像的作用也消除了,所以不需要truncated potential. 你做固体计算,你的盒子中的Kohn-Sham轨道会与周围的镜像盒子中的轨道作用(Hartree potential),但是很远的原子间的轨道作用很弱,因此要截断,因此要用truncated potential。或者short range类型的potenital.

For example in a HSE06 calculation:
    &SCF
      EPS_SCF 1.0E-6
      MAX_SCF 200
      ADDED_MOS 20    //增加空轨道数量。注意如果是LSD,就要写 20 20分别是up和down的空轨道
      SCF_GUESS ATOMIC
      &SMEAR            // 对金属和半导体要smear轨道的占据
         WINDOW_SIZE  0.1
      &END SMEAR
    &END SCF
    &XC
      &XC_FUNCTIONAL
      &XWPBE
        SCALE_X -0.25
        SCALE_X0 1.0
        OMEGA 0.11
      &END
      &PBE
        SCALE_X 0.0
        SCALE_C 1.0
      &END PBE
      &END XC_FUNCTIONAL
      &HF
        &SCREENING
          EPS_SCHWARZ 1.0E-10
          EPS_SCHWARZ_FORCES 1.0E-5
        &END
        &INTERACTION_POTENTIAL
          POTENTIAL_TYPE SHORTRANGE  //标准HSE06使用shortrange potential,截断系数omega=0.11
          OMEGA 0.11
        &END
        &MEMORY
          MAX_MEMORY  512
        &END
        FRACTION 0.25
      &END HF
    &END XC
    &PRINT
        &PDOS     // 打印density of states
        &END PDOS
        &MO        // 打印轨道占据
           OCCNUMS .TRUE.
        &END MO
    &END PRINT
6楼2013-04-02 03:49:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiazf1406

铜虫 (初入文坛)


月只蓝: 金币+1, 鼓励交流! 2014-03-19 20:28:33
请问cp2k中可以计算出pdos代表什么意思?如何利用cp2k中的东西做态密度图和前线轨道分析?
7楼2014-03-19 11:31:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ll385113236

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by sunshoutian168 at 2013-01-03 13:38:16
可以算带隙, 参考manual里的  DFT 部分的 print, 可以输出 DOS。
但是 band structure 可能不能计算, 因为cp2k是单K点的,没有多K点功能
想算band structure的话可以用 CPMD, VASP 等等 很多软件都可以

2.6版本的cp2k也不能设置K点么?
8楼2014-11-11 19:22:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 乌洛托品 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见