24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3760  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 乌洛托品 的主题更新
信息提示
请填处理意见